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PHYSICS OF THE ALIVE 


Synergetics and Phase Transitions: Mounting the Quantum Ladder of Nature 


A.V.CHALY, I.S.DOBRONRAVOVA, S.P.SIT'KO 


Scientific Research Center “Vidhuk”, Kiev, Ukraine 
Ukrainian state Medical University, Kiev, Ukraine 


The problem of relationships between two approaches to living systems — synergetics and 
quantum mechanics — is discussed. Special attention is paid to the common methodological 
aspects of these approaches based on the fundamental ideas of the modern phase transition 
theory. 


(Received August 29, 1994) 


The problem of entire stability of macroscopic systems is one of the crucial points in natural 
science, being a subject of unremitting attention (see, for example, [1—9]). Stability of solids as an 
evidence, was starting point in constructing of classical mechanics, both in justifying its 
measurements (rigid rods) and in establishing the applicability limits of its theoretical principles 
(reference systems). However, due to the high sensitivity of dynamical systems to variability of the 
initial conditions, classical physics by itself was not able to substantiate the observed stability of 
solids. That circumstance was recognized by the creators of quantum mechanics. At the same time, 
the possibility of explaining (within the quantum mechanics framework) stability of solids, at least 
on their elements’ level, was discovered. 

Weisskopf [10] was probably the first who noticed that most of the authors of quantum 
mechanics books treating the methodological problems concentrated mainly upon the probabilistic 
character of quantum-mechanical results. That could be justified at the stage of forming the 
probabilistic way of scientific thinking to comprehend the micro-cosm processes. However, 
emphasizing the statistical nature of the quantum-mechanical laws often grew into the idea of 
uncertainty of the quantum mechanics results, in contrast, for example, to the exact ones of 
classical mechanics. This is incorrect. As a matter of fact, just the objects which the quantum- 
mechanical principles of identity and discreteness can be applied to are the truly integrable systems 
whereas the conventional mechanics, as a rule, describes ergodic systems whose infinitesimal 
uncertainty of their initial conditions results in homogeneous filling the phase space (regime of the 
strange attractor type). This circumstance has been completely cleared up only in the course of 
modern revolution in natural science when applying to physics the results of elaboration of the 
nonlinear dynamics mathematical basis. However, in this new context the problem of solids 
stability looks somewhat different from what could be considered before. It is important to 
emphasize that in this case it could be preserved the meaning of idea that namely the applicability 
of quantum-mechanical principles to such microobjects as nuclei, atoms and molecules determines 
their stable existence in the capacity of subjects of corresponding sections of physics as a 
fundamental science. It is these objects that form the three steps of the Weisskopf “quantum 
ladder” [10]. In terms of physics of the alive [9], the fourth step of this ladder is the alive. 

As a product of nonlinear thinking, the modern evolutionary outlook shifts the stress from 
the statement of the stable systems existence to the problem of their formation and dynamical self- 
sustention. In the passage of such systems through the critical (bifurcation) points, owing to an 
anomalously high sensivity to the external factors (even of extremely weak intensity) at this instant. 
the fluctuation effects give rise to a new ordered phase, often having qualitatively new properties. 


Physics of the Alive, Vol. 2(1), © 1994. “Vidhuk” . З 


A.V.CHALY, I.5.DOBRONRAVOYVA, S.P.SIT'KO 


It is especially important to consider the quantum systems which form the quantum ladder 
steps from the viewpoint of self-organization. The necessity of such treatment of the alive is 
obvious. Considering self-organization of quantum systems of nuclei, atoms and molecules offers 
basic scope for uniting physics as a science on a new nonlinear basis. For the time being, there 
exist two parallel worldviews in physics, namely the quantum-relativistic and nonlinear (or 
synergetical) ones. Such a condition is maintained by the linearity of quantum-mechanical 
equations and the fact that the dynamical chaos propagation described within nonlinear dynamics is 
confined by the applicability limits of quantum mechanics [11]. However, this situation cannot be 
accepted as satisfactory one, since the conceptions of reality within the abovementioned worldviews 
are quite different and cannot consistently coexist. Besides, dynamical chaos does not exhaust all 
the diversity of nonlinear phenomena. Among these, as the most important ones can be considered 
self-organization and self-reproduction of relatively stable structures, hierarchy of dissipative 
structures in living organisms in particular. Thus, initiated by physics of the alive, studying the 
relationships between synergetics and quantum mechanics is of great importance for the whole 
physics. ў 

The present paper is aimed to a considerable degree at the methodological problems of 
physics, physics of the alive as well. Our main purpose is to show the connection between the 
modern theory of phase transitions and i) the theory of formation of ordered structures 
(synergetics) and ii) the quantum theory of dynamically stable systems beyond the self-organization 
threshold. 

In spite of the fact that the connection between nonequilibrium phase transitions and 
dissipative structures is a subject of numerous investigations [3, 5], we believe that there are new 
ideas and results based upon applying the fluctuation models to the processes of self-organization. 
In the first section of the paper we show how the fluctuation models of self-organization processes 
and formation of stable ordered structures can be constructed in accordance with the modern 
theory of phase transitions. We briefly discuss the main results of studying of these models. 

In the second section we examine an analogy between the renormalization group method, 
and quantum mechanics. The analogy originates from solving virtually the same. Sturm-Liouville 
problem of seeking eigenfunctions and eigenvalues [12]. In that way the theory of phase transitions 
forms a basis on which it is possible to consider in a new fashion the problem of formation of 
stable ordered structures in its connection with the basic principles of quantum mechanics. 


SYNERGETICS AND PHASE TRANSITIONS: FLUCTUATION MODELS OF SELF- 
ORGANIZATION PROCESSES 


In papers [13—16] and monograph [17] new models of the self-organization processes were 
proposed and studied. The main idea of the models was to coordinate the initial principles of 
model construction and the corresponding equations of motion (kinetic equations) with the 
fundamentals of modern physics of phase transitions. Let us note that the known kinetic models of 
self-organization (such as "brussellator", *oregonator", Girer-Meinhardt's model, etc.) do not 
satisfy this requirement, since they cannot be obtained by making use of the Landau-Ginsburg 
Hamiltonian for several interacting order parameters, whereas the new models can. The structure 
of the fluctuation models obtained in [13—17] is directly connected with the following form of the 
free energy fluctuation part: 


AF uc = >. Ньз[@ + Hoo] (1) 


Here 


1 1 
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is the Landau-Ginsburg Hamiltonian for the order parameter ф(х); а, bp c; are coefficients, given 
that aj ~ T — Te (T, is critical temperature for the i-th order parameter) d is the space 
dimensionality, and A, is the Hamiltonian (free energy) of interaction between two (several) order 
parameters. 
According to (1-2), the motion equations have the form of the known kinetic equations of 
Landau-Ginsburg's time-dependent theory: i 
др; oH, 
E 3 ан * bigi + Я + Д; до; (3) 
One can easily see that motion equations (3) coincide in their form with the kinetic 
equations of reaction-diffusion models widely used in synergetics [5]. Indeed, the nonlinear terms 
(those containing g and Hj) assure the feedback (interaction) among the order parameters, 
whereas the nonlocal (Ornstein-Zernike) contribution (containing Ло) is responsible for the 
diffusion processes (L;C; = D, are diffusion coefficients). Moreover, beyond the phase transition 
threshold, when self-organization causes the formation of dynamically stable structures, D « 0 [17], 
and, therefore, Eq.(3) does not differ in its form from motion equations in standard synergetic 
models 


4 = -kq TA kaks. (3a) 
associated (at k, < 0 and k, > 0) with the Landau-Haken potential and having the limit-cycle 
solutions in phase plane [9] (one should bear in mind that in general the interaction Hamiltonian 
can be explicitly taken into account during the standard procedure of proceeding from (1) to (3) by 
using the methods of calculus of variation. 

For definiteness, let us consider a fluctuation model with two interacting order parameters 
and nonlinear coupling of the general form 


E = Hop) + ВА, " 
2 = Яд, да) + Р.А), 


where ‘¥(g,, p) and G(gi, p) are arbitrary nonlinear functions. Analyzing the type of singular 
points according to Poincare's classification and the stability of stationary solutions 9,,(x, f) with 
respect to small perturbations 


õp; (x,t) = ул) - eu Xt) = | Z pilk, Ae" at (5) 
k 
implies examining the characteristic equation 
Det[Ü - Ôk? - | = 0 (6) 
which gives the dispersion law and relates the damping coefficient 4 to the wave vector k. Here 
[e cone 
^ |e 90р) ae, 1, 
ds eu 0 
p, др, 


and 7 is 2 x 2 unit matrix. 
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In [14—17] a series of statements (theorems) determining the conditions of the emergence of 
Hopfs and Turing's bifurcations was formulated for fluctuation model (4). For example, the 
necessary condition of the emergence of uniform-in-volume temporal oscillations reads. as the 
following inequality for the coefficients of dispersion equation (6): 


(E, - E) * AE, <0, 
where 
Dk? - ди др, 
E, = Dk? - 2/ ap), 
E, = eY ap ‹ 22/69, . 
In the specific case E, = Е, = 0 the periodic limit-cycle solution of model (4) exists for 


Е; < 0 only. 
In [18] a numerical modelling of ordered structures was made for the fluctuation model with 


Im 
D 


quadruple interaction, i.e., with Hi, = 44 9; (x) o2 (x) ах. : 


Special attention should be paid to the fact (confirmed by both analytical and numerical 
7 calculations) that the ordering and self-organization within the framework of the considered models 
take place beyond the boundary of system's stability, that is on the nonthermodynamic branch of 


states, where the determinant of inverse susceptibilities Det Ü and/or the determinant of diffusion 


coefficients det D change their signs. 

Concluding this Section, we should note that in the vicinity of bifurcation points the 
diffusion coefficients D, (in the general case — all kinetic coefficients) and the inverse relaxation 
times c;! = Га, ~ y;! (related to system's susceptibilities y;'s) included in motion equation (3) 
cease to be the local functions of spatial coordinates and time. Consistent taking into account this 
basic fact needs using the results of papers [19—21] in which the following long-range (in the 
vicinity of critical, or bifurcation, points) effects were simultaneously taken into account: i) spatial 
nonlocality and “memory” in the equation of state which are related to the value of the system 
anomalous susceptibility y (7,7) to the external factors; ii) spatio-temporal dispersion of the kinetic 
coefficients, such as diffusion coefficients D,’s. 

Applying the results of papers [19—21] to examining the fluctuation models of self- 
organization processes leads to the following equations linearized for quantities бік, f) from (5): 

dsp; (k, t ^ l Г 
T Y ue My (köp (К, 1) X № (Kk, - t") óp; (Kk, t") ar (7) 
j j 
Here the complete matrix of inverse relaxation times 


M,(k,t- гу = MPO- r) + №, n) 


is decomposed into the following matrices: a) M Ni awhich corresponds to hydrodynamical 
equations describing only the nonlocality of the equation of state in the Ornstein-Zernike 
approximation; b) N(k,t - t") which describes the effects of gonlocality and memory in the kinetic 
coefficients. 


In view of the similarity of Eqs.(7) and the time-dependent Schrödinger equation Й is 
convenient to introduce the new variables 


gi (kt) = pue бр (Kk, t), 
y 


sagi 
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which are’ the litéar combinatiotis of Fourier transforms др, 4's. This'procedure is equivalent to 
proceeding to the tritéraction ‘representation in the ‘quantuny’field theory and-enables ‘us’ to-use the 
technique ‘of е perturbation’ tHéory to éaleülate 'the’ quantities "pK; ‘ава consequently “the 
corresponding pair correlation futictions б», DAR Ww а ! 

As noted in [22] when considering the theoretical basis of morphogenesis; 'such попібсаї 
effects can initiate the break of symmetry and instabilities of Hopfs and Turing’s types. 02 
PHASE TRANSITIONS AND QUANTUM MECHANICS 


The latest advance in physics of phase transitions and critical phenomena has resulted from 
employing the concepts and methods of scaling theory [23], with the renormalization group 
approach being its macroscopic justification [24]. The renormalization group (RG) transformation 
is complex nonlinear transformation, usually applied to the “seed”, or “effective” Landau- 
Ginsburg Hamiltonian Ніс (see Eq.(1)). The coefficients а» b, c; subjected to the RG- 
transformation are described as the point и in so-called parametric space, that is ш = (ai вусі). 


Within the RG-method framework the pair correlation function of the order parameter 
G(r, и) = (ot) - (009) 100) - (oc D 


r=lh~A, (8) 


(in terms of which, by virtue of the known relations of the condensed matter statistical theory most 
of the equilibrium and kinetic properties of substance can be expressed) satisfies the following 
equation 


G;(k, п) = 92 10, (sk, Ё, u) (9) 


Here R, is the RG-transformation operator, s is a certain great parameter which is usually related 
to the fluctuation correlation radius, 7 is the critical index of anomalous dimensionality of the 
correlation function G». If the point д of parametric space is close to the fixed (immovable) point 
4* which is invariant under the RG-transformation (R,utu*) and plays the role of a critical 
(bifurcation) point, then the following relation takes place: 
Жи=и*+У tse, (10) 
i21 

where e; and 5” are eigenfunctions and eigenvalues of the RG-transformation operator Ri 
linearized in the vicinity of the fixed point, and г; is the coefficient of the linear expansion of the 
quantity ди = д — y* into a series in eigenfunctions ег. 

Let us note that the power laws of scaling theory (scaling laws) are determined by the 
greatest positive eigenvalue у, of the operator RE , given that y; = 1/v, where v is the critical index 
of temperature dependence of the radius of the density fluctuation correlation, namely: 


2 = Sor °. (11) 


Here до is the correlation radius amplitude which coincides in the order of magnitude with the 
intermolecular interaction radius, г = (T — То)/ Tc — is deviation of temperature T from its value 
Те in the critical point, and v is the critical index which equals 0.63 for d — 3 in the systems with a 
scalar (one-component) order parameter. Other eigenvalues y(/ > 2) specify different scaling 
corrections (see review [25]), making it possible to extend considerably the applicability limits of 
the phase transition theory with respect to the thermodynamic variables (temperature, density, 
pressure, etc.). 
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From all the aforementioned it definitely follows that the methods and mathematical tools 
employed in the modern fluctuation theory of phase transitions, on the one hand, and in quantum 
mechanics, on the other hand, are very similar. Indeed, both the problem of finding critical indices 
like the index v and that of finding the quantum numbers are in essence Sturm-Liouville's 
problems. Moreover, this analogy (see Table) seems to be even closer in view of not only the 
common mathematical tools but also the basic symmetry reasons. The very fact of the existence of 
quantum numbers and stationary states has appeared as a manifestation of the fundamental 
properties of space and time symmetry as well as the processes in the microcosm. The existence of 
critical indices and universal commonality of properties pertaining to the systems of completely 
different. nature displays a new class of symmetry, namely that defined by the renormalization 
group transformations. This new type of symmetry turns out to be inherent in all, not only micro- 
but also macroscopic, systems in the vicinity of their critical (bifurcation) points. A new phase, new 
ordered structures included, emerges in such systems due to fluctuation intcraction which can be 
taken into account with use of the RG-method only. 


Table 
Correspondénce between quantum mechanics and physics of phase transitions (synergetics) . 


Phase Transitions (Synergetics) 


y" 
^ 


Linearized RG-operator К/ 


Eigenfunctions V, | Eigenfunctions of operator RL -@; 
Eigenvalues Е, Eigenvalues of operator R- s" 
Quantum numbers Critical indices 


| 
| 


Thus, based upon the concepts of scaling and renormalization group, the contemporary 
development of the phase transition theory makes it possible 
i) to evaluate in a new fashion the validity and profundity of the known principle “via fluctuation 
to ordering”; 
ii) to give a new interpretation of the steps of the “Weisskopf quantum ladder”. 


СИНЕРГЕТИКА TA ФАЗОВІ ПЕРЕХОДИ — ПІДЙОМ ПО «КВАНТОВИХ СХОДАХ" 


О.В.ЧАЛИЙ, I.CJIOBPOHPABOBA, С.П.СІТЬКО 


Аналізується проблема співвідношень двох підходів до живих систем -- синергетики та 
квантової механіки. Особливу увагу приділено загальним методологічним аспектам у цих 


підходах, що грунтуються на фундаментальних ідеях сучасної теорії фазових переходів. 
СИНЕРГЕТИКА И ФАЗОВЫЕ ПЕРЕХОДЫ — ВОСХОЖДЕНИЕ ПО “КВАНТОВОЙ · 
ЛЕСТНИЦЕ” 

А.В.ЧАЛЫЙ, И.С.ДОБРСНРАВОВА, С.П.СИТЬКО 


Обсуждена проблема соотношений между двумя подходами к живым системам — 
синергетикой и квантовой механикой. Особое внимание уделено общим методологическим 
аспектам в этих подходах, основанных на фундаментальных идеях современной теории 
фазовых переходов. 
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Corrective Action of Millimeter Wave upon Biological Systems of Different 
Hierarchical Levels | 


L.S.BUNDYUK, A.P.KUZMENKO, N.N.RYABCHENKO, G.S.LITVINOV, S.A.BOBROVNIK, 
L.I.MATSEIKO, I.V.GORBENKO, O.R.FEDOROVICH : 


Scientific Research Center "Vidhuk", Kiev, Ukraine 


The data on normalizing action of low intensity millimeter wave (MMW) radiation upon isolated 
cells, damaged by ionizing irradiation, are represented. It is shown that MMW applied to 
organism's biologically active points exert a corrective and immunomodulating influence upon 
functional states of animals and patients with different pathologies: MMW effect has a 
pronounced resonant character. The decrease of mm-radiation power density values approaching 
107^ W/cm? increases the therapeutic effect ої MMW. 


The mechanisms of functional state recovery of biological systems of different hierarchical levels 
under MMW action are suggested. 


( Received June 6, 1994) 
INTRODUCTION. 


Interest to low intensity millimeter wave (MMW) electromagnetic radiation interaction with 
biological objects has been recently intensified. The achieved theoretical and experimental results 
allowe1 investigators to lay the groundwork for resolving the fundamental problems of physics of 
the alive, which culminated in the creation of a new trend in medicine Milliwave Resonance 
Medicine (MRM) [1—3]. MRM takes into consideration the fact that each organism posesses a 
coherent electromagnetic field and that in pathology this field is disturbed. This disturbance can be 
corrected with an external electromagnetic field of intensity approaching quantum threshold [4]. 
The method consists in controlled exposure of biologically active zones of a human body to 
extremely low intensity electromagnetic action with individually chosen parameters. 

The therapeutic effect of MMW in general could be explained by the fact that energy of 
quanta of these waves passing through the skin in the area of biologically active zone is absorbed by 
the definite biological microparticles, and this absorbtion has a bioresonant character. As a result 
the bioinformational signals are generated which regulate the organism's physiological functions as 
a whole. The frequency and intensity dependence of MMW action on bioobjects are very important 
for realization of the MMW therapeutic effect. 

Unlike the so-called "thermal" effects, the response of biological systems to MMW radiation 
may show a sharp resonance frequency dependence with relative width of action spectra less than 
107—107". Threshold dependence оп the intensity and exposition time is specific to resonance 
effects [1, 3, 5, 6]. 

It has been established that MMW at the range from 40 to 140 GHz may affect a wide range 
of properties of different biological objects: alter proliferation, differentiation, content and 
cytochemical indices of blood cells [7, 8], change the pathogenicity of bacterial cells [9], stimulate 
the immune system [10—12], normalize cells' viability damaged by ionizing irradiation [13]. 

Taking into account the increasing levels of ionizing irradiation in the enviroment which 
produce the damaging influence the ability of MMW radiation to correct the organism's functional 
state has become very important. Until now the mechanism of MMW action upon bioobjects has 
not been explained, therefore, new experiments which can help to acquire such knowledge are 
quite appropiate. 

The data on MMW corrective action on isolated cells, damaged by ionizing irradiation, as 
well as on experimental animals and patients with different pathologies, are represented in this 
paper. In order to find the frequency and intensity dependencies of MMW effect, different sources 
of irradiation were used. The attempt to clear up the mechanisms of functional state recovery of 
biological systems of different hierarchical levels under MMW action was made. 
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CORRECIIVE ACTION OF MILLIMETER WAVE UPON BIOLOGICAL SYSTEMS OF DIFFERENT HIERARCHICAL LEVELS 


A. COMBINED EFFECT OF MILLIMETER WAVE AND GAMMA IRRADIATION UPON 
L-LINE CELLS 


'The influence of MMW radiation upon the viability of L-line cells exposed to different doses 
of y-quanta "Cs was examined. The frequency dependence of the combined effect was studied. 


MATERIALS AND METHODS 


L-line cells (mice transformed fibroblasts [14]) were grown as a monolayer in the Eagle 
tissue culture nutrient with glutamine and 10% foetus calf serum. The cells grew on the coverslides 
placed on the bottom of penicillin flasks in the thermostate at 37°C. Into every flask 2 ml of the 
culture medium were added, which contained 2x10* tumor cells. Two days later, when L-cells 
entered the stationary phase of growth, the cultures (except for controls) were exposed to ;-quanta 
(radiological unit with the irradiation source "Св, rate of dose — 0.19 Gy/min). The irradiation 
doses were: 0.25, 0.5, 1.0, 6.0 and 10.0 Gy. Low dose effects were provided by 0.25 and 0.5 Gy. 
The use of 1.0 Gy dose resulted in death of 10% of cells (LD), 6.0 Gy — LDs, 10.0 Gy — 10%. 


The first set of experiments 


Immediately after the exposition to yquanta 50% of the cells’ cultures were exposed to 
ММУ radiation. Noise generator "Porog-1" (output power density — 10 ^—10 ^ * W/cm/, frequency 
range level of 54—76 GHz at 0.5 of power) was used as a source of radiation [15]. Exposition time 
was 7 min. Before the exposition the cultural medium was removed from the flasks, and the emitter 
was set inside the flask directly 2 cm above the cells' layer. After the exposition fresh cultural 
medium was added, the flasks were incubated at 37°C. Experimental cultures were divided into the 
following groups: 1) controls, intact cultures (without impact), 15 samples; 2) exposed to MMW 
action, 15 samples; 3) exposed to different doses of ;-quanta, 15 samples per each of 5 doses; 4) 
exposed to different doses of y-quanta and then to MMW, 15 samples per each dose of y-quanta; 
5) exposed to MMW with the following action by »-quanta (the dose was 6.0 Gy), 15 samples. 


The second set of experiments 


Forty eight hours after the cultivation, when L-cells entered the stationary growth phase, the 
flasks were divided into eight groups, 20 flasks in each: the Ist group — controls (the cells were not . 
exposed); the 2nd group — L-cells were exposed to y-quanta in the dose of 6.0 Gy, the 3rd, 4th 
and Sth groups — L-cells were exposed to MMW radiation in the frequency range from 35.9 to 
42.3 GHz, from 42.3 to 48.9 GHz and from 48.9 to 55.1 GHz, respectively; 6th, 7th, 8th groups — 
the cells were subjected to the combined effect of y-quanta and MMW in the frequency range from 
35.9 to 42.3 GHz, from 42.3 to 48.9 GHz, and from 48.9 to 55.1 GHz, respectively. 

P2- 68 Scanning generator was used as a MMW radiation source (with output power 
1 mW/cm! in the range of 35.9 to 55.1 GHz). The time of exposure was 7 min. 

The flasks with cultures were placed on a special horn and MMW exposure was carried out 
through the*flasks" bottoms. At the level of the horn radiation power was equal to 1 mW/cm’, the 
MMW radiation acted on cellular monolayer at flux power of 0.5 to 0.6 mW/cm’ (taking into 
account absorption by glass). 

Twenty four hours after the y-quanta, MMW or combined exposure the coverslides with cells 
were taken out of flasks, fixed by alcohol, painted with hematoxylin-eosin, lucidated with xylene, 
and placed in balsam. In micropreparations the number of cells per 0.1 mm’ (cell layer density) 
and mitotic index — MI (the number of dividing cells per 1000 cells, %0) were determined. 

In the first and second sets of experiments a total of 540 and 480 cultures were used, 
respectively. In statistic processing of quantitative data Student's criterion was used. 
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RESULTS AND DISCUSSION 
The data of the first set of experiments showed that MMW action had no influence upon L- 
line cell growth (certain variations were not reliable and remained within the limits of the 
calculation error), and the exposition by ;-quanta (6.0 Gy dose) resulted in almost double 
reduction both in cellular density, and mitotic index (Table 1). 
Table 1 
The influence of y-quanta (dose 6.0 Gy) and MMW radiation combined effect upon the growth of L-line 
cultures 
Mitotic index, %o 


Type of action 
Days of growth Days of growth 


| 
: + 6 3 : 


MMW + y-quanta | 11+2.5 20x2.1* | 20-22" | 5422.3" 60+ 1.8* 
y-quanta + MMW 


Note: * — variations are Statistically reliable compared to control (p < 0.05); 
** — variations are statistically reliable compared to cultures exposed to y-quanta only (p < 0.05). 


In cases when MMW was applied after у-диапіа irradiation, normalization of the growth in 
L-line cell cultures (p < 0.05) was observed: cellular density and their MI were almost identical to 
those of controls. If MMW radiation was employed prior to y-irradiation, cellular density and MI 
were almost similar to the corresponding data in the irradiated cultures, i.e., the protective action 
of MMW in this combination was not observed. Apparently it was levelled by the subsequent great 
dose of ionizing radiation (6.0 Gy, LD5). 

In Table 2 the data on MMW effect upon the density of cells and their MI changes in 5-day 
L-line cultures with preexposition by various doses of y-quanta are shown. Cellular density was : 
almost identical in controls, in specimens exposed to small y-quanta doses (0.25— 1.0 Gy), and in 
specimens under MMW action. Mitotic index in the y-exposed cultures was a little higher than in 
control and under MMW action, which complied with the data found in literature: small 
irradiation doses enhance the activation of cell proliferation, the large ones inhibit it [7, 9]. It was 
confirmed by the exposure of the culture to 10 Gy dose of у-диапіа (LDso, see Table 2). 


Table 2 
Combined effect of y-quanta and MMW radiation on the data of growth of 5 day L-line culture 


Type of action Cellular density 
* 


34223 102102 
38+2.0 133=11.5 
y-quanta, 0.25 Gy 33+4.6 160+12.1 


10.0 Gy + MMW = 
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Observation of giant cells in the micropreparations showed rather high content of giant cells 
(13+1.1%о) in 5-day culture of L-line cells, which pointed to disturbance of mitotic cell division. 
As is known, giant cells are formed as a result of partial division: nuclei divide, but cytoplasm 
remains undivided. Such cells are not viable, they eliminate from the cultures, i.e., perish.. 
Exposition by "Porog-1" showed a double reduction of their number (61.3960, p < 0.05). Ionizing 
irradiation increased the giant cells number, i.e., resulted in the disturbance of mitotic division of 
these cells, and the higher radiation dose, the more of giant cells were observed. MMW action 
normalized this ratio bringing it closer to the control data. It is especially noticeable under the 
exposition to large doses of y-quanta. Thus, at the dose of 10.0 Gy giant cell ratio has grown to 
3041.1%0 from 13+1.1%о in control (Fig.1). The subsequent MMW action upon this da. 
reduced the giant cell ratio to 122: 1.7960, that was close to the control level. 

Thus, in the L-line cell cultures 
exposed to ionizing irradiation MMW 
radiation of  nonthermal intensity 
normalizes the process of mitotic cell 
division reducing giant cell ratio, which 
stimulates cell culture viability, i.e., 
produces a normalization effect in case 
of radiation injury of cells. 

Study of the slides obtained in the 
second set of experiments showed that 
the dynamic of growth in the control 
series was as follows: the cells’ density 
on the 4th day of growth equalled 
45+2.7 cells, on the 5th day — 60+3.0 
and on the 6th day — 57+3.2; and MI was 
88+3.2%0, 128+3.7%0 and 77+3.6%o, 
respectively. In the cultures exposed to 0,25 0.5 
quanta in a dose of 6.0 Gy the density _ : : j 
of cells was lower than in the control Fig.1. The number of giant cells in 5-day L-line culture under MMW 

and y-quanta irradiation, %o: 


series: 23+2.3, 41 222 and 20+3.1, 4 - in the culture exposed to y-quanta with subsequent MMW action; 
respectively to the day, and MI was @ — in the culture exposed to different doses of у-диапіа; 
51+3.1%0, 73+3.5%о and 56+3.2%0 on ~~ — level, which determines the number of giant cells in the control culture 


Ist, 2nd and 3rd day after exposure. As is obvious from these data, y-quanta in dose of 6.0 Gy 
almost twofold reduced density of cells and from 1.5 to 2-fold — MI of L-line cells. 

Table 3 
L-cells growth on n the 3rd day after exposure to MMW radiation і in ади Шла aic Ai ан 


Cellular density 
BEY SUE NOM £ POE EM Days of E 
сна єння 


we 
> 
Т 


t2 
cn 
zT 


— > 
м © 
T 


= 
> 


Number of giant cells, %o 


CA 


1.0 6.0 10.0 


Dose, Gy 


MMW 35.9 = 42.3 GHz | 50220 59223 13829.37 RE 10829.17 
MMW 42.3 — 48.9 GHz 3 Sui 


Note: * — р < 0. 05 as eS to T 


As seen in Table 3, MMW in the frequency range of 35. 9 to 42.3 GHz almost did not affect 
cells’ density as compared to the control cultures, though MI of L-cells was somewhat higher than 


15 


LS .BUNDYUK, A:P-KUZMENKO, N.N. RYABCHENKO et al.  \ у 


that of the control series. After exposure to MMW. in other frequency. ranges. cells': density 
reduction was noticed alongside with their unreliable MI increase. à 

In cases of.combined action, і:е.; when the cell cultures were: exposed to y-quanta and then 
exposed to MMW, a normalizing MMW effect was noticed: the indices under study approached 
those of the control set and were reliably higher than those of the cultures exposed to y-quanta 
(p € 0.05). These data are given in Table 4. As is shown in Table 4, the best indices of cells' density 
and their MI normalization (that is, the closest to the control values) were noticed in the cultures 
exposed to a combination of ;-quanta and MMW exposure in the freguency range of 35.9 to 
42.3 GHz. 

Table 4 


Combined effect of y-irradiation and MMW on 3 days old L-cells' cultures 


Type of exposure Cellular densit 
Days of growth Days of growth 


4 
Control 
dua i OC 
y-irradiation + MMW. 
35.9 — 42.3 GHz 42+2.8** 


y-irradiation + MMW. 
42.3 — 48.9 GHz 41+2.5** 
y-irradiation + MMW 

3343.3 


48.9 — 55.1 GHz 


Note: * — p < 0.05 as compared to control; ** — р < 0.05 as compared to cultures damaged by ;-irradiation. 


B. MMW CORRECTIVE ACTION UPON DISTURBED FUNCTIONAL STATES IN 
EXPERIMENTAL ANIMALS 


B.1. ANTITUMOR EFFECT OF MMW 

The impact of MMW radiation upon tumor process with regional and distant metastatic 
spread by exposition of the acupuncture points in animals was studied. The assessment of the 
possible mechanisms of the impact was suggested. 


MATERIALS AND METHODS 


Investigations were carried out on 73 4-week-old mongrel rats weighing 150—200 g each and 
68 6-week-old mice of C57BI line weighing 22-28 g each bred in R.E.Kavetsky Institute for 
Oncology and Radiobiology Problems of the National Ukrainian Academy of Sciences. Two strains 
of inoculated tumors were studied: Guerin's carcinoma in rats and Lewis lung carcinoma in mice 
with metastatic spread into regional lymphatic nodes and lungs, respectively. Transplantation of 
Lewis and Guerin's carcinoma cells was carried out by subcutaneous injection of 2x10* and 10° 
tumor cells per animal in the right shin surface. Industrial generator with output power of 
3 mW/cm' served as a source of MMW radiation. Impact was performed upon the point located іп 
the same anatomic region as Zusanli in humans 24 hours after the inoculation. Choice of the point 
was conditioned by its profound impact upon both humans and animals. Locus of the point on the 
left extremity was exposed to radiation, thus preventing the possibility of direct MMW action upon 
tumor. Position of the points was checked by means of a device for electropuncture diagnosis and 
therapy "Electronica ELITA-4M". Animals were divided into 4 groups according to frequency 
characteristics of the exposition: Г group — 35.9—42.3 GHz, II group — 42.3—48.9 GHz, Ш 
group 48.9—55.1 GHz, IV group — controls (inoculated animals, non-exposed to MMW). In given 
frequency ranges the impact was carried out in sweeping regime at 12.5 Hz frequency. The course 
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of treatment included 10 daily sessions 3—7 min each, depending on the magnitude of skin electric 
resistivity in the point. In the first set of experiments animals were decapitated under ether 
narcotization on the 16th day (Guerin's carcinoma) and on the 23d day (Lewis lung carcinoma) 
after the inoculation. Antitumor effect was estimated by inhibition of primary tumor growth by the 
mass of regional lymphatic nodes, number and size of metastases in lungs. Energy status of the 
animals was estimated by the amount of adenine nucleotides in liver tissue by means of thin-layer 
chromatography on siluphol plates [16]. Cytochemical investigations of lactate dehydrogenase 
(LDG, 1.1.1.27) and succinate dehydrogenase (SDG, 1.3.99.1) in peripheral blood lymphocytes 
were conducted by the method described in [17]. Biochemical examinations were performed before 
and immediately after the course of MMW. Intact group — 8 rats and 10 mice was of the same 
line as the animals with implanted tumors. In the second set of experiments the percentage of 
complete recovery was detected as well as the percentage of Guerin's carcinoma metastatic spread 
into regional lymphatic nodes. Average life time (ALT) of rats was determined. Animal observations 
had been conducted for 4 months since MMW exposure. Statistical data processing was performed 
with the use of Student-Fisher and X criteria. 


RESULTS AND DISCUSSION 


It was established that the most substantial inhibition of tumor growth in the animals with 
Guerin's carcinoma was noted in group I (Table 5). Maximum primary tumor regression and 
minimal lymphatic node metastatic disturbance were observed in the same group. After 10 MMW 
sessions inhibition of Lewis lung carcinoma growth in mice from groups 2 and 3 made up 15,8 and 
35.5%, respectively, and in group I stimulation of growth by 27.3% were observed (Table 6). 
Metastatic spread in lungs was almost the same both in experimental groups and controls, though 
the volume and the number of metastases in group I and 3 were essentially smaller than those of 
controls (p < 0,05). It should be noted that the application of MMW radiation at 35.9—42.3 GHz 
frequency resulted in metastatic inhibition, though stimulated the growth of primary tumor, while 
the use of 42.3—48.9 GHz frequency produced the reverse effect. It makes possible to assume that 

. the amount of therapeutic effect depends on impact frequency characteristics. And the ranges of 
frequencies that have a therapeutic effect upon primary tumor and metastatic spread may vary. 


Table 5 


Effect of acupuncture point exposure to mm-range electromagnetic radiation on tumor process in rats 
with Guerin’s carcinoma 
Groups of animals 


Tumor process data 
2 


n 2 a 
Nuwberoteniina о | 9; Е o a | 
5.2+1.8* 12.5+3.4 9.1=2.4* 16.9 1.8 


4 


Tumor growth inhibition, % 


Mass of regional lymphatic 
nodes, mg 112.4+18.0* | 158.3+19.6 | 119.1+16.9* | 239.0+24.6 


Lo Seaan — | 3 | 8$. |. 9. o 7 

Frequency of regional pM 
metastatic spread, % 0 222^ 42.5 
Frequency of tumor complete B n d 
regression, 96 55.5* 22.0* 


40.0+4.2 37.0x 1.5 30.0x 3.5 35.8x1.9 


Note: * — р < 0.05 as compared to group 4. 


ALT of animals, days 
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Table 6 


Effect of acupuncture point exposure to mm-range electromagnetic radiation on tumor process in animals 
with Lewis lung carcinoma 


у а 


Groups of animals 


em ELLE 
Number of animals 15 14 
Primary tumor mass, g 2.3+0.4** 1.5=0.1* 1.2+0.1* 1.8+0.3 
Inhibition (-) or stimulation не 
+) of tumor growth, % oro 215899 -35.5* 

18.4+2.9** 10.5+1.9* 19.6+1.8 
29.1+8.9** TIREI. 7AT 15.7+6.9* 


Note: * — p < 0.05 as compared to group 4; ** — p < 0.05 as compared to group 3. 


Frequency of metastatic 
spread, % __ 

Average number of metastases 
per mouse 

Average volume of metastases 
per mouse, mm: 


100 


10.8=2.3* 


To estimate the possible mechanisms of discovered antitumor effect of MMW action some 
data of the animals’ energy metabolism were studied: ATP and AMP content in liver tissue, SDG 
and LDG activities in peripheral blood lymphocytes. Considering that lymphocytes are classified as 
cells with mostly aerobic metabolism and their major demand for energy is provided by respiration 
[18], correlation between respiration enzyme (SDG) and glycolysis enzyme (LDG) is of primary 
importance. These data are shown in Table 7. Thus, when the tumor is growing, the shift in 
lymphocytes towards glycolysis takes place and it should affect their functional activity. After 10 
MMW sessions partial rehabilitation of SDG activity was observed. The amount of MMW radiation 
correcting influence on lymphocyte energy metabolism in peripheral blood depended on frequency 
characteristics. 


Table 7 
SDG/LDG activity ratio in blood lymphocytes in tumor implanted animals after MMW action 


SDG/LDG 


Number of Guerin's Number of Lewis lung 


animals carcinoma animals і carcinoma 


1.62+0.14* 10 
0.96+0.08* ~ Ù 
0.760.225" 

0.34+0.19%* 
Intact 


vet m" m 


Note: * — p < 0.05 as compared to group 4; ** — p < 0.05 as compared to intact animals. 
p p 


No. 


In different groups of animals the analysis of adenine nucleotides content in liver tissue 
showed that the presence of tumor process led to an essential decrease in ATP content (Table 8). 
Such a decrease pointed to the deterioration of organism's energy status, probably, due to the 
disturbance of ATP synthesis processes in mitochondria. It was also confirmed by the decrease of 
SDG activity in peripheral blood lymphocytes. Besides, the growth of malignant tumors was 
responsible for authentic increase in AMP content in liver tissue what could point to an activation of 
ATP decomposition process. In tissues of the animals with Guerin's carcinoma from experimental 
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group 1 the highest ATP and the lowest AMP contents as compared to those of control group were 
detected immediately after the course of MMW. Among mice with Lewis lung carcinoma the 
highest ATP and the lowest AMP contents were observed in group 3. From the above data it was 
evident that adenine nucleotide content in group 1 of rats and in group 3 of mice approached the 
same data in intact animals which may be indicative of energy metabolism rehabilitation in these 
groups of animals due to MMW radiation effect. In the other experimental groups corrective 
MMW radiation action on energy metabolism was less pronounced. 


Table 8 
Changes in adenine nucleotide content (micromole/g of tissue) in tumor implanted animals' liver after 
MMW action 


Number 
of 
animals 


2.35:-0,2** 10 
0 


EM 
2.60+0.09* м. = 


Thus, it may be concluded that in animals the exposure of the acupuncture point corres- 
ponding to Zusanli one in humans to MMW radiation in a number of cases resulted not only in 
primary tumor inhibition, but also in anti-metastatic effect in case of regional and distant 
metastatic spread. Corrective influence on energy metabolism disturbed by the tumor process, and 
.hence the increase of organism's antitumor activity may be regarded as one of the possible 
mechanisms of the observed antitumor effect of MMW radiation. Obtained data showed that the 
MMW frequency characteristics were specific for different animal species. 


B.2. IMMUNOMODULATING EFFECT OF MMW 


The MMW influence on cellular and humoral immunity of CBA line mice exposed to 
gamma-irradiation was studied. 


MATERIALS AND METHODS 


30 mice of CBA line were exposed to y-quanta in a dose of 4.0 Gy. On the third day after 
irradiation (during recovery of immunoreactivity) the one half of the animals was subjected to 
MMW action by generator "Рогоз-1” [15] on the nasal point for 5 minutes during 5 days. The 
second half of the animals was kept in the analogous conditions but without MMW action. Then 
all the irradiated animals and the control group (10 intact animals, devoid of y-irradiation action) 
were immunized by intravenous injection of the suspension of bacterial cells of staphylococcus 
aureum Wood 46 inactivated by heating. The titres of antibodies in the experimeftal animals’ blood 
serum were determined on the peak of immune response (on the 6th day after immunization) with 
the method of solid phase immunofermentative analysis (ELISA). 

The same animals’ blood was taken from the tail vein before the immunization and 6 days 
after it, and the blood smears were prepared and fixed for 1 minute in 6095 aqueouz acetone. Then 
the activity of succinate dehydrogenase (SDG, 1.3.99.1) and lactate dehydrogenase (LDG, 1.1.1.27) 
was defined in lymphocytes by the method [17] modificated by A.S.Janovskaya. After the 
dehydrogenase's definition the smears were fixed again in the vapors of 10% neutral formaline, and 
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the acid nonspecific phosphatase (AP, 3.1.1.6) was determined by the method of simultaneous 
azocombination [19]. Using the peculiarity of AP distribution in the cytoplasm of lymphocytes [20] 
the contents of T- and B-lymphocytes, B-activated lymphocytes as well as the subpopulation of T- 
lymphocytes: helpers (Th), suppressors (Ts), and killers (Tk) — were found. 


RESULTS AND DISCUSSION 


In Fig.2 the results of investigations of this set of experiments are shown. As it can be seen, 
the irradiation with y-quanta dose of 4.0 Gy led to significant reduction of serum antibodies titre 
comparing with the control (intact animals). The using of MMW resulted in an increase in 
antibodies titre, but after 5 sessions it did not achieve the indices value of control group. 


E E 
2.0 2.0 
1.5 1.5 2 
1 
1.0 1.0 
0.5 
d 
Еа тозаест Re Ce a ME СИЛ, 
1/300 1/1000 1/3000 
/ / Dilution van Ne а 


Fig.2. The titres of antibodies acting against staphylococcus.  . Fig.3. The titres of antibodies acting against staphylococcus 
in the CBA line mice under combined action of y-quanta and in the CBA line mice under MMW action: 1 — control; 2 — 
ММУ: 1 ~ control; 2 — y-quanta (4.0 Gy); 3 - y-quanta + MMW radiation 

MMW radiation t 


Thus, the data obtained bear witness that the MMW action on animals damaged by ionizing 
irradiation enhanced their humoral immunity. 

It was interesting to clear up if MMW acted similarly upon the intact animals or upon the 
animals with pathology only. 

For this purpose 15 mice of CBA line were irradiated by MMW during 5 days (for 5 mini) 
on the nasal point, then the animals were immunized according to the methods described earlier. 
After immunization the animals were treated with MMW again Юг 2 days. Then a 2 days’ 
interruption was made, and again MMW were applied during 2 days. 6 days after immunization the 
titres of antibodies acting against staphylococcus were defined in blood serum in 15 irradiated and 
10 intact animals. 

In Fig.3 the data obtained are shown. Thus, as in the previous experiment, a stimulation 
(average 4095) of antibodies synthesis against staphylococcus was obtained. 

Simultaneously the influence of MMW on the cellular immunity in animals damaged with 
ionizing irradiation was also under investigation. In the blood smears taken from the experimental 
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animals the correlation between SDG and LDG activities in the lymphocytes characterizing their 
functional activity, and the subpopulation of T- and B-lymphocytes: T-helpers, T-suppressors, T- 
killers, and B-activated lymphocytes were determined. The data are presented in Table 9. 


Table 9 
The changes in the blood indices in mice of CBA line under combined action of y-quanta and MMW 


Groups of Intact animals y-quanta of 4.0 y-quantat MMW 
animals Gy/mouse 
Blood indices Before After Before After Before After 
immunization | immunization | immunization | immunization | immunization | immunization 


B-activated 
lymphocytes, 95 3.20.2 9.90.6 6.6+0.3** | 18.4+1.1** | 12.6+3.2* 14.6+2.3** 


8.6+0.4 6.9+0.2 6.2+0.1** | 3.6+0.9** 


Note: * — p < 0.05 as compared to animals damaged with y-irradiation; 
** — p < 0.05 as compared to intact animals. 


As seen from the data in Table 9, the MMW action after y-irradiation has significantly 
improved the indices of mice blood compared to the group of animals subjected to the ionizing 
radiation damage only: the correlation of Th/Ts was 5 times higher due to an increase in helpers' 
amount; the content of B-activated lymphocytes responsible for antibodies synthesis was as much 
as twice higher. t 

It should be mentioned that in this group of animals after preceding immune system 
stimulation induced by the MMW action the immune response to the injection of staphylococcus 
was not so clearly expressed as that of the group subjected to y-radiation only. The Th/Ts ratio and 
the amount of B-activated lymphocytes increased only thrice as much. At the same time the killer's 
content had decreased twice as much what could be evidently connected with their enhanced 
activity and their death resulted from the destruction of target-cells. 


C. IMMUNE STATUS CORRECTION UNDER MMW ACTION IN PATIENTS WITH 
PRETUMOR PATHOLOGIES 


The influence of low-intensity MMW radiation on immune status of patients suffering from 
polyps and polyposis gastrica was examined. 


MATERIALS AND METHODS 


The therapeutic potencies of MRM method were investigated by the model of polyps and 
polyposis gastrica treatment. The patients’ group numbered 162 persons and consisted of 84 men 
and 78 women of the age of 28—76 years. The method of MRM was applied to 142 patients: in the 
lst group (62 patients) — after an endoscopic polypectomy, in the 2nd group (24 patients) — 
without any other treatment. In the 3rd group (20 patients) only a surgical treatment (an 
endoscopic polypectomy) was conducted. 

MRM treagnent was applied to 86 patients with an industrial tuning generator (53.37— 
78.33 GHz, 10° /ст') and to 56 patients with the same generator using a power attenuator 
(1071 W/cm’). The treatment was applied at the resonance frequency individually chosen according 
to the patients’ sensor reactions. MMW action, was carried out through the acupunture point loci, 
which are specific to these pathologies. тн the energy of MMW radiation to a given 
acupunture zone the flexible teflon waveguide 5 cm long and with 1.8x3.6 mm? cross-section was 
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used, whose loose end was placed at a distance of 4—6 mm above the skin surface. The duration of 
one session of MMW treatment was equal to 20 min. During one session from 2 to 4 points with 
the lowest indices of electrodermal resistance were irradiated: CV12, St21, CV6, CV10, CV9, Liv3, 
SP4, SP9, BI18, St36, ВІ20. 

The MRM treatment consisted of 8—12 sessions. The treatment was conducted daily or every 
other day, and during the process the endoscopic control of epithelium recovery processes of the 
injured surface was made. 

The MRM treatment efficiency was evaluated by clinical condition of the patients, 
gastroendoscopic control and blood cytochemical indices. For cytochemical investigations of 
peripheral blood the smears were fixed in 60% water solution of acetone for 30 s. The activity of 
succinate dehydrogenase (SDG), lactate dehydrogenase (LDG) [20] and acid phosphatase [19] in 
lymphocytes were determined. Using distribution peculiarities of acid phosphatase in Iymphocytes' 
cytoplasma they were divided into T- and B-lymphocytes and also undifferentiated cells, which did 
not have markers of T- and B-cells [20]. The content of T-, B- and undifferentiated lymphocytes 
was determined before treatment and after 5th and lOth sessions. As control, the blood smears 
taken from 10 healthy donors without the treatment were used. A total of 540 micropreparations 
was studied. Data processing was made by Student criterion. 


RESULTS AND DISCUSSION 


After analyzing the condition of patients of the first group it was established that MRM 
method was an aid to effective stopping of pain in polypectomy zone during the first two days after 
the surgical treatment. Besides, this method substantially facilitated the epithelium recovery 
processes of the injured surfaces of gastric mucosa. The. duration of recovery period depended on 
the size of injury. The complete epithelization of postpolypectomy wounds ‘during 14 days was 
noticed among 90.4% of the first group of patients. Among 9.6% of patients wound healing was 
registered close to the end of one month. Under MRM treatment the functional state of gastric 
mucosa injured by polyps was substantially improved, as it was seen from decreasing of mucus 
amount, stopping of inflammation signs. When MRM treatment was used after polypectomy the 
recovery of gastric motility which is usually lowered in the patients with polyposis gastrica was 
noticed by means of endoscopic control. Analogous changes of motility and functional state of 
gastric mucosa took place under the influence of MRM treatment given to the patients of the 
second group. In the third group (surgical treatment only) healing of the injured areas was noticed 
after 10—40 days. 

The results of cytochemical investigations of the blood cells in patients undergone MRM 
treatment are represented in Table 10. 


Table 10 
Cytochemical indices of blood cells under combined MRM treatment of polyps and polyposis gastrica 


P er E Before treatment 


17.6218 
8.8204 
144209 


Note: Dehydrogenase activity is expressed in number of diiBhormasan granules per cell [17]: 
* — reliable results (p < 0.05). 
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It was shown that among the patients suffering from polyps and polyposis gastrica the 
metabolic processes in blood lymphocytes were substantially disturbed, and in lymphocytes an 
energetic shift towards glycolysis took place. Considering that lymphocytes are mostly: classified as 
the cells with aerobic metabolism and their major demand for energy is provided by respiration, the 
abovementioned shift leads to disturbance of lymphocytes’ functional state [18]. Besides, according 
to our results not only the disturbance of lymphocytes protective function occured, but the relative 
content of T-lymphocytes was lowered as well. This led to a T-system deficit and thus enhanced 
the possibility of polyps’ malignancy. 

According to our investigations the discussed indices for healthy donors were: SDG activity 
— 10.9+0.4, ГОС activity — 4.0+0.1, T-lymphocyte content — 69.93.2066. By comparison with 
Table 10 it follows that after 10 sessions of the treatment the indices in the patients and healthy 
donors came closer. 

The problem of minimal and maximal dose settings with the use of new methods of 
electromagnetic treatment has a special interest. Regarding to MRM treatment this problem is 
concluded in finding the optimal values of radiation power and exposure time. We partially solved 
this problem by investigation of MRM efficiency when the MMW radiation power density was 
lowered to 10 ^W/cm/, what is close to the values of the human body's electromagnetic field. It 
was established that the decrease in MMW power density from 107 W/cm’ to 107“ W/cm’ led to 
the increase in therapeutic effect from 71.9% to 93.7%. It was manifested in a decrease in recur- 
rence rate of polyps after polypectomy, in the improvement of gastric motility and secretory functi- 
ons. These improved gastroendoscopic results correlated with blood cytochemical indices (Table 
11). 

Table 11 
The increase in immunocorrective action of MMW treatment of patients with polyps and polyposis 


gastrica by decrease in radiation power 
Blood cells indices Before treatment After 5 sessions 
1.8+0.15* 


MMW radiation power density 10^ W/cm” 
55.5=3.4* 
27.2=4.5* | 


17.3+4.2* 


SDG/LDG activity 
T-lymphocytes, 96 38.2=4.2 
B-lymphocytes, % 31.3 4.5 


Undifferntiated cells, % 2 Ро a0 Sa 8.9 


MMW radiation power density 10^ W/cm’ 

0.40.01 1320.12" 
T-lymphocytes, % 38.624.5" 

B-lymphocytes, % 40.1+3.7 36.7+5.9 28.3+4.0* 

Undifferntiated cells, % о 4.5x0.2* 

Note: * — reliable results (p < 0.05). 


As is seen from our experiments, glycolysis in lymphocytes was prevalent among all the 
patients, but after the 5th session of MRM treatment with the use of attenuator cell respiration 
predominated over glycolysis. After the complete course of 10“ W/cm’ MRM treatment 
SDG/LDG ratio was 6 times higher than at the starting level. In the group of patients treated by 
MMW radiation with power dénsity of 10° W/cm^ SDG/LDG ratio was only 3 times higher than 
the starting one. After the 10th session of MRM treatment with attenuator blood the cytochemical 
indices of blood célls were practically equal to ones of the healthy donors (SDG/LDG — 2.7+0.27; 
T-lymph — 67.9+3.2%; B-lymph — 30.0+2.3%; undifferentiated cells — 2.1+0.7%) 

Under MMW influence the changes of metabolic enzymatic activity took place, the 
maturation of undifferentiated lymphocytes was accelerated. In both cases this led to an increase of 
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T-lymphocytes content, but for attenuated power these changes took place much earlier (after 5th 
session) and were more pronounced. Thus, the immune status was basically activated and the 
imbalance of some enzyme systems was corrected. These data are correlate with the results of other 
investigations [8, 9]. 

Regardless of primary reception we could assume the following mechanisms of MMW effect: 

— correction of the human organism's specific electromagnetic fields results in self-consistent 
renormalization of its physiological functions; 

— direct action on the blood cells in biologically active zones lead to the subsequent chain 
reaction of neuro-immuno-endocrine system. 


CONCLUSIONS 


Our experimental data allow us to make the follówing conclusions: 

— MMW radiation exerts a positive corrective influence upon functional states of biological 
systems of different hierarchical levels; 

— MMW effect has a pronounced resonant character; 

— decrease in mm-radiation intensity to the values approaching quantum threshold increases 
the therapeutic effect of MMW applied to organisms’ biologically active points at selected 
frequencies; 

— normalizing action of MMW upon bioobjects damaged by ionizing irradiation could be 
used as universal preventive treatment in the environment with the constantly increasing levels of 
ionizing irradiation. 

In terms of physics of the alive the normalizing action of mm-field upon biological organisms 
could be interpreted as follows: nonadditive interaction between the human eigen electromagnetic 
field and the external field results in creation of new configuration of electromagnetic field which 
corrects in its turn the direction, velocity and the other parameters of biochemical status of the 
living organisms of different hiererchical levels due to'self-consistent mechanism. 


: sole qae ДІЯ МІЛІМЕТРОВИХ ХВИЛЬ HA БІОСИСТЕМИ РІЗНОГО РІВНЯ 
ІЄРАРХІ 


Л.С.БУНДЮК, О.П.КУЗЬМЕНКО, Н.М.РЯБЧЕНКО, Г.С.ЛИТВИНОВ, С.О.БОБРОВНИК, 
Л.І.МАЦЕЙКО, І.В.ГОРБЕНКО, О.Р.ФЕДОРОВИЧ 

У роботі показано, що міліхвильове випромінювання нетеплової інтенсивності діє 
нормалізуюче на ріст клітин лінії L, пошкоджених іонізуючим випромінюванням. 
Встановлено, що міліхвильове випромінювання при дії на біологічно активні зони справляє 
коригуючу та імуномоделюючу дію на функціональний стан людини та тварин з різними 
патологіями. Терапевтичний ефект підвищується. при зменшенні потужності міліхвильового 
випромінювання аж до величин 10° Br/cm’. 

Запропоновано механізми відновлення функціонального стану біосистем різних рівнів 
ієрархії під дією міліхвильового випромінювання. 


КОРРЕГИРУЮЩЕЕ ДЕЙСТВИЕ МИЛЛИМЕТРОВЫХ ВОЛН НА БИОСИСТЕМЫ 
РАЗЛИЧНЫХ УРОВНЕЙ ИЕРАРХИИ т 


Л.С.БУНДЮК, А.П.КУЗЬМЕНКО, Н.Н.РЯБЧЕНКО, Г.С.ЛИТВИНОВ, С.А.БОБРОВНИК, 


Л.И.МАЦЕЙКО, И.В.ГОРБЕНКО, О.Р.ФЕДОРОВИЧ 


В работе показано, что милливолновое излучение нетепловой интенсивности оказывает 
‚ нормализующее действие на рост клеток линии L, поврежденных ионизирующим 
излучением. 
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Установлено, что милливолновое излучение при воздействии на биологически активные 
зоны оказывает коррегирующее и иммуномодулирующее воздействие на функциональное 
состояние человека и животных с различными патологиями. Уменьшение мощности 
милливолнового излучения вплоть до величин порядка 10 ^ Вт/см” усиливает его терапевти- 
ческий эффект. 

Предложены механизмы восстановления функционального состояния биосистем различных 
уровней иерархии под воздействием милливолнового излучения. 


2. 


3. 
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Effects of Ion- Conformational Interaction in Biomembrane Channels 
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The existence of discrete levels of conductivity is one of the fundamental properties of ion channels. 
It is generally accepted that the presence of stable conformational states of channel-forming protein 
is determined by its structural features. Every conformational state has its own conductivity level. 
From this standpoint, the existence of the mentioned conformations is not determined by ion flux, 
though current-voltage dependencies and kinetics may depend on it. 


In this paper we show that the interactions between channel charged groups and passing through it 
ions (ion-conformational interaction — ICI) determine kinetics of the channel closing and opening 
(gating kinetics). Thus it depends not only on voltage but bulk concentrations as well as the ion type 
and charge of gates. This approach is approbated on experiments where the open state probability- 
voltage dependencies shift when the potassium concentration is changed. 


Sometimes near reverse potential channel functions in the regime of single conformational state 
when the closed state probability is practically equal to zero. These experiments may be explained by 
the model of channel self-organization in nonequilibrium conditions when the current exceeds a 
certain critical value and the ion-conformational interaction is strong enough. 


(Received August 8, 1994) 


DESCRIPTION OF ION TRANSPORT AND CHANNEL FUNCTIONING 


Ion channel is a protein macromolecule piercing membrane and forming an aqueous pore in it 
(see [1]). It is known that a channel-forming molecule includes up to 4000 amino acid residues and 
it is generelly accepted that the pore is formed by four homologous subunits. Each of them consists of 
five (sometimes seven) hydrophobic segments and one positively charged segment (S4) which serves as 
a voltage sensor. This sensor “feels” the external electric fields and determines the channel kinetics 
(probability of open (closed) state, lifetime, etc.) as a function of the membrane potential. 

On the other hand it was shown that the properties of ion selectivity, permeability, conductivity 
are determined by the interactions between limited number of charged amino acid residues and 
passing ions [1, 7—8, 11—12]. There is an important example, the shift of the open state kinetics- 
voltage dependency under different factors such as: 

— alterations of the permeant ion type [2, 3]; 

— glterations of the ion concentration in potassium [4] and calcium channels [1]; 

— alterations of the gating charge value [5]. 

If the external electric field (membrane potential) acts directly on voltage sensor, interpretation 
of these experiments can hardly be found. We have to take into account the passing ion charge actions 
on the sensor. Indeed, the intensity of the external electric fields in the channel run up to 10° V/cm 
whereas the ion fields in its vicinity run up to 10° V/cm. 

Thus, the ion channel functioning may be treated in the following way. There are some charged 
groups in a channel which move under the action of thermal fluctuation and may close the pass for 
ion flux or remain it free. On the other hand the ion movement may be treated in a similar way 
because it move under the similar fluctuations. 

Let us describe the Hamiltonian for ion transfer in a channel 


H(r,x) = Hont) + Ha (x) + Vi 1 Xx) (1) 
where r is the ion coordinate and x is the channel conformational vector determining conductivity 
levels. The first term determines the energy of ion in the channel. The second term determines the 


energy of the channel without passing ions, and the third term determines the energy of ion-channel 
interactions depending on both ion localization and channel conformation [22]. 
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As we mentioned above ion transport may be presented as random processes described by a 
System of stochastic equations [22]: 


d? S ах, _ дН(х,г) г) 


М ар +k dt m OX, on T EE ex 1), (2.a) 
2 
n eh e UL gy | (2.5) 


— + шо 
а "а & 
where 7, and 7, are viscosity coefficients, ¢, and £ are random forces, m is the ion mass and M, are 
the conformational groups masses. 


Let us proceed on the assump- ^P. 
tion that the conformational motion is д s X \ 
slower in comparison with the charac- 1 


teristic times of single ion motion. In 
this case, the displacements made by Eg 
separate ions sum up and the position 7 
of the molecular (conformational) out 
groups will be determined by average 
values of ion concentration in the 
channel. On the other hand the change 
of molecular group position in the 
channel may produce a substantial ef- 
fect on the magnitude of ion flux. 
Often the ion motion in a chan- 
nel is described with the help of ap- 
proach based on Eyring’s absolute re- 
action rate theory. According to this > 
approach the ion channel may be rep- 
resented as an energy profile with al- r 
ternating minima and maxima, Fig.1. 
The process of ion transport in such a | Б.І. Energy profile E(r) of the channel with one binding site and two barriers 
system is the sequential random jumps 
of ions from one minimum to another under action. of the thermal fluctuation of the environment. It 
is assumed that ions are localized in so called binding sites - minima of potential energy E. The 
balance equation may be written for the occupancy № 


dN 
dt 


= М) + М - N (Wio + Wi) (3) 


where 
E-E 
Й, = Оу ex EE ^em 
и 4) 
is the probabilities of ion jump from well / to m over barrier lm (т = / + 1); T — the absolute 
temperature; К — the Boltzmann constant; ø — the external electric field potential; e — the electron 
charge; z — the ion valency. Stationary ion flux passing through the channel may be found from 
expression: 
МИ о = NoWoM12 i» 
Wio +Wi2 0) 
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In assumption of one-dimensional slow conformational mobility it can be shown that the 
system of stochastic differential equations (2) may be replaced by the balance equation (3) and the 
stochastic differential equation for the movement of the conformational group in the effective confor- 
mational potential U(x) (see Eq.(1)): | 


d dU | 
[eee uero. 
A (6) 
009 zHj(X)* | Євою. Ñ (x)dx. 


where c is the time of conformational damping, and N (x) is the normalized quasistationary value of 
the occupancy М. Thus 

Хо М +, =1. (7) 

On the last stage of description it is necessary to define the connection between conformational 

changes and energy profile. Evidently, the latter is to be changed as a result of moving of the charged 

groups of the channel-forming protein. We will be restricted only by the energy barriers mobility 

because the value of the well depth in a channel (except the outermost positions) does not affect the 


value of the ion flux (5). Basing on the aforesaid the mobility of the energy barriers have to be defined 


in the following way: 


Ej; = Ер + Vis (X) (8) 


Here Ej) is a constant. 


The kinetics of ion channel func- 
tioning is determined by the effective con- 
formational potential U(x) (see Fig.2). 
Indeed the stationary disctribution func- 
tion for conformational variable has the 
following form: 

U(x) 

FAD CA Bs 9) 
where C is the normalization coefficient 
and D is the noise intensity. The form of 
this potential, its voltage dependencies may 
be found from experimental data the way 
Hodgkin and Huxley done it. 

Тһе system of balance equation (3), 

Xi X: X, X stochastic differential equation (6) and the 

equation for conformational mobility (8) and 

Fig.2. The ion channel functioning in effective conformational poten- the expression for potential function Н NE) 

tial U(x) terms. completely determines the ion transport 
.. through the channel. 

According to the generally accepted point of view the rate constants and the open state prob- 
ability are determined by expressions: : 


Plc 10 
De -n|--.n, (10) 
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[ U(xx) - (хз) [U(x,)) -U(x3)] 
аер 09200). p ~ exp) VE) | (11) 


THE INFLUENCE OF ION CONCENTRATION ON CHANNEL KINETICS 


In developing approach the influence exerted by passing through channel ions on its function- 
ing may be taken into account. This influence was investigated by some authors [7—8, 15—19] but 
usually the authors followed discrete model of the channel where the conductivity levels corresponded 
to a set of discrete channel conformations. Our approach differs from the previous ones by continuous 
conformational variables and gives us the opportunity to study channel transformations including the 
number of conformational states as well as to describe some experimental effects. 

Let us study the effect of the shifts of dependencies between the open state probability and the 
membrane voltage at the alterations of potassium concentrations [4]. The stationary open state prob- 
ability can be determined from Eqs.(10)-(11) we obtain using the effective conformational potential 
definition (6): 

* 
U(xy) -U(xy) = А + [Pi язд (12) 


x 


where A is the difference H (х )-Н. (х,). It can be a function of the voltage and can determine the 
channel kinetics. The voltage Sensor "hypothesis i is based on this point of view [1]. But we will see that 
this value may be a constant as to function of the voltage because of ion-conformational interaction 
presence. It allows us to describe not only the voltage dependencies but the concentration dependen- 
cies too. 

Using the perturbation theory we can determine the deviation of the stationary state caused by 
ion-conformational interaction. Indeed expanding the channel potential H (х) near the unperturbed 
positions we obtain the expression for new stationary points: 


1 dVQ(x) Ñ 
L á N Cio) (13) 


Ч XENI 


Xi = Xoi — 


where x, are the unperturbed points, k, are parameters of channel potential (constant within our 
approximation). 

Thus we obtained the channel conformational states (Eq.(13)), its effective conformational 
potential and now we can determine the channel kinetics dependencies against voltage and concen- 
trations. The stationary value of open state probability can be found from Eq.(10) 


1 
nos 
oo 
“Le exp) GU ES) (14) 


On the other hand, the barrier mobility (8) alters the ion flux value (5). This approach allows to 
make the following expression for ion flux: 
-y/2 N 2 e" -N ое DE 


£z Joe eh int x) + ne A 


(15) 


, a ; | Ej - Ey | 
where J, is the normalization coefficient, 7 = SPP y = zeg/kT. 
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Fig.3. Comparison of numerical simulations with the experimental data observed in [4]. 
Parameters: С, = 100 mM, 1 — С, = 100 mM, 2 — С, = 325 mM, V „7 a/x, а = - 80, 
Rep — 20,2; eddy A= = 615,7 40 


The results of numerical simulations of the experiments mentioned above and approximated 
within approach proposed are presented on Fig.3. It could be considered as evidence that the ion- 
conformational interaction can serve as a tool of the channel detection of membrane voltage and 
determine its kinetics. 


ION-CONFORMATIONAL INTERACTION CAN INDUCE APPEARANCE OF ADDITIONAL 
CONDUCTIVITY LEVELS 


There are some additional examples of effects similar to the mentioned above: 

— for the isolated membrane fragments from molluscan neurons it was established in [14] that 
independently on potential, temperature, and concentration the rate constant (1/7,) approaches zero 
when the current is approaching the zero value (see Fig.4a); 

— it was shown in [13] that the lifetimes of the open and closed states ATP-sensitive potassium 
channels of myocardial cells take respectively the maximum and minimum values near the reverse 
potential and decrease or, correspondingly, increase during the deviations from it to more negative as 
well as to more positive potential (see Fig.5a). This effect takes place when outer potassium concen- 
tration is varied. 

The shifts of the kinetics-voltage dependencies could be explained by the approach described in 
the previos parts of the paper. But there is a question which remains unanswered: why the probability 
of closed state approaches to zero value in the equilibrium state (reverse potential) in spite of the 
alterations of any parameters (temperature, ion concentrations, etc.)? 

Such stability of the effect may be explained in the framework of notions about channel self- 
organization under non-equilibrium conditions. It was shown earlier [6, 10, 21] that the strong ion- 
conformational interaction can induce арау of the new conformational state. Let us study a 
simple situation when interaction potential V. (x) are linear functions of the conformational variable. 
And the channel without ICI has a single conformational state. The dynamically stable states are 
detemined by the following nonlinear equation (see (6)): 


x 2x, N(x) (16) 
(x, is an ICI parameter) which may have from one to three solutions. It means that the system 


functions either as the monostable one (one solution) or the bistable one (three solutions). The 
corresponding current-voltage dependencies are plotted on Figs.4b and 5b. One can see that the 
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Fig.4. The disappearance of the closed state observed in the experiments [14] (a) and its theoretical interpre- 
tation (b). The current in the experiments was changed by temperature at different fixed membrane potential: 
1—9-7-20mV,2—9-0mV,3—9-20mV,4—9-40 mV. p, < o < o, — region of bistability 


-80 0 80 У Мету 


[K*],, 7300 ам 
o-[K *], =10 mM a-[K*], =150 mM 
»-[K ^], =50 mM a-[K *], 2300 mM 
Fig.5. The disappearance of the closed state observed in experiments [13] (a) and its theoretical interpre- 


tation (b). 9, < 9 < 0, 9, < p < p, — regions of bistability 
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channels with one and two bistability regions are possible. The corresponding channels are described 
by one. and two binding sites models correspondingly [10]. The stability of the above effects may be 
explained by the limitations for the existence of bistable regime (three solutions of Eq.(16)). One of 
them is the limitation for the driving force and ICI intensity. In a single binding site model it means: 


Ix ( Nae" г №") > [геном (17) 


Thus the self-organization may take place only if the ion flux exceeds some threshold value and 
the channel functions under strong nonequilibrium conditions. Hence, for any varying of control 
parameters C,, С. , V for the approaching of the current to the zero value, one of the channel state 


disappears and the system turns to the monostable regime of functioning. 


CONCLUSION 


The effects of ion influence on channel functioning was investigated by some authors [6—11, 
15—22]. In this paper we showed that the ion-conformational interaction can produce the voltage 
depencies of the channel kinetics describing gating processes in some experiments. In this case the ion 
flux is the tool of channel kinetics voltage control, not the direct action of membrane potential on 
voltage sensor. Such approach gives us opportunity to take into account the ion concentration influ- ' 
ence as well as its type and the charge of the gate also on channel kinetics, voltage-current dependen- 
cies. The appearence of an additional conductivity levels may be induced by interactions mentioned 
above if the ICI intensity is strong enough. All the above effects have experimental corroborations. 
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ЕФЕКТИ ІОН-КОНФОРМАЦІЙНОЇ ВЗАЄМОДІЇ Y КАНАЛАХ BIOMEMBPAH 


Г.Є. ВАЙНРЕБ 


Одним із фундаментальних властивостей іонних каналів є наявність дискретних рівнів 
провідності. Загальновизнано, що наявність стійких конформаційних станів білка- 
каналоутворювача зумовлена його структурними особливостями. З цього погляду існування 
вказаних конформацій не визначається іонним потоком, хоча вольт-амперні характеристики 
та кінетика канала можуть від нього залежати. 

У цій роботі показано, що взаємодії між зарядженими групами канала та іонами, які проходять 
крізь нього (іон-конформаційна взаємодія), визначають кінетику відкривання та закривання 
(кінетику воротних процесів). Таким чином, вона залежить не тільки від напруги, але й від 
концентрації іонів у розчині, їх типу та заряду воротних груп. Цей підхід було апробовано на 
експериментах, де залежність ймовірності відкритого стану від напруги зсувається при зміні 
концентрації калію. 

Тноді поблизу потенціала реверсії канал функціонує у режимі єдиного конформаційного стану, 
коли ймовірність закритого стану дорівнює практично нулю. Ці експерименти можуть бути 
пояснені самоорганізацією іонних каналів у нерівноважних умовах, коли струм перевищує 
деяке критичне значення, а іон-конформаційна взаємодія є достатньо сильною. 
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ЭФФЕКТЫ ИОН-КОНФОРМАЦИОННОГО ВЗАИМОДЕЙСТВИЯ B КАНАЛАХ 
БИОМЕМБРАН 


Г.Е ВАЙНРЕБ 


Одним из фундаментальных свойств ионных каналов является наличие дискретных уровней : 
проводимости. Общепризнано, что наличие устойчивых конформационных состояний 
каналообразующего белка обусловлено его структурными особенностями. С этой точки зрения 
существование указанных конформаций не определяется ионным потоком, хотя вольт- 
амперные характеристики и кинетика канала могут зависеть от него. 

В этой работе показано, что взаимодействия между заряженными группами канала и 
проходящими через него ионами (ион-конформационное взаимодействие) определяют 
кинетику открывания и закрывания (кинетику воротных процессов). Таким образом, она ; 
зависит не только OT напряжения, HO и OT концентрации HOHOB в растворе, их типа и заряда 
воротных групп. Этот подход был апробирован на экспериментах, где зависимость вероятности 
открытого состояния от напряжения сдвигается при изменении концентрации калия. 
Иногда вблизи потенциала реверсии канал функционирует в режиме единственного 
конформационного состояния, когда вероятность закрытого состояния равно практически 
нулю. Эти эксперименты могут быть объяснены самоорганизацией ионных каналов B 
неравновесных условиях, когда TOK превышает некоторое критическое значение, а ион- 
конформационное взаимодействие достаточно сильно. 
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Nonlinear Effect of Dynamical Self-Organization in Macromolecular Systems 
Which Is Induced by Photocontrolled Electron Flux 


A.O.GUSHCHA, A.A.DOBROVOLSKY, M.T.KAPUSTINA, A.V.PRIVALKO, 
V.N.KHARKYANEN. 


Scientific Research Center "Vidhuk", Kiev, Ukraine 


The bifurcation of optical transmission coefficient dealing with monostable-bistable functioning 
regime transition was found for the first time in experiments with the reaction centers (RC) of 
purple bacteria Rb.Sphaeroides. The phenomenon was observed both in suspensions of isolated RC 
and Lengmuir films on their basis under the conditions of acting optical pumping. The theoretical 
interpretation of the effect in terms of dynamical self-organization of structure with photoinduced 
electron flux was developed. The phenomenon can be considered as nonequilibrium phase 
transition of the first order. 


( Received August 2, 1994) 


INTRODUCTION 


Of the known self-organization processes in the realm of nonlinear physics are those of 
radiation-material interaction in the laser which represent, after H.Haken, an example of a new 
class of phenomena — nonequilibrium phase transitions [1]. The essence of the laser phenomenon 
is that electromagnetic radiation, interacting with atoms, induces a redistribution of populations of 
the ground and excited states which, in its turn, affects the coherent irradiation-absorption 
processes, ensuring the system's feedback. When the external pumping reaches a certain critical 
value, the system goes over into a qualitatively new state, representing. the dynamical self- 
organization effect. 

In Ref.[2] it was theoretically shown that the interaction of the charged particles’ flux with 
the macromolecular structure which forms a biomembrane ionic channel could initiate similar 
effects of self-consistent partial rearrangement of the macromolecular structure, bringing about the 
change of the transmitted ionic flux intensity. 

In the present work we have undertaken an experimental investigation into the effects of 
dynamical self-organization in the structure of the purple bacterium photoreaction center (PRC) 
macromolecule through which the photoinduced electron flux passes. 


I. EVIDENCE OF THE ELECTRONIC-CONFORMATIONAL INTERACTION IN THE 
(PIGMENT — PRIMARY ACCEPTOR) PRC SECTION 


The PRC macromolecule is a molecular structure which consists of the protein matrix with 
embedded molecules of bacteriochlorophyll, pheophytin, ubiquinone and the Fe atom [3-4]. The 
charge photoseparation scheme is pictured in Fig.l. Under the actinic light, the pigment Р 
(chlorophyll dimer) becomes exited. Subsequently, in 2—4 ps the charge separation takes place, 
resulting in the electron transfer of nearly one hundred percent efficacy to the pheophytin 
molecule Ph (not depicted in Fig.1). At the next fast stage (150—200 ps) the electron comes to 
the primary quinone acceptor Q,, with almost the same efficacy. The subsequent slower stage of 
the electron transfer to the secondary quinone acceptor Ов and i ; possible recombination with the 
hole at the pigment is of particular interest since their efficacy depends on the positions of certain 
labile structural elements. The very rearrangement of these elements, in its turn, is initiated by the 
photoinduced redistribution of electron density [5]. Thus, under stationary illumination, the time of 


О Ов > О „Оз transfer increases by several orders during a decrease in the light exposition time; 
while the time of P*Q, — PQ, recombination, diminishes 3—5-fold under the same conditions 
I6]. | 
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x Treating the PRC specimens with 
РОО, orthophenanthroline makes it possible to remove 
the acceptor Оз. Such being the case, the 
P*O} з РО, recombination stage is the only 
transfer section controlled by the labile structure, 
with the populations of states P'Ph Q, and 
P PhO, being, under condition of stationary 
illumination, neglected as compared with those of 


PPhQ, and P* PhQ%, states (henceforth we label 
the latter by np and пл, respectively). The 
corresponding transition scheme is 


p 
РО ОВ 


DM S ро 
PQ, QUIS 04 (1) 
Fig.1. Simplified scheme of the charge photoseparation and 1 


transfer іп PRC of Rb.Sphaeroides. Кгѕ are the correspon- B : 
ding rates of the electron transfer; np, n4, пв ате the electron Under real conditions, the chlorophyll dimer 


population densities of the pigment and the primary and photoexcitation constant J (proportional to the 
secondary acceptors, respectively actinic light intensity) limits the electron transfer 
in the chain PBPhQ, > P'BPhQ, > РВРИ Q, — PBPhQj and the rate constant k; є J. The 
recombination constant k., [6] is a function of the structural variables which specify the positions 
(c.nformational coordinates) of mobile PRC elements. The electron recombination from the 
primary quinone acceptor to P is of a tunnel nature, so in Ref.6 the conformational coordinate is 
identified with the distance of this tunneling. For & ;(r) the dependence 


r 


k ,(r) = ne є (2) 
is employed, where r is the change (related to photoinduced reconstruction of the system) of the 
distance between P and Q, 7 is the rate of recombination in non-light-adapted PRC, 


h 
п = Вх is the characteristic parameter of the electron tunneling, m is the electron mass, and 
m, 
A is the tunneling barrier height (in Ref. [6] т, was assumed -1А ). 
Evolution of the donor and acceptor populations, np and ил, respectively, is governed by the 


set of balance equations 


dn v 
Ё. = [пр + пе "n, 
dt 
; (3) 
dn, n 
mene Un In 
| dt 7] A P 
which hold the normalization condition 
Ф np(t) - ny(t) 21 (4) 
The motion equation for the conformational variable r can be written as 
dr dr 
М = +9 = -kr + foo), 5 


where M is the generalized mass, q is the friction coefficient, and k is the elastic c8nstant of the 
conformational coordinate. 


36 


NONLINEAR EFFECT OF DYNAMICAL SELF-ORGANIZATION IN MACROMOLECULAR SYSTEMS 


In the state P'Q' the conformational coordinate is affected by an additional force Л, 
emerging due to the charge redistribution, so that 


1, in the state P*Q^ 
@(t) = 
0, in the state PQ 


Set (3) taken together with Eq.(5) determines the self-consistent evolution of the electronic 
and conformational variables. However, Eq.(5) is of stochastic character, though thermal 
fluctuations was not directly introduced in it. 

The most proper way of describing such a system is to proceed to Fokker-Planck’s formalism 
for the distribution function which determines the joint probability of the electron and 
conformational coordinate localization. In the present work we are interested to determine the 
system behavior under stationary light action only. Therefore, we restrict ourselves by ensemble- or 
time-averaged values of the conformational coordinate which are given by (5) regarding that 
©(t) = ny (the line stands for averaging): 


FF (7) 


(6) 


where п, = fo is the maximally possible displacement. As follows from Ref.6, the quantity Г 


increases under photoexcitation, that is why we set / > 0 and, consequently, г, > 0, what is of 
defining importance for our theory. 


Let us assume that the quantity 7e " ~ ge . This should be valid in the case when the 
electronic transitions proceed much faster than the conformational rearrangements. Henceforth it is 
handy to omit the averaging line of the conformational variable. 

Under stationary illumination, from Eqs.(3) we have 

I 
ат | (8) 
Із пе ^ 

The joint solution of Eqs.(7) and (8) as a function of the actinic light intensity J is showed in 
Fig.2. One can readily see that in case of r, < 4r the dependence is of monotonic character. If 
г, > 47, then the plot takes the characteristic S-form, manifesting the bistability area appearance, 
see Fig.2a. 


9.2.4 4 8 
lg I, au. —— r, ац — 


Fig.2. а) The dependence of stationary donor-acceptor distance г on optical excitation intensity J for 
different parameter values: у = 1, љ = 0.1 (curve 1); п, = 4 (curve 2); re = 8 (curve 3); 

b) Nonequilibrium conformational potential о of the system for the following parameter values: 7 = 1, 
то = 8, I = 0 (curve 1); Z= Д = 0.015 (curve 2); J = 0.25 (curve 3); J = Ду = 0.062 (curve 4) 
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The curve of stationary values of the conformational coordinate has two bifurcation points: 
= (Л, n) and My = (Д, rg), with the lower branch (0, Му) corresponding to stable states of the 


r 


system with п, ««1 and ge ” >> Г, and the upper one { М, œ} — to those with пу «1 and 
r 
пе З ««I. 

Relative thermodynamical stability of stationary states of the S-curve upper and lower 
branches can be determined by introducing the nonecuilibrium conformational potential of the 
system, depended on the actinic light intensity: ; 

SUP feu. zl 


е E гаю о "| 
p Il m o o 


Fig.2b shows a collection of curves v(r5 Г) with the actinic light intensity as a control parameter 
whose values correspond to zero intensity, bifurcation points and the middle of the bistability area. 

It is obvious that in this area the stationary distribution function of conformational 
coordinate 


(9) 


ud dz A 


ptr: -4 хе (10) 


has two maxima; with the actinic light SR growing, density redistributions take place among 
them. 
Applying the abovementioned theory to the case specified in Scheme (2), one can evaluate 


Г а 
the quantity 5, proceeding from the fact that illumination of the PRC without the secondary 
0 


acceptor increases the recombination rate by the factor of 3 + 5 [6], i.e. 


Гр 
k (0 m 
£u coquus (11) 
Kk о) 
Y 
Therefore, —^ ~ Ш + 16 < 4, so the bistable behavior does not emerge in the given section of the 


To 
PRC electron transport chain. Another situation can be observed in the section between the 
primary and secondary quinone acceptors, Q,Q0, > 0,05, where the electron transfer rate 
changes by several orders of magnitude. 


II. EXPERIMENTAL OBSERVATION OF BISTABLE BEHAVIOR OF THE TRANSMISSION 
COEFFICIENT OF PURPLE BACTERIA PRC SPECIMENS 


Employing optical methods, we have undertaken an experimental investigation in order to 
clear up the possibility and conditions of the emergence of bistability effects in electron transfer 
within the Rb. Sphaeroides PRC Q4-Qs section. 

The optical set-up scheme is pictured in Fig.3. As a probing radiation, that of санок 
lamp KGM-150 was employed, after passing through attenuating neutral light filter NS-3, 
monochromator MDR- й and filter KS-17 which isolated the radiation of wavelengths of operating 
diffraction orders. Lense Ll served for shaping of the probing radiation front. Varying the supply 
voltage of the incandescent lamp, it was possible to vary the probing radiation (4 — 865 nm) 
intensity within a wide range; however, in all experiments its value did not exceed 
10° quantum/ (cm? sec). " 
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Fig.3. Scheme of the set-up 1 — cuvette with the specimens under study; 2 — block of controlling 
the thermostat temperature; 3 — photomultiplier FEU-157; 4 — preamplifiers; 5, 9 — lamps 
KGM-150; 6 — laser LG-79-1; 7 — polarization variator of light intensity; 8 — rotating mirror 


As a radiation collector, a photoelectronic multiplier of FEU-157 type was used, having 
almost the constant spectral sensitivity characteristic in the whole visual range up to 4 — 880 nm 
and being relatively to its noise characteristics on a par with such devices as FEU-136. In front of 
FEU-157 filter KS-19 was placed so as to prevent the possible coming of scattered exciting 
radiation to the photocathode. Recording signals was made in the dc conditions. * 

Supplementary stationary excitation of the specimens under study was carried out into their 
absorption bands 550—650 nm with He-Ne laser (4 = 632,8 nm) or incandescent lamp KGM-150. 
In the latter case filter SZS-22 extracted the operating range of the exciting light spectrum 
(A, = 450—650 nm). Pumping radiation intensity was varied from 0 to 10° quantum/(cm'sec) by 
neutral light filters or controlled polarization attenuator and was monitored by reference detection 
channel in all experiments. Excitation beam was pointed to a specimen at angle = 30° with respect 
to probing radiation. Beam aperture was maintained to be sufficient for complete overlapping of 
both beams on the surface and in the volume of specimens under study. 

Specimens (suspensions of isolated RC's of 10? cm? concentration) were placed into cuvette 
of 2 mm in thickness. Measurements were made at room temperature. In particular cases, to 
maintain constant temperature when making prolonged experiments, cuvette of 1 mm in thickness 
was placed into a specially designed thermostat on Peltier elements. Precision of temperature 
maintenance in the cuvette volume was + 0.5°C. Controlling thermostat temperature was carried 
out by an independent power supply. All measuring, controlling and monitoring signals were 
connected with CAMAC system interfaced with PC/AT-286. 

As is known, photooxidation of the bacteriochlorophyll RC pigment leads to modification of 
the absorption spectra [9]. The most striking changes take place in the spectrum region 700— 
900 nm (see Fig.4). Under saturating illumination, i.e. when the density of exciting light quantuma’ 
flux per second is close to that of RC's in the specimen volume, the absorption band 4 = 865 nm 
bleaches almost completely. The characteristic isobestic point for the specimens under 
consideration is in the region 4; = 808 nm. Under any exciting light intensfties the absorption 
coefficient at 4; does not change. We measured the absorption spectra of RC's from Rb.Sphaeroides 
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at various levels of optical excitations. In one 
series of such experiments the excitation 
intensity gradually grew from zero .to its 
maximal value ІЗ = 5-10? quantum/(cm'sec) 
(Fig.5a), whereas in another series — it 
diminished from Г“ to zero (Fig.5b). Every 
Spectrum was recorded during a 20 sec time, 
and the specimen's time delay at a given 
excitation intensity until recording the spectrum 
was 4 min and more. This guaranteed the 
completion of conformational transitions related 
to redistribution of the electron density in RC 


A, a.U. 


М-и 


2 


700 800 900 molecules and ensured adequacy of the 
recorded spectra to specimens’ stationary 
2 , hm excitation pre-history. Let us take notice of the 


identity of curves 1 and 5 in Figs 5a and 5b and 
Fig.4. Characteristic region of the Rb. Sphaeroides PRC of the pronounced shift of all curves shown in 
bsorpti t; in th ited state (1 d und . : : : гі 

маковий Nomine d) que (D and under Fig.5b and picked at intermediate excitation 

; intensities relatively to the corresponding curves 

in Fig.5a. This indicates the presence of the memory in the RC electronic-conformational system, 

with this memory manifesting itself on reaching certain threshold levels of specimen's optical 
excitation. 


a) 


A, 3.l. — > 
"Pu 
n 
un + wn = 
А, аль > 
їз 
) 
A 
> Ww Ne 


850 860 870 850 860 870 
a, nm ——» a, nm ——> 


Fig.5. The dependence of the optical absorption spectrum region from 850 to 880 nm on the pre- 
history of the optical excitation (A, = 450 + 650 nm) during its rising from zero to a fixed value J (a) 
and lowering from the maximum level to J (Ы). The curves in both figures correspond to the following 
values of J (in quantum/(cm’sec)): 1 — 10'°; 2 101; 3 — 5.101; 4 — 10; 5 — 3-10 


Below we present experimental data on studying the dependence of the RC optical 
transmission coefficient at ^ = 865 nm on the excitation pre-history, with the excitation increasing 
continuously from zero to a given maximal value /"* and decreasing subsequently to the initial 
zero value. The rate of changing excitation was selected thus as to not violate the quasistationary 
conditions. Measurements of the optical transmission А at every J were made after passage of the 
time sufficient to establish the system equilibrium (1—5 min). In every experiment the reached 
maximal excitation /" was different. The results for suspensions in a cuvette of 2 mm in thickness 
are presented in Fig.6a. For /"* < 2- 10!! quantum/(cm'sec) the downgoing and upgoing branches 
of the R(/)-dependence form a hysteresis loop. The intersection of both curves in the vicinity of 
I = 4. 10" quantum/(cm’sec) allows us to assume this value to be a bifurcation point of the (labile 
molecular structure + photomobilized electron flux) system. One of the reasons which prevented 
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Fig.6. The dependencies of optical transmission coefficient R at 2 = 865 nm on pre-history of 
the excitation (4, = 633 nm) for suspensions (a) and PRC Lengmuir films (b). 

The arrows indicate the direction of the excitation intensity varying І is: 1 — 103; 2 — 
3: 103; 3—8. 108 (quantum/(cm'sec)) 3 


recording the intersection of both branches of the A(J)-dependence бог 1» 2 · 10! quantum/(cm'sec) 
was absorption of the exciting radiation in the specimen's volume, according to Buger's law. We 
attempted to diminish this factor by conducting the research with Lengmuir's films in RC. The 
films of 10 and 30 monolayers were obtained by Sheffer's method in a bath with bidistillated water 
under the surface pressure of 25 din/cm’. The experimental data recorded under the same 
conditions that those for RC suspensions are shown in Fig.6b. The first bifurcation point 
corresponds to J = 4: 10 uantum/(cm'sec), i.e. to the same value as in the experiments on RC 
suspensions. For J > 5-10" quantum/(cm’sec), one can see the narrowing of the hysteresis loop 
and the two branches tendency to merge. Apparently, the excitation intensity range 
(2° 10! + 5-10") quantum/(cm'sec) may correspond to the second bifurcation point in the 
electron transport chain between О, and Оз in Rb.Spaeroides. 


ПІ. EVIDENCE OF ELECTRONIC-CONFORMATIONAL INTERACTION AT THE PRIMARY 
QUINONE — SECONDARY QUINONE STAGE 


To interpret the obtained data, let us consider the following three-level scheme of electronic 


transitions: 
РО Ов, РОО, > Р'0,05 
n "E (12) 
T k; | 


Note that, according to Ref.[12], k..<< Кз, and we henceforth set k., = 0. The constant of 
the nonradiative electron transfer between the quinone acceptors can be written as 


pt 2 у 
(х) = бе“ (13) 
Е 48 - A AB | ) 
where x = —= and ЕЧ is the conformational reorganization energy, A45 is the energy 
AJ EU) k,T 


difference, К; is Boltsmann's constant, T is absolute temperature, and c is a certain constant. 

Obtained by Marcus [10] for electronic transitions, the expression (13) corresponds to the 
electron tunnelling during thermally activated adjusting of the conformational variable, in 
accordance, with Franck-Condon's principle. 
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Under a conformational rearrangement induced by the electron coming to the secondary 
quinone acceptor, the quantity x changes from ж to the value x, which corresponds to the 
permanent electron localization at Qg (x) > x,). Let us assume that the conformational. variable 
stationary displacement related to the actinic light intensity / сап be read as 

x(I) = хо +(X — Хо) hg(1) (14) 
where пв(1) is the stationary Оз population. With np +n, + ng = 1, the balance equation set 
corresponding to scheme (12) results in the following stationary populations of quinone acceptors: 


I 
ою (у * ке Й + ce) 
Ic on (15) 
xs (у + се) Й + cer I 


where 


"ONUS о Wen 
a Кз mm ne й UN ks (16) 


Joint solving (14), (12) and (7) allows us to find out the stationary values of conformational 
variables ғ, x and populations п, and пв as functions of the control parameter 7. However, in the 
room-temperature case the problem gets somewhat simplified, since it follows from the experiment 
that n, << 1, and we can neglect the change of the distance between P and Q}, i.e. r ғ 0. As can 
be readily shown, theoretical estimations prove this conclusion. Therefore, we further set у = const. 
The existence of the bistability area follows from examining the transcedent equation 


Ice” 


о Wr aa o (17) 
n сеї у се? и 
The dependence x(/) is depicted in Fig.7. The case (a) conforms to the situation in which 
Сен (18) 
жи 
where £1 = а , 


lg] —— : lg] —— 


Fig.7. The dependence of stationary values of conformational coordinate x on the logarithm of actinic light 
intensity / Depending on values of Structure parameter ©, the three types of dependencies may be harnessed: 
(а) 2 = 105 (b) 2 = 101; (с) С = 10” 
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(Хо +X) + xe eae Xo +X, 
ее. 
3 
Such being the case, the curve х(Г) is simply connected and monotonic. The conformational 
coordinate change is little. 

In case (b), when ¢ > С i. there emerges the second part of the curve x(/), separated from 
the first one by a forbidden values' interval. The bistability area is characterized by only one 
bifurcation point. As the parameter © grows to a certain value 2.1, the unstable branch of stationary 
values of x from the lower part merges with the stable one of the upper part, giving the infinitely 
distant second bifurcation point. Lastly, in case (c) there takes place a singly-connected curve with 
two bifurcation points, allowing us in principle to observe a hysteresis behavior of the object under 
photoexcitation intensity changing. Experiment [11] shows that in dark-adapted specimens, after a 
short pulse of actinic light population of the primary and secondary quinone acceptors takes place 
almost equally. This yields ¢ e = 10. However, according to the data of Ref.[7], the electron 
localization at О, initiates fast (~ 1—10 ms) relaxation processes. They reduce the energy gap 
between electronic levels at О, and Q5 consequently increasing the range of xy. Let us assume 
Се 210! + 10-4. 


In light-adapted specimens, according to [6], ¢ e* ж 10. Using these data and also 
assuming that ¢ x 10Ps? (molecular collision frequency in condensed matter), we have: хо = 5,5; 
Xo я 4,29; С 10°. Of course, we can slightly vary these parameters. 

The experimentally measured quantity is 


(19) 


n té ех du 


rod 


n+ ce sce) 


R-R, ‚4.01. —— 


0 1 
lg I(a.u.) ——> lg I (a.u.) —> 
Fig.8. Theoretical value of optical transmission coefficient A, Fig.9. Theoretical value of optical transmission coefficient 
corresponding to the stable branches (solid line) of the R, taking into account the normal distribution of the PRC 
bistable curve of stationary PRC states for ¢ = 10°, п = 10, parameter xo about the mean 5.5 with the variance 0.07. 
y= 10. The arrow indicates the hysteresis behavior of the Different hysteresis curves relate to different maximal values 
curve A(7)- Ry under quasistationary increasing / іо а certain of actinic intensity 


maximal value and subsequent vanishing 
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where A, is the background value of the transmission coefficient, а is the parameter which gives the 
contribution of photomobilized electrons to the specimen transmission coefficient. Fig.8 shows 
quantity (20) calculated for the parameters corresponding to the curve in Fig.7c (contribution of 
stable branches of the bistable curve x(/) only is taken into account). 

Comparison of the theoretical (Fig.7c) and experimental (Fig.6a, b) curves can be made only 
after averaging the theoretical result over different values of parameters Хо, Хо» ©, / Such an 
averaging leads to a dispersal in positions of the bifurcation points, i.e. to smoothing the hysteresis 
curve. Our computation has shown that the normal distribution of the parameter x, with the 
variance 0.07 about the mean 5.5 leads to the hysteresis curves presented in Fig.9, what 
qualitatively agrees with the experimental data. Apparently, for complete quantitative agreement 
one should take the conformational coordinate anharmonicity into consideration. Besides, a notice- 
able role can be played by additional states with the electron at the secondary acceptor which 
manifest themselves by the appearance (under prolonged illumination) of long-time recombination 
components. ) \ 
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НЕЛІНІЙНИЙ ЕФЕКТ ДИНАМІЧНОЇ СОМООРГАНІЗАЦІЇ В МАКРОМОЛЕКУЛЯРНИХ 
СИСТЕМАХ, ВИКЛИКАНИЙ ФОТОКЕРОВАНИМ ПОТОКОМ EJIEKTPOHIB 


О.О.ГУЩА, А.ОДОБРОВОЛЬСЬКИЙ, М.Т.КАПУСТІНА, А.В.ПРИВАЛКО, 
В.М.ХАРКЯНЕН 


У експериментах з суспензіями ізольованих реакційних центрів (РЦ) пурпурових бактерій 
Rb.Sphaeroides i виготовленими на їх основі ленгмюровськими плівками в умовах оптичного 
збудження вперше виявлено біфуркацію коефіцієнта оптичного пропускання при переході 
від моностабільного до бістабільного режиму функціонування. Запропоновано теоретичну 
інтерпретацію одержаного ефектуд Як явища динамічної самоорганізації структури та 
фотоіндукованого потоку електронів, що може розглядатись як нерівноважний фазовий 
перехід І роду. 


НЕЛИНЕЙНЫЙ ЭФФЕКТ ДИНАМИЧЕСКОЙ СОМООРГАНИЗАЦИИ В 
МАКРОМОЛЕКУЛЯРНЫХ СИСТЕМАХ, ВЫЗВАННЫЙ ФОТОУПРАВЛЯЕМЫМ 
ПОТОКОМ ЭЛЕКТРОНОВ 


А.О.ГУЩА, А.А.ДОБРОВОЛЬСКИЙ, М.Т.КАПУСТИНА, А.В.ПРИВАЛКО, 
В.Н.ХАРКЯНЕН 


В направленных экспериментах с суспензиями изолированных реакционных центров (РЦ) 
пурпурных бактерий Rb.Sphaeroides и приготовленных на основе РЦ ленгмюровских пленок 
в условиях действующего оптического возбуждения впервые обнаружена бифуркация 
коэффициента оптИческого пропускания при переходе от моностабильного режима к 
бистабильному. Дана теоритеческая интерпретация полученного эффекта как явления 
динамической самоорганизации отоиндуцированной электронной плотности, что может 
рассматриваться как неравновесный фазовый переход І рода. 
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Calcium Signal: Registration in Cells and Quantitative Estimation 


V.L.ZYMA, O.M.DYACHOK 
Department of Biophysics, Kiev University, Kiev, Ukraine 


Microspectrofluorimeter is described. This setup allows to record excitation and emission spectra 
of fluorescent dyes, such as fura-2 and indo-1, in different types of cells as well as to measure 
Са?" — transients in the muscle. Recording of the spectral characteristic of fluorescent dyes in 
single-beam and dual-beam configurations is shown. 


Results suggest that a quantitative estimation of intracellular calcium in intact cells by means of 
fluorescence spectra of calcium-bound (LM) and calcium-free (L) indo-1 forms only could not be 
sufficient. The quantitative estimation also requires to consider the protonated form (LH) and 
fluorescence related to binding indo-1 with some proteins in cytoplasm. The calcium estimation 
in human platelets were realized by computer analysis of complex fluorescence spectrum indo-1 
on "reference" spectra of established indo-1 forms. 


Fluorescent properties of indo-1 under binding with Са?" in protein solutions has been 
investigated. No changes in indo-1 fluorescent parameters have been recorded in solutions of such 
proteins as trypsin, calmoduline, papain, lysocyme. Significant changes of indo-1 fluorescent 
parameters have been observed in solutions of histones, bovine and human albumins. Obtained 
results allow us suggest that indo-1 binds well with protein sites containing positively charged 
amino-acid residues and hydrophobic *pockets". Destruction of the compact structure of serum 
albumin after addition of 8 M urea shifted indo-1 fluorescence spectrum and it approached to 
spectrum indo-1 in water. Fluorescence spectra of indo-1 in dioxan suggested that binding of 
indo-1 takes place in hydrophobic microenvironment of protein globule. 


(Received September 17, 1994) 


INTRODUCTION 


Calcium ions are the most important second messengers which transduce information from 
excitable cell membrane to different cell structures. Thus, a lot of cell functions become activated: 
contraction, motility, secretion, the processes of synthesis and others [1, 2, 3]. Application of Са?" 
sensitive fluorescent dyes (fura-2, indo-1, fluo-3) [4, 21, pH-sensitive dyes [6] carried out radical 
improvements in methods of intracellular calcium | Са? ^], and pH [H ]; measuring and monitoring 
of changes during cell's activitity. In commonly used sets for registration of Са" -transients in 
nervous and muscle cells double-wave fluorescent parameter В = IA,/IA, has been measured. It is 
the ratio of fluorescence intensities at the wavelengths selected by filters [4, 7, 8, 9]. But spectral 
origin of Ca^'-signal is more complex due to binding of fluorescent dyes, apart from Ca^, to 
cytoplasmic proteins, and also due to the changes of dyes spectral properties after changes in 
cellular pH [6, 9, 10, 11, 12]. Analysis of fluorescence spectra testifies the existence of four 
individual forms of indo-1: a free form of indo-1 (L) pH > 7.5, a calcium-bound form (LM), a 
protonated form (LH) pH < 7.0 and a form of indo-1 bound with some proteins [10, 13]. The 
mechanism of indo-l-proteins interaction is not quite understandable. There is a question as to 
how this interaction can influence calcium determination in cells. 

Thus, more universal sets to measure fluorescence intensity at a determined wavelength as 
well as spectral characteristics of calcium indicators are needed. 
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The aim of this study was to investigate influence of medium polarity and microenvironment 
of proteins on indo-! fluorescent parameters and to clear up the nature of protein sites binding 
indo-1. 


MATERIALS AND METHODS 


I. FLUORESCENCE SPECTROSCOPY 


Fluorescence spectra were digitally recorded in a 1-cm-width quartz cuvette with а SDL-2 
(LOMO, Russia) luminescent spectrometer interfaced. with a microcomputer. For the scanning of 
emission wavelengths the microcomputer drove stepping motors with a 1 nm step. Fluorescence 
spectra were recorded between 360 and 600 nm (excitation 350 nm). Both excitation and emission 
slits were 5 nm. The fluorescence was detected by a photomultiplier tube (PMT) operated in 
photon count mode. Photon counts collected from PMT through photon counter were integrated 
over 0.1 s at each wavelength and transmitted to the computer. 

. In order to analyze indo-1 fluorescence we used the technique based on computerized 
resolution of numerized fluorescence spectra [14]. 


Il. MICROSPECTROFLUOROMETRY 


A. OPTICAL COMPONENTS 


The setup has been designed using  microspectrofluorimeter MFTX-2M (LOMO, Russia) 
(Fig.1). The optical device consists of a microscope (objective O, dichroic mirrors ОМ, ОМ, and 
DM;) апа monochromators M, and M.. Both excitation and emission spectra are recorded dué to 
rotation of difraction gratings by means of two stepper motors SM, and SM; carried by a special 
device SMC. Dyes fluorescence is excited by xenon lamp 150 W (I). Fluorescence output is 
monitored by means of photomultiplier PMT; and PMT.. 

Equipment can work in single- or dual- beam designs. 


Fig.1. Schematic diagram of the experimental set-up: I — light output; О — an objective; F,, Fə and 
Ез — filters; DM,, DM; and DM; -- dichroic mirrors; M; and M». — monochromators equipped 
with difraction gratings; SM; and SM» — stepper motors; PMT, and РМТ. photomultiplier tubes; 
MD — mirror disk; D — a photodiode; PC, and PC, — photon counters; SMC —‘a Stepper motor 
controller; АРС — analogous-to-digital convertor; К — a flow chamber with a biological object; En 
— channel for entering fluid; Ex — channel for exiting fluid; T — a force transducer, t — a 
thermoresistence; P — an external loading; AT 286 — an IBM compartible microcomputer 
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1. Single-beam design: 


1.1. Measurement of fluorescence integral intensity. In this design excitation filter Е, passed 
light from lamp I and cut off filter F; transmitted fluorescence but not transmitted excitation light 
are used. Fluorescence is monitored by means of photomultiplier PMT). Dichroic mirror DM, is 
entered in the light path and DM, ОМ, are not used. 

1.2. Recording of excitation spectra of fluorescence. In this design xenon lamp I is mounted 
instead of РМТ,. Dichroic mirrors DM, and DM, are not used. Light from lamp pass through 
monochromator М, with Е; and objective О onto a biological object. Fluorescence from the object 
was collected by the same objective passing through DM; and filter Е; to PMT. Automatic 
scanning of excitation spectrum is carried out by stepper motor SM). 

1.3. Recording of fluorescence spectra. In this design DM, and DM; are used. Excitation 
light from I passes through filter F;, objectives O excites the fluorescence of biological object. The 
emission light then passes through a monochromator M, and is focussed onto the photocathode of 
a РМТ). In this design, mirror disks MD are not entered in the light path way. 


2. Dual-beam design: 


2.1. In the first mode dichroic mirrors DM,, DM;, ОМ», mirror disks are entered in light 
path. Fluorescence is excited by light I through filter F,. Fluorescence output I4, and 12», selected 
by monochromators M; and М; is alternately referred to PMT, due to rotation of MD. 
Dualwavelength fluorescent ratio R = L4;/I4; is measured. The data is analyzed from fluorescence 
spectra. This method is preferable to indo-1. 

2.2. In the second mode lamp I is used instead of PMT, and MD, DM, and DM; are 
entered, but DM, is removed. For selection 4; and 2, in the excitation spectrum monochromators 
М, and М; are employed. The excitation light passes through the objective O and illuminates the 
biological object. Fluorescence emitted from object passes through F; and is monitored by PMT). 
In this method two excitation wavelength fluorescent ratio R = I4;/IL4; is measured. This 
registration is preferable to fura-2. 

The bandwidth of the incident light beams L4, and I4; аге 10 nm. The mirror disk revolves at 
150 Hz, allowing to measure fluorescent ratio R at intervals of about 6.6 ms. 

The microscope is equipped with UV transmissive optics. Using UV objectives magnification 
58x with numerical aperture (NA) of 0.8 (water immersion) and 10x, 0.2 NA allows to excite the 
fluorescence of biological objects without considerable losses. To increase setup sensitivity and 
decrease the noise, photomultiplier tubes are operated in photon count mode. 


B. TENSIOMETRIC DEVICE 


A flow chamber K with muscle strip is placed on the stage of the microscope. Isometric 
tension is measured using force transducer 6MXIC (Russia). Initial force of the muscle strip is 
provided by loading P. Temperature in the chamber is measured by thermistor probe t. 


C. DATA COLLECTION 


The PMTs output is monitored by photon counting circuitry PC, and PC; which is 
interfaced with an IBM compatible AT 286 microcomputer. The microcomputer reads and then 
clears the counts accumulated in the photon counter at the time when the opaque sectors of the 
mirror disks are in the light path, and so the count reflects emission or excitation at one 
wavelenght only. Such analogous signals as tensiometric force, temperature and etc., monitored in 
the experiment are digitized to a 12 bit resolution and transmitted to the microcomputer via a 8- 
channel analogous-to-digital convertor (ADC). Synchronisation of data acquisition from the 
photon counters, ADC and stepper motors control is carried out by the entry of signals from 
photodiode with the frequency of 150 Hz. Programs for the system written with use of the *C" 
programming language allows to measure and display optical and tensiometric information for any 


48 


CALCIUM SIGNAL: REGISTRATION IN CELLS AND QUANTITATIVE ESTIMATION 


setup design. Data collected in the experiment saves in ASCII format and can be analyzed and 
calculated by commercial software. : 


III. MATERIALS 


Indo-1 (pentapotassium salt), indo-1/AM (pentaacetoxymethyl ester of indo-1) and EGTA 
(free acid, purity 97%) were obtained from Sigma, USA. Pluronic F-127 (BASF) was obtained 
from Molecular Probes, USA. HEPES and bovine serum albumin were from Fluca Chemical, 
Switzerland. Human serum albumin was from Reanal, Hungary. Trypsin was from SPOFA, 
Czechya. Papain, lysocyme and dioxan were from Reachim, Russia. Histones and calmodulin were 
kindly supplied by Prof.S.Khrapunov, Dept. Genetics, Kiev University and by Dr.V.Danilova, 
Inst.Physiology, Kiev University respectively. All solutions were prepared with double-distilled 
water (resistence 2 MQ:cm). All measurements in. protein solutions were carried out in the 
presence of 5mM HEPES, pH - 8.0 at the room temperature. 

Human platelets were obtained from Kiev Inst. of Cardiology. Cells were suspended, at 
about 7x10? cells/ml in solution containing (in mM): 150 NaCl, 2.7 KCl, 2 CaCl, 1 MgCl, 0.37 
NaH;PO, 5 glucose, 0.35% (w/v) bovine serum albumin, 50 unit/ml heparine, 5 HEPES; 
pH = 6.55. Cells were incubated for 1 h in the similar medium, without heparine, pH = 7.4 
contaning 5 иМ indo-1/AM and washed thrice by the same solution. Cells were maintained in 
suspension by frequent gentle shaking of the cuvette. 

Fluorescence spectrum of ап indo-1 treated human platelets was obtained in a SDL-2 
spectrofluorometer with a 1 cm cuvette at room temperature using the excitation wavelength of 
350 nm. 

Smooth muscle thin strips were dissected from guinea pigs of either sex, weighing 0.200 to 
0.250 kg. The muscle strips (with the length of 3 mm, width of 1 mm, thickness of 200 шт) were 
loaded with indo-1/AM (5 uM) by 1 h incubation at 36°C. The noncytotoxic detergent, pluronic 
F-127 (0.0596 w/v) was premixed in the loading physiological salt solution (PSS) to help dissolve 
indo-1/AM into PSS. Then strips were rinsed 
several times with PSS, mounted in flow chamber 
K between two stainless-steel hooks. One hook 
was attached to a force transducer and from the 
other end it was connected to external load P to 
create the initial tension of muscle strip. Previosly 
calibrated force transducer measured isometric 
tension development in Newtons. Indo-1 fluo- 
rescence excited through Е, with Ax = 350 nm. 

The composition (in mM) of PSS was as fol- 
lows: 135 NaCl, 5.9 KCl, 2.5 CaCl, 1.2 MgCb, зим 
11.5 glucose, 5 HEPES; pH = 7.4. 


RESULTS AND DISCUSSION 


Ca’*-signal has been monitored in dual- 
beam design 2.1. Ca-transient signal and isomet- 
ric tension response of thin strip of caecum to K^ 
depolarisation are shown in Fig.2. Ca-transient E =] ыы 
signal has a complex shape. In the beginning, a Vom M ECI 50 c 
rapid increase in calcium signal was observed 
which corresponded to phasic tension of caecum. Fig.2. Са — transient signal and isometric tension response of 
The second nome in Ce] and ong oni jiu om ages, She 
tension of caecum were followed by a slow Asc 350 nm JOS cius : 
decrease in Са" -signal. 


.3.5 
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Fig.3. Spectral analysis of the fluorescence of indo-1 — 
treated human platelets (without cell autofluorescence) 
(1) using the spectra of the calcium-free (2) and the 
calcium-bound (3) indo- 1. Aex = 350 nm; [+ — arbitrary 
units. Inset: WR — weighted residuals calculated as the 
ratio of the residuals versus the absolute value of the 
random noise superimposed on the signal [14] 
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Fig.4. Fluorescence spectra of different forms of indo- 1: 
1 — calcium-free, (L-form); 2 — calcium-bound, 
(LM-form); 3 — BSA-indó-1 complex (LP form of 
indo-1) without calcium ions; 4 — ВЅА-іпао-1 
complex with calcium ions; 5 — BSA-indo-1 complex 
without calcium ions ingghe presence of 8 M urea 
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The whole fluorescence spectrum of indo-1 
treated human platelets was recorded. An attempt to 
resolve this fluorescence spectrum using the spectra 
of the calcium-bound and calcium-free indo-1. was 
unsuccessful, as illustrated in Fig.3. The results of the 
resolution should be rejected since both the graphic 
and numeric estimators of the fit between the 
resolution and the experimental spectrum indicated 
severe distortions in the determination. This overva- 
luation could occur due to the fluorescence, about 
450 nm, of one or several fluorochromes related to 
indo-1, whose spectra have not been introduced in 
equation set. These observations are in accordance 
with the findings of Bancel et al. [10] and Owen et 
al. [11]. The intracellular calcium concentration in 
platelets was calculated according to the equation by 
Grynkiewicz et al. [4]: 


[Ca^*], = Ка. B -(R- Rmin) / (Rmax – R) 


where К = LjyLs- 1.304 — fluorescence ratio recor- 
ded from the platelet suspension; В= Iss /T,455— 2.218 
— the ratio of the fluorescence of indo-1 in the Ca’*- 
free form to the Ca?'-saturated form recorded at 495 
nm; Kd = 282 nM [13] — calcium dissociation 
constant of the indo-1. 50 JM digitonine was added 
to the cell suspension to record Rmax (2.109) followed 
by the addition of 10 mM EGTA to obtain Rmin 
(0.115). The intracellular calcium concentration in 
the platelets was 920 nM. 

The results obtained suggest that the created 
setup allows us to monitor Са" -signal from cells and 
tissues in time. This equipment can be used for 
measuring the spectral characteristics of fluorescent 
dyes, for investigation of interaction between proteins 
and probes and for quantitative estimating of [pH], 
and [Ca^]; changes. 

The fluorescence spectra of calcium free form 
of indo-1 (L) in the solution without Ca^' in the 
presence 2 mM EGTA and calcium bound form 
(LM) in the solution with the saturated calcium 
concentration (2 mM) are shown in Fig.4. Emission 
spectrum maximums at 2, = 477 nm for L and 
Am = 415 nm for LM forms were detected. Data ob- 
tained in this study are similar to those presented in [10, 
13], despite of some differences for L form, which 
according to these works has position 4,, = 485 nm. 

Fluorescence spectra of indo-1 were measured 
in solutions in the presence of various proteins. In 
the calcium free solution in the presence of 40 uM 
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bovin serum albumin the emission peak maximum 150000 г зни т pcm 
shifted towards short wavelengths (emission peak 
Am = 445nm). It is testified that BSA-indo-1 complex 
(LP- form of indo-1) is constituted (Fig.4, spectrum 
3). After adding of 2 mM Са" in this solution, 
fluorescence spectrum is considerably shifted towards 100000 | 
Short wavelengths at Amn = 412 nm (Fig.4, spectrum 4). 
Destruction of the compact structure of serum 
albumin after additon 8 M шеа shifted шдо-1 
fluorescence spectrum and. it approached the | 
spectrum without BSA (Fig 4, spectrum 5). Some зоо | 
shift of the indo-1 fluorescence peak (Am = 454 nm) 
was observed in the solution with low ionic strength 
и = 0.004 in the presence of 50 uM histones (Fig.5, 
spectrum 1). Increase in ionic strength to д = 0.5 led 
to shift of indo-1 fluorescence spectrum towards long А 
wavelengths and approached the L-form indo-1 350 400 450 500 550 600 
spectrum. In solutions of such proteins as calmodulin, Anm 

lysocyme, trypsin, papain the changes of indo-l Fig.5. Fluorescence spectra of of indo-1: 
fluorescent parameters were not observed. This shows 1 —. in the solution with histones (50 uM), ш = 0.004; 
the absence of binding indo-1 by these proteins. Our 2 — in the solution with histones (50: ИМ), д = 0.5; 
data rather differ from the results of study [6], where 3 — in the water without any proteins and calcium ions 
indo-l spectral parameters changes were established зососо, rs 2n 


in trypsin solution. Fl 
Obtained data suggest that polarity of microen- f A 
viroment of indo-1 bound with BSA and histones can 120909 F | 
be determined through cyclic hydrophobic structures | \ 
of indo-1 as well as through electrostatic interaction... | | \ | 
between four carboxyl groups with positively charged — із \ 
lysyne and arginine amino asids residues. In case of | 
indo-1 binding by histones the electric charges in- 150000 
teraction are predominant, because under д = 0.004 ; 
histones have disordered conformation [15, 16] and PAPA 
under low ionic strength the positively charged amino 170300 | N 
acids residues localized predominantly at the N- | / \ 
terminus of histones are accessible to interaction. This... ju A 
indo-1 form (called LP" in the following) has typical ^ |  / / Й | 
fluorescence spectrum at 2, = 454 nm. In BSA Lh bs 
solution the two forms of indo-1 are possible — LP*- à LIU о ; р. 
form (charge interaction) and LP-form due to strong 
binding in hydrophobic sites. It has been shown [10] Anm 
that two types of binding centers has been established  ,. dni 
at Scatchard diagram after binding with serum albu- adn уні осо I hae 
min. We defined little changes in indo-1 fluorescence 2 — in the dioxan with calcium ions 
spectrym followed the alteration of ionic strength 
from 4 = 0.004 to џ = 0.5. It suggests that charges interaction occuring for LP'-form is dimi- 
- nished. But predominant part of indo-1 is bound by hydrophobic sites. LP-form with an emissiom 
peak at An = 445—450 nm corresponds to that binding. 
Our model investigations suggest that the polarity of the environment exerts influence on 
fluorescent parameters and the indo-1 spectrum of LP-form is determined by hydrophobic nature 
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of binding sites. Indo-1 fluorescence spectrum in dioxan (= = 2.2) was greatly shifted towards short 
wavelengths at 4,, = 443nm (Fig.6, spectrum 1). This spectrum is very similar to BSA bound indo-1 
spectrum. Chelating of Са" by indo-l in dioxan led to strongly increasing quantum yield with 
small spectral shift towards short wavelengths (Fig.6, spectrum 2). 

LO ETE ж The strucural similarity between indo-1 
| но | and EGTA molecules should lead to the similar 
іно ^u acid-base properties, which could affect the 
calcium — binding efficiency of the fluorescent 
dye at acid pH. However, the effects of 
decreased pH on indo-1 fluorescent properties 
have not been studied in detail. Acidosis leads 
to a shift in indo-1 fluorescence spectra toward 
short wavelenghts. The spectrum consists of L- 
and LH-forms (Fig.7). We have established that 
LH-form fluorescence spectrum is the sum of 
fluorescent spectra of LH1- and LH2-forms with 
different values of protonation. Effect of acidosis 
on the fluorescence spectra indo-l was 
investigated by fluorescent ratio H at 2= 412 
nm and 4-446 nm, corresponded to the 
emission maxima of LH1- and LH2-forms. Н 
increased with acidosis in the range from pH 7.8 
to pH 5.9 (Fig.7, inset). 

Taking into account the existence of such 
forms of indo-1 as a calcium-free (L) and 


100000 
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Fig.7. Effect of pH on fluorescent properties of indo-1. 
Fluorescence spectra of different forms of indo- 1: 

1 — integral fluorescent spectrum at pH 6,5 2 — L-form; 
3 — LH-form; 4 — LH1-form; 5 -- LH2-form. Inset: effect 
of pH on the H values 


calcium-bound (LM), a protonated form (LH) 
and a form interacting with proteins (LP and 
LP") the intracellular calcium concentration in 
platelets was estimated with use of the technique 


based on computerized resolution of complex 
numerized fluorescence spectra [13]. From contributions of L-, LM-, LP- and LH-spectra of 
indo-1 and ratio of their quantum yields the calcium concentration in the platelets was 238 nM, 
calculated according to the equation by Grynkiewicz et al. [4]. 
Our investigations as well as the data given in some other papers [3—7] suggest that when 
Са? is measured in intact cells, besides L- and LM- forms, there exist intermediate LH, LP- and 
LP" -forms which can significantly complicate indo-1 fluorescence spectra in cells. Thus, numerical 
resolution of complex indo-1 fluorescence spectra is needed. 
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КАЛЬЦІЄВИЙ СИГНАЛ: РЕЄСТРАЦІЯ B КЛІТИНАХ І КІЛЬКІСНА ОЦІНКА 


В.Л.ЗИМА, О.М.ДЯЧОК 

Описаний мікроспектрофлуориметр, який дозволяє вимірювати спектри збудження i 
випромінювання флуоресцентних барвників фура-2 та 1ндо-! в різних клітинах і 
реєструвати кальцієві транзієнтні сигнали в м'язових тканинах. Показано, що для кількісної. 
оцінки внутрішньоклітинного кальцію в інтактних клітинах не достатньо використання 
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спектрів флуоресценції кальцій-зв'язаної "LM" і кальцій-вільної "1" форм індо-1. 
Необхідно ще враховувати протоновану “LH” форму та флуоресценцію індо-1, яка зв'язана 
з певними білками цитоплазми. Оцінка кальцію в тромбоцитах людини проведена за 
допомогою комп'ютерного розкладання складного спектру флуоресценції індо-1 на реперні 
спектри виявлених форм індо-1. Досліджені флуоресцентні характеристики індо-1 при 
зв'язуванні Са?" в білкових розчинах. Встановлено, що в розчинах таких білків, як трипсін, 
папаїн, лізоцим i кальмодулін, не відбувалася зміна спектру флуоресценції індо-1. Значні 
зміни флуоресцентних параметрів індо-1 спостерігалися в розчинах гістонів сивороткових 
альбумінів людини і бика. Отримані результати свідчать, що 1ндо-1 добре зв'язується з тими 
ділянками білку, де знаходяться гідрофобні "кишені" і позитивно заряджені амінокислотні 
залишки. Порушення компактної структури сивороткового альбуміну в 8 М сечовині 
викликало зміщення спектрів флуоресценції індо-1 і наближення їх до спектрів 
флуоресценції індо-1 в воді. Про зв'язування індо-1 в гідрофобному мікрооточенні білкової 
глобули свідчать спектри флуоресценції індо-1 в діоксані. 


КАЛЬЦИЕВЫЙ СИГНАЛ: РЕГИСТРАЦИЯ В КЛЕТКАХ И КОЛИЧЕСТВЕННАЯ ОЦЕНКА 


В.Л.ЗИМА, О.М.ДЯЧОК 

Описан микроспектрофлуориметр, который позволяет измерять спектры возбуждения и 
излучения флуоресцентных красителей фура-2 и индо-1 в разных клетках и регистрировать 
кальциевые транзиентные сигналы в мышечных тканях. Показано, что для количественной 
оценки внутриклеточного кальция в интактных клетках не достаточно использования 
спектров флуоресценции кальций-связанной “LM” и кальций-свободной *L" форм индо-1. 
Необходимо также учитывать протонированную “LH” форму и флуоресценцию индо-1, 
связанную с определенными белками цитоплазмы. Оценка кальция в тромбоцитах человека 
проведена с помощью компьютерного разложения сложного спектра флуоресценции индо-1 
на реперные спектры установленных форм индо-1. Исследованы флуоресцентные 
характеристики индо-1 при связывании ионов Са" в белковых растворах. Установлено, что 
в растворах таких белков, как трипсин, папайн, лизоцим и кальмодулин, не наблюдалось 
изменение флуоресцентных параметров индо-1. Значительные изменения флуоресцентных 
параметров индо-1 наблюдались в растворах гистонов, сывороточных альбуминов человека 
и быка. Полученные результаты свидельствуют, что индо-1 хорошо связывается с теми 
участками белка, где имеются гидрофобные “карманы” и положительно заряженные 
аминокислотные остатки. Разрушение компактной структуры сывороточного альбумина в 
ЗМ мочевине вызывало смещение спектров флуоресценции индо-1 и приближение их к 
спектрам флуоресценции индо-1 в воде. О связывании индо-! B гидрофобном 
микроокружении белковой глобулы свидетельствовали спектры флуоресценции индо-1 B 
диоксане. 
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Biological Mechanoluminescence (Review Article) 


V.E.OREL 
The Ukrainian Research. Institute of Oncology and Radiology, Kiev, Ukraine 


This paper gives the review of the problem of biological mechanoluminescence (ML) emerging as 
the result of mechano-electrochemical activation of biological objects. It describes the principles 
of the experimental techniques initiating ML of biologic objects with the help of electrical, optical 
and hydrodynamic tribogenerators. ML measurements of liquid blood, larvae, skin of distal 
phalanges, vessels and lymphoid cells are presented and analyzed. Attention is paid to the 
dependence of ML intensity on electrolytes, amino acids and lipids concentration. ML Spectra of 
DNA in the A and B forms are shown. The explanatory models are offered. 


(Received May 16, 1994) 


Abbreviations: 

CL — chemiluminescence 

E — electrostatic field 

ECL  — electrochemiluminescence 
I — electric current 

MD — — malone dyaldehyde 

ME — mechanoemission 

ML  -- mechanoluminescence 
PM = — photomultiplier 

RH — relative humidity 


Key words: mechanoluminescence, mechanical activation, mechanoemission, biological object. 


INTRODUCTION 


Mechanical kinds of energy such as gravitation, pressure, oscillation and friction are 
indispensable conditions for existence of every hierarchical unit of the host and have undoubtedly 
become important factors in the evolution of our planet [1]. Mechanoemissional (ME) processes 
occurring during the mechanical work are the parts of metabolism and are accompanied by 
` mechanoluminescence (ML) or sometimes triboluminescence, observed in the form of mechanical 
deformation and contact electrization of the epidermal surface, of osseous and soft tissues at the 
time of chewing food, during friction in joints of vertebrates, sexual intercourse, blood circulation, 
etc. [2, 3]. The earliest data on ML of sugar can be found in ancient Indian and Persian 
manuscripts. One of the oldest records in Europe, already referred to, comes from Francis Bacon 
in the Advancement of Learning (1605). Rubbing was the easiest way to produce electricity and the 
main characteristic of electrification was "electric light”. It is not surprising that later the students of 
electricity should turn to the study of ML. 

Let us dwell in greater detail on observations of biologic objects ML. Even the first scientific 
records reported that in darkness when a person grinds sugar one can see the light emitting from 
his mouth. A description can also be found of a generally known phenomenon which lies in the 
fact that if you stroke the cat's fur in dry weather it gets electrized and we can observe 
luminescence due to electric discharges in the air [3, 4]. In 1858 Swedish botanist Fries observed 
“glittering” of Oriental poppy in the evening. According to the botanist Molish, this was a pheno- 
menon similar to St. Elmo's fire luminescence, а form of luminous corona discharge in the shape 
of luminescent clusters appearing at the sharp edges of high objects at high intensity of atmospheric 
electrical field [5]. 


* 


"This review was presented at ihe International School of Biological Luminescence (Ksiaz Castle, Wroclaw Poland, 
1989) and I International Symposium on Mechanoemission and Mechanochemistry (Liberec, Czechosldvakia, 1991). 
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The mechanical excitation developing as a result of hydraulic stress causes a significant 
increase in the background marine bioluminescence of bacteria [6—8]. However, one should also 
take into consideration the fact that during hydraulic stress hydrodynamic resistance and cavitation’ 
accompanied by an electric discharge pure luminescence is observed [9, 10]. There are reports of 
ML during bone tissue destruction [11] and an emission in the range of 225—230 nm out of fibers 
of frog nerve and heart is recorded in the process of their mechanical deformation [12, 13]. 

Nowadays the scientists again pay attention to investigation of the biological ML. Soybean 
root and stem segments, wheat seeds produce luminescence upon mechanical injury. The results 
are consistent with peroxidative stress [14—16]. Research on this subject started mainly in East 
European countries and the results were not widely known in the West. 

This is the first Review on the topic published in Physics of the Alive. ML from biological 
objects was under discussion mainly in the papers published in 1983— 1993. 


EXPERIMENTAL TECHNIQUES 


Let's consider the design and operation of the techniques for biological ML analysis. 
Mechanical activation of biological samples is realized in three various types of producing ML 
[17—23]. The in vivo study of biological ML was conducted with the help of the photomultiplier 
(PM) FEU-130 (spectral sensitivity 200 to 650 nm), a photon counter and the high voltage supply 
unit of chemiluminometer HLMIC-01. The ML of circulating blood was studied with use of the 
hydrodynamic tribogenerator techniques in an optically transparent cylindrical cuvette with spiral 
canals through which the liquid under study was pumped (Fig.1a). The cuvette was connected with 
silicone tubes to the closed circuit of a peristaltic pump with controlled pumping velocity or to the 
large circuit of canine blood circulation. In the latter case the heart of the animal served as a 
pump. The device fór ML study in liquid phase biosamples has a highly curved concave surface of 
the canal which allows to generate the high density of electrical charges exciting ML. 


Fig.l. Various methods of producing ML: a) Techniques with hydrodynamics tribogenerator: 1 — an 
optically transparent cuvette, 2 — a spiral canal; 3 — tubes; 4 — РМ; b) Techniques with optics 
tribogenerator: 1 — shaft; 2 — clamp; 3 — paper; 4 — metallic platform; 5,— rubber, 6 — cells; 7 — 
double electric layer; 8 — glass; 9 — PM. c) Techniques with electric tribogenerator: 1 — paper with fixed 
and dried biosanfples; 2 — electric engine; 3 — electretic probe; 4 — electrode; 5 — PM; 6 — shaft 
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ML of cells was measured with the device composed of an optical tribogenerator (Fig.1b), 
PM FEU-39A (spectral sensitivity, 160 to 600 nm), photon counter USSM, high-voltage source 
B5-24A, and quartz glass (8 KY-1) which had a coefficient of passage in the range of 200—600 nm 
not lower than 80% and low-intensive self-luminescence. The tribogenerator operated as follows. 
Using an electric motor and a device for regulating mechanical pressure through shaft 1, we 
generated ML in a light-tight chamber at 45% relative humidity (RH) with angular velocity of 1200 
rpm and mechanical pressure of 0.5 kPa per 6 s. This resulted in triboactivation of the sample 
relatively to the quartz glass. ML was detected on the opposite side of the quartz glass with the 
PM. Reproducibility of luminescence measurements was assessed with radioluminescent indicators 
(brightness EY-1). Reproducibility of ML measurements was controlled by measuring ML intensity 
generated by the "void" pad following neutralization of the static charge. 

Other apparatus for ML evaluation have used the device with techniques of electrical 
tribogenerator shown in Fig.lc. Measuring electretic probe 3 is set up on the earthed current 
conductive shaft 6, connected with the electric engine, and earthed electrode 4 is fixed in front of 
photoreceiver 5. Mechanical activation of the  biosample was carried out with a 
polytetrafluoroethylene probe at the mechanical tension of 125 kPa and the rotational speed of 
1000 rpm per 5 s in the atmosphere of highly purified nitrogen at pressure of 990 hPa. After the 
even pressure of the frame with investigated sample 1 to the surface of the electretic probe there 
occurs surface triboelectrization on the touch tangents that initiates electric charge, whose 
parameters depend on the properties of materials under investigation. Then, turning at the definite 
angle, the excited part of the probe touches the operating electrode, and in the place of their 
contact occurs the maximum luminescence of the electretic probe surface and its environment 
occurs. The light flow is registered by PM using FEU-140 (spectral sensitivity, 200 to 650 nm). 

The spectral measurements were made with a panel of 26 interference light filters 242—649 
nm, each with an average penetration margin of 9 nm. 

Liquid biosamples in the amount of 0.02 ml are micropipetted on the cellulose 
chromatography paper FN-2 (Filtrak) and dried at 37°C and 45% RH for 30 min. The techniques 
with optical and electrical tribogenerators allow to study not only the liquid phase of the biological 
samples layered on paper, i.e. blood, urine, saliva, tears, lymph, Sperm, sweat, etc., but also the 
frozen biological objects, such as tissue sections and homogenates. 


MECHANOLUMINESCENCE OF BIOLOGICAL SYSTEMS 


In view of the relaxation of tribocharge, observed on the surface of corneal layer of skin 
epidermis [24], it may occur that the emergence of temporary electric fields on the skin is also 
accompanied by ML. That is why we had to carry out a comparative investigation of skin 
luminescence on the distal phalanges of the human finger. 

Before measuring the luminescence, the skin was washed with a 2% solution of boric acid 
and bidistilled deionized water. After skin drying total luminescence of forefinger distal phalange 
skin area, limited to 177 mm?, was measured in a special light-protected thermochamber with 
5596 RH. 

A series of experiments was carried out to compare the chemi- and mechanoluminescence 
from the skin of distal phalanges of human fingers at various ambient temperatures. At first, skin 
chemiluminescence (CL) was registered, followed by the measurement of total CL and ML effect 
produced as a result of graded triboelectrization at 0,01 kG/mm’ when the skin of the thumb and 
the forefinger was in contact. ML value was substracted from the total signal value and the ML- 
conditioned effect was obtained (Fig.2). 

As is seen, the observed ML value depended on the ambient temperature. Thus, at +22°C 
and -15°C it exceeded the CL level 9.5 and 2.5 times, respectively, compared with a 3.3-fold 
increase at +50°С. 
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Fig.3. Chemi- and mechanoluminescence of larvae CALLIPHORA VICINA 
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As another object of study, we used the larvae of 
Calliphora Vicina which possess a capacity for 
peristaltic movements under comfortable temperatures 
[25]. 35 larvae were kept in darkness for one day and 

: placed in a thermostatic cuvette made of organic glass 
um “Tepe ^ Vo Te over the PM. The techniques of the study comprised 3 
i series of measurements: the background of the empty 

Fig.2. Change in ML of the skin of distal phalanges cuvette, the background plus the signal of larvae 
of human fingers after their gradual triboclectrization ly minescence in native conditions, the background plus 
the luminescence signal at the mechanical tension of 66 Pa limiting the peristaltic motion of the 
insects. The values of CL and ML were determined by calculation. The rythm of peristaltic motion 
as related to the temperature was assessed by the velocity of larva passage through a glass tube with 
a diameter of 5 mm covering a 10 cm distance. The analysis of the findings shows (Fig.3) that 
when the temperature is increased from -4°C to +38°C, larval CL increases 6-fold. The value of 
ML at -4°C equalled 50% of CL intensity, while at +38°C, when larval velocity increased to 5 mm 
per second, the value of ML was 6196 of CL. 

The biological liquids and tissues are known to be the source of CL in the visible, ultraviolet 
and infrared regions of the spectrum. Whole blood, unlike plasma and serum, is not luminous 
because of the absorption of luminescence quanta by hemoglobin. However, Backai E.A. et al. [26] 
reported the CL in the circulating arterial blood. In our opinion, luminescence of circulating blood 
can be hardly explained to the full on the basis of the generally accepted concepts of the CL 
generation mechanisms of biological systems. It is rather can be conditioned by the contribution of 
ML effect to the emission of photons. In this context, we had carried out the following expe- 
riments. 

In the experiments we studied the relation between the normal canine blood luminescence 
and the one under the effect of ionizing radiation with the electrization and peroxidation processes 
[27—29]. Such an experiment was chosen in view, of the ionizing radiation-induced sharp decrease 
in the negative membrane potential of erythrocytes and the disturbance of stable imbalance 
between the peroxidation processes and the activity of bioantioxidant systems. The experiments 
were performed on 16--20 kg mongrel dogs under thiopental anesthesia (30 mg/kg, i/v). The 
optically transparent flowing cuvette was connected by means of cannulas and silicone tubes 
directly or through peristaltic pump to the closed circuit with the blood taken from the femoral 
artery and returned back to the canine femoral vein. Heparin (500 U, i/v) was administered to the 
femoral vein 10—15 min prior to the experiment. In the experiments on dogs the ML registration 
procedure was similar to that of the first experiment. Blood ML was accompanied by electric 
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discharges which were simultaneously registered as the mean values of electric current (I) using the 
platinum electrodes on the apparatus F-30. 

Products of the peroxidation were determined photometrically by the concentration of 
malone dyaldehyde (MD) [30] and the integral estimation of peroxidation kinetics was performed 
by the method of spontaneous luminescence and electrochemical one (ECL) in the alternating field 
of electric current [23]. Blood plasma samples for determination of MD concentration and 
spontaneous luminescence and ECL were taken before the measurements of ML and I. The two 
were measured directly in the closed circuit: flowing cuvette-dog with/without pump, after 
initiating the flow rate, before the irradiation and 15 min after semitotal gamma-irradiation of the 
animal (abdominal cavity) (Fig.4). 
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Fig.4. ML, I, CL, ECL and MD of canine blood after semitotal gamma-irradiation Со, 5 Gy: 
а) 1 -- ML; 2 — I (intact dog); 3 — ML; 4 — I (irradiated dog); b) 5 — ML; 6 — I; 7 — CL; 
8 — ECL; 9 — MD 


It was revealed that relatively to the intact dogs an increase in pumping speed from 40 to 
250 ml/min resulted in a nonlinear increase in ML and I intensity, the correlation coefficient being 
close to 1. 15 min after gamma-irradiation the average ML intensity and I decreased by 3996 and 
8-fold, respectively. The correlation coefficient decreased to 0.75. After irradiation in the absence 
of peristalic pump within the system, ML and I intensities decreased by 44% and 10%, 
respectively. Peroxidation indices of blood plasma of gamma-irradiated animals also decreased: MD 
concentration by 6095, spontaneous CL — by 3396 and ECL — by 4796. 

The data obtained show that during the interaction between circulating blood components, 
alongside with hydrodynamic, electrostatic, thermodynamic, acoustic, electromagnetodynamic 
phenomena there also takes place ML, which correlates with processes of triboelectrization and 
free-radical oxidation. 

Alongside with ML of whole blood we studied its mechanostimulated emission from the 
thoracic aorta endothelia stimulated with a 5-fold longitudinal stretchings of blood vessels from 25 
to 40 mm per every one second. The results of comparative study of the changes in spontaneous 
CL of thoracic aorta endothelium and ML influenced by 0.25 Gy gamma-irradiation are shown in 
Fig.5. The data analysis shows that ionizing radiation results in an increase in CL by 11896 and a 
decrease in ML by 9996 in the thoracic aorta endothelium. 
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Fig.6. Correlation between mechanoluminescence and electrostatic field strength of 
donor lymphocytes 


MECHANOLUMINESCENCE OF CELLS 


This section presents the results of the investigation on the influence of immunocompetent 
lymphoid cells mechanical activation obtained using the techniques with optical tribogenerator [31, 
32]. : 

We performed a series of experiments to test the relation between the optical emission and 
the mechano-electrochemical parameters of cells. In Fig.6 the results of measurements of the 
electrostatic field (E) and ML of chromatographic paper FN-2 with and without peripheral blood 
lymphocytes are compared. The E and ML values of chromatographic paper FN-2 with 
lymphocytes are 2.2 and 16.6 times higher than without lymphocytes, resrectively. Also the E value 
on the glass is positive, while on FN-2 it is negative, what corroborates the mechano- 
electrochemical mechanism of ML generation, due to discharge events in the double electrical 
layer [33]. 

This can be supported by the experiments which show that the effect of the gamma-emitter 
(80Co) reduces cellular ML by 24%. This is due to the well established fact that the E value 
decreased as a result of ionizing radiation action [34]. 

We studied the effect of the fluorescent probe 5-phenyl-2-(thienyl)oxazole on lymphocyte 
ML. It caused a significant reduction of lymphocyte surface ML for the lymphocytes incubated 
with the probe, versus that of the lymphocytes incubated in saline (72% reduction). The ML 
inhibition may be attributed to the chemical interaction of fluorescent probe with lymphocyte 
membranes influencing the electric charge of the cells [35]. ; 

We studied the comparative changes in ML from murine lymphoid cells to experimental 
effects producing a certain effect of cellular shape deformation: at different intensities of the 
acoustic field initiated by ultrasone irradiation, the influence of additional gravitational forces due 
to centrifugation, as well as indirect influence of osmotic pressure of the NaCl solution acting on 
the plasmic membrane during varying molarity of sodium chloride solution. Following ultrasound 
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exposure in the range of 62—94 dB at the constant frequency of 3.75 MGz we observed the 
increased ML intensity of splenic lymphocytes, and what is especially significant, by 288% at 78 
dB. Various regimes of centrifugation in the range of 160—2800 g during sedimentation of splenic 
lymphocytes in the suspension resulted in changes of cellular ML relatively to that of lymphocytes 
centrifuged at 400 g, a regime is most often used in laboratory practice, only at 900 g a 35% 
decrease was observed and a 51% increase at 2800 g. An increase in molarity of NaCl solutions 
from 0.07 to 0.23 M caused a gradual 4-fold increase in splenic lymphocyte ML. Our studies 
suggest the dependence of cellular ML on putative deformation of cellular surface and shape in 
general as a consequence of artificial effects based on a mechanochemical factor. 


MECHANOLUMINESCENCE OF ELECTROLYTES AND BIOPOLYMERS 


This section describes the results of studying the ML of some electrolytes and biopolymers 
recorded with use of the techniques with electric tribogenerator [36—38]. 

It was revealed that with a rise in the concentration of NaCl, NaBr, KCl electrolytes from 
10? to 5-107 mol/l the ML intensity increased 10-fold on the average (Fig.7a). Amino acids 
glycine, phenylalanine, tyrosine and tryptophan are distinctive of their low initial ML level and 
with the rise of concentration their luminescence grows insignificantly (Fig.7b). Lipids — р 
апа oleic acids, cholesterol have high initial ML and with the concentration rise from 3- 10° 
1.5* 10! mol/l their luminescence increases 10-fold, on the average (Fig.7c). 
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Fig.7. ML: a) 1 — NaCl; 2 — NaBr; 3 — KCl; b) 1 — glucine; 2 — tyrosine; 3 — tryptophan; 4 — 
phenylalanine; c) 1 — palmitinic acid; 2 — oleic acid; 3 — cholesterol 


We have also carried out the comparative determination of ML intensity of human blood 
plasma samples and various classes lipoproteins isolated from plasma in concentrations 
corresponding to their status in human blood plasma (Fig.8). The highest ML intensity was noted 
in low-density lipoproteins and then, in descending order, follow very low-density lipoproteins and 
high density lipoproteins. It is essential that the ML intensity of each lipoprotein fraction exceeds 
that of the whole blood plasma, while the intensity of low-density lipoproteins is three times higher 
than that of blood plasma. 
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luminescence. Protein components Fig.9. The spectrum of DNA mechanoluminescence: (А), DNA, RH 75%; 
of the blood, particularly aromatic (0. ОМА, RH 92% ў 

amino acids, may play a certain role as light energy converters. Having their own low 
luminescence, these amino acids participate in emission of the energy of freeradical oxidation of 
lipids, correcting its spectral composition and intensity. Presence of antioxidants in blood plasma 
conditions the ML stabilization at the significantly lower level than luminescence of isolated 
lipoproteins, which are the part of it. Lipoproteins, isolated from plasma, deprived of antioxidant 
protection, are bound to slight oxidation with the increase of ML intensity. As for the reasons, 
which condition the enhanced ML of electrolytes, they may be primarily related to the ability of 
the latter to electrify actively, what in turn, may initiate intensive luminescence. 

In view of the fundamental role of DNA macromolecules in biological systems and their 
electric properties which are dependent on their hydration [39—49], I studied the light emission 
characteristics during mechanical activation of DNA (Fig.9) [50, 51]. The highest maximum of light 
emission for DNA with RH 7596 (A-form) was recorded at 422 nm and exceeded the ML value of 
nucleic acid at this wavelength with RH 9296 by 35-fold. Maximum light emission from DNA with 
RH 9296 (B-form) was observed at 252 nm and exceeded 6-fold the ML of nucleic acid with RH 
75% at the same wavelength. 

The main contribution to the structure of ML spectra is due to the degree of atomic 
ionization and that of molecules of the biopolymer environment, initiated by the process of charge 
recombination during mechano-electrochemical generation of the electric discharge [52—59]. The 
pattern. of ML spectra describes the interaction of electronic transitions with local atomic 
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fluctuations in the emission center. This is determined by the energy of DNA chemical bonds and 
the kinetics of electronically excited states in nitrogen and oxygen [60, 61]. This is also confirmed 
by the findings of our experiments. Since the gas discharge events during mechanical activation 
took place in the nitrogen ambient, the ML spectra of DNA, while their intensity values for light 
emission at fixed wavelengths were essentially different. 


ANALYSIS OF BIOLOGICAL MECHANOEMISSION PROCESSES 


Various existing theories which strive to explain the ML of objects of non-biological origin 
[62] could only partially interpret the biological ML mechanisms. Let us therefore consider the ML 
mechanisms of biological structures. The mechanism of ML emission is based on the streamer 
theory, according to which the weakly luminescent gatherings of ionized particle, i.e. streamers, 
arise between the oppositely charged surfaces. Streamers create electrically conductive and 
luminous bridges in the form of Lichtenberg's figures. 

Let us examine now the elementary processes underlying ML of objects of biological origin 
[63]. Mechanical activation of biopolymer results in emergence of electric field, a rise in the 
temperature, activation of chemical reactions and photon emission. All these processes accompany 
the phenomenon of tribo-induced luminescence. This mechanism may be given by the following 
equations 

505 


kT 
с & [о pu n; + не), (1) 
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Equation (1) describes the process of biopolymer electrization, where: o is the surface density 
of the charge; & is the electric constant; є is a relative dielectric penetrability; ЛФ — work 
function; k — Boltzman's constant; Т — thermodynamic temperature; e — electron charge; и; — 
particle concentration; Н — distance between the electrization surfaces; E — electrostatic field 
strength. Equation (2) describes the process of microdischarge between the dielectric barriers, 
where W; — drift velocity of electrons, ions, phctons, neutral atoms and molecules. 

Mechanical activation of bioobjects does not lead to photon emission. Excited states to which 
the laws of tribo-induced electro- and chemiluminescence are applicable from the position of 
radically-chained oxidative reactions, occur in accordance with which the optical emission value in 
the processes under study is conditioned by recombination of free radicals, as shown in equation 


1= 01,1, (3) 


where wœ — speed disproportion reaction of free radicals, п, — excitement exit, 7, — emission exit. 
The, above chemo-physical basis is used for the model analysis of ME processes in biological 
objects, particularly in DNA macromolecules and cells. It is known that a chromosome actually 
consists of two long DNA molecules wrapped about one another in the shape of a "double helix". 
As shown in Fig.10a the two strands are held together by electrostatic forces — і.е. by the 
attraction between positive charges and negative charges. We see in Fig.10b, that an A (adenine) on 
Опе strand is always. opposite to a T (thymine) on the other strand; similarly, a G (guanine) is 
always opposite to a C (cytosine) [64]. As the result of mechanical forces F, and F, action on the 
DNA structure between 2 nucleotide bases there may be changes in the electrostatic field of 
bordering regions accompanied by separation and redistribution of electrical charges, as well as 
polarization, re- or depolarization of certain structures, within which the charge can be conveyed 
by electrons or ions. The disruption of the contact between bases is accompanied by a sharp 
increase in potentials which leads to electrical microdischarges in the double electric field at 
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streamer current (0, light emission (hv), electromagnetic fields (E, Н), acoustic fields (A), 


temperature changes (Af) and which also initiates the radical chain reaction (-Н--№----:; 
~H-+O~.--). : 
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Fig.10. Mechanoemission model for DNA: H 
a) Section of a DNA double helix: A — 
adenine; C — cytosine; С — guanine; T — 
thymine; b) Chemical physics reactions 
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Our data of biological ME in the optical spectrum can be interpreted in the terms of cellular 
model of the ME of a bioreactor (Fig.11) [65—66], The theoretical model is based on the following 
well established assumptions: (1) the cytoplasm of the cell and such contents as organoids and 
colloid particles аге in constant movement and are subject to mechanical activation; (2) the cell is 
subject to external contact effects from the adjacent cells or colloid particles; (3) during transport 
events, cellular membranes are subject to local mechanochemical effects in enzymes, channels and 
pumps. 

The effect of mechanical properties, viz. protein deformation, liquid crystal parts of the 
membranes, DNA acentric structure and other cellular elements, initiate the radical chain reactions 
(RH..., A..., B...) in individual cellular elements. UCM LONE 

This is concurrent with light emission (hv) and heat (AT), electromagnetic (E, Н) and 
acoustic (4) field effects. Mechanochemical reactions, accompanied by ME, may influence the 
matrix synthesis of nucleic acids, protein synthesis in ribosomes, and modification of proteins and 
lipids in Golgi apparatus, including immunochemical reactions. 

The phenomenon of dissipation, through various channels, of cellular mechanical energy, 
apparently accounts for the generation of two major effects. The first is an exogenous emission of 
the energetic excess of cellular metabolic products to prevent,: for instance, events of excessive 

. intracellular peroxidation; the other effect is endogenous, accounting for the physicochemical 
events of collective interactions of various cellular components and concurrent influence on 
' intracellular interactions. 
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Fig.11. Cell Model of a Mechanoemission Bioreactor: Fig.12. Schematic Model of Mechanoemission Cellular Modu- 
1 — Plasmic membrane; 2 — Nucleus; 3 — Mitochondrion; lation: 1 — Universe; 2 — Earth; 3 — Organism; 4 — Organ; 
4 — Rough endoplasm; 5 — Cytoplasm; 6 — Golgi appara- 5 — Cell; 6 — Cellular nucleus; 7 — Chromatin; M — Mitosis; 
tus; 7 — Nucleolus; 8 — Nuclear membrane; 9 — Trajectory G, — Presynthetic period; S — DNA synthesis period; G5 — 
of ionizing particle (quantum) Postsynthetic period; Со — Quiescence phase 


Electromagnetic and ionizing radiation produce charged particles, ions, electrons and excited 

molecules in the substances. The transformation of primary radiolysis products results in the 
‘formation of free radicals. Such changes lead to alterations in the parameters of the membrane 
structure and surface, what on the whole, might influence postirradiation release of ME 
characteristic, ML parameters in particular. 

Cellular function is subject to appreciable ME modulation in living systems, under the 
immediate influences of the environment via endogenous and exogenous regulation factors. Let us 
conventionally subdivide them into seven regulatory steps (schematically shown in Fig.12). 

ME cellular modulation takes place at: (1) sub- and supra-molecular levels because of the 
mechanochemical properties of biopolymers; (2) the level of the cell nucleus — information 
channels for this were considered above; (3) the cell level during the cell cycle — cellular 
differentiation at the molecular level consists in the regulation of gene functioning — channels of 
communication at this level utilize plasmatic membranes and the cell's cytoskeleton; (4) the organ 
level — information is transmitted by cellular interactions involving both contact and non-contact; 
(5) the organism level — information is transmitted via the circulatory system, lymphostasis, central 
nervous system and meridians; (6) at the planetary level — the Earth — via direct contact with the 
habitans and non-contact communication through the electromagnetic and gravitational currents of 
the Earth; (7) the level of the Universe via non-contact channels — probably electromagnetic and 
gravitational waves. 


CONCLUSION 


Certainly, this review of biological ML could nog provide a generalized picture of all 
experimental and theoretical aspects of this new approach to the study of the MF phenomenon in 
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biological objects. No doubt, further laboratory studies will provide a deeper insight into the 
mechanisms of enzymatic and free radical reactions in biological membranes as well as into the 
genetic apparatus of cells, hemostasis, etc. At the same time, the use of mechanoluminometry in 
clinical laboratories allows to make diagnoses of free-radical pathologies in oncology, radiology, 
pediatrics, cardiology, toxicology and some other areas of medicine [67— 73]. 


БІОЛОГІЧНА МЕХАНОЛЮМІНЕСЦЕНЦІЯ 


B.E.OPEJI 

B статті розглянуто проблему біологічної механолюмінесценції — оптичного випромінюван- 
ня, що виникає в результаті механоелектрохімічної активації біологічних об'єктів. Описано 
принципи  експериментальної техніки, ініціюючої механолюмінесценцію біологічних 
об'єктів з використанням електричного, оптичного та гідродинамічного трибогенераторів. 
Приведено та проаналізовано експериментальні вимірюваання механолюмінесценції рідкої 
крові, личинок, фаланги пальців, судин та лімфоїдних клітин. Звернено увагу на залежність 
інтенсивності механолюмінесценції електролітів, амінокислот та ліпідів від їх концентрації. 
Продемонстровано спектри механолюмінесценції ДНК в А- і В-формах. Запропоновано 


біофізичні моделі механоемісійних процесів. 
БИОЛОГИЧЕСКАЯ МЕХАНОЛЮМИНЕСЦЕНЦИЯ 


В.З.ОРЕЛ 

В статье дан обзор проблемы биологической механолюминесценции — оптического 
излучения, возникающего в результате механоэлектрической активации биологических 
объектов. Описаны принципы экспериментальной техники, инициирующей механолюми- 
несценцию биологических объектов с использоваанием электрического, оптического и 
гидродинамического трибогенераторов. Приведены и проанализированы эксперимен- 
тальные измерения механолюминесценции жидкой крови, личинок, фаланги пальцев, 
сосудов и лимфоидных клеток. Обращено внимание на зависимость интенсивности 
механолюминесценции электролитов, аминокислот и липидов OT MX концентрации. 
Продемонстрированы спектры механолюминесценции ДНК вА- и В-формах. Предложены 
биофизические модели механоэмиссионных процессов. 
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Deterministic chaos has been proposed as an underlying conception fo: analyzing calcium 
oscillating signals. The possibility of chaotic behaviour has been discussed for the dynamics of the 
calcium oscillations in different cells. This description is based on the theory of dynamical 
systems. In order to consider possible deterministic behavior of physiological systems, we describe 
some types of attractors by using simple examples. In an attempt to differ chaos from stochastic 
noise, we investigate the characteristic features of deterministic chaos in intracellular calcium 
oscillations using return maps. It is shown that irregular oscillating calcium signals can be 
characterized in terms of the low-dimensional deterministic chaos. On the basis of the approach 
developed, it was analyzed some experimental data of calcium oscillations in secretory cells 
presented in literature. This data allow to make the conclusion that the irregularity in the 
oscillating calcium signals in pancreatic acinar cells is caused by appearance of deterministic 
chaos at low [ACh]. At high [ACh] this irregularity is also connected with the stochastic noise. 
The influence of the another agonist - vasopressin on the calcium signal in hepatic cell is 
characterized.by appearance of deterministic chaos only. In the case of gonadotropin-releasing 
hormone for gonadotrophs and of carbachol for rat parotid cells the irregularity is caused by 
deterministic chaos and random noise, respectively. 


( Received August 6, 1994) 


INTRODUCTION 


Deterministic chaos has been proposed as an underlying principle of organization in many 
systems, including biological ones [1, 2]. The phenomenon of deterministic chaos was discovered 
by Lorentz [3]. The mathematical foundations of the analysis of the chaotic states existence in 
dynamical systems were developed by Shilnikov and Mel’nikov in the 1960s [4—7]. In spite of the 
many different definitions of low-dimensional deterministic chaos, it can be refered to aperiodic 
behavior of a dynamical system. Deterministic chaos is characterized by high sensitivity to initial 
conditions and is observed under strictly defined experimental conditions. The time dependence of 
system parameters (e.g., concentrations of intracellular calcium) is like the stochastic behavior due 
to the noise. Chaotic states usually arise from periodic or complex-periodic states when the 
experimental conditions (controlling parameter) change. 

In physiology the deterministic chaotic processes have been put into the sphere of scientific 
interests only in the last years [8]. The deterministic chaos can be identified in information transfer 
of the endocrine regulation. Parathyroid hormone secretion seems to be the first example of a 
chaotic hormonal rhythm in human physiology [9]. 

The aperiodic oscillations of calcium signal had been found during 1980s [10—11]. Many 
investigators discovered that neurotransmitters or hormones could evoke calcium oscillations in a 
number of different exocrine glands [12], fibroblasts [13], liver cells [14], endothelial cells [15], 
HeLa cells [16], macrophages [17], oocytes [18], megacariocytes [19], neurones [20]. Reviews of 
the last achivements in this subject are given in [10, 11, 21]. 

In the resting cells cytosolic free calcium concentration ([Ca?']) is tightly controlled to the 
levels of 50—100 nM. The calcium concentration outside the cell is approximately 10* times higher 
than that inside the cell. The control of [Ca?^'], in resting and stimulating cells is achieved by 
multiple chanells, pumps, and exchangers. [Ca]; is ultimately controled by the proteins in the. 
plasma membrane. Calcium enters throuhg channels that are gated by voltage, neurotransmitters, 
G proteins, and cytosolic messengers. The Са" -АТРазе pump and the sodium calcium exchanger 
extrude Са"! from cells. | Са: |, is also continuously regulated by fluxes between the cytosol and the 
internal calcium stores. Intracellular calcium Tunctions as an important intracellular signaling 
mechanism whereby many agonists regulate different cellular processes [22]. 
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DETERMINISTIC CHAOS IN CYTOSOLIC CALCIUM OSCILLATIONS 


Hormones, neurotransmitters, and growth factors activate cellular responses by mobilizing 
the cellular Ca?'. In general, this Ca mobilization response is composed of release of Ca^ from 
intracellular stores, as well as of increased entry of Ca * into the cytoplasme from the extracellular 
space. Changes in [Са?*],; are moderated by an extensive system of endogeneous Ca * buffers. Ca?* 
reuptake into intracellular stores and Ca * extrusion across the plasma membrane prevents an 
excessive rise in [Ca?*]; and could play a role in the termination of cytoplasmic Са?" signals. These 
signals that last more than a few minutes after the onset of the stimulation depend critically on the 
influx of extracellular Са" [23]. Such calcium influx can be triggered in many types of cells by 
depletion of the intracellular calcium stores [24]. The mechanism by which store depletion can 
control plasma membrane Ca?* permeability remaines contoversial. It was shown that the emptying 
of intracellular Са?" stores releases a novel small messenger (calcium influx factor) which stimulates 
Ca influx [25]". 

External signals arriving at the cell engage surface receptors to initiate signalling pathways 
whereby the information flows from one component to the next until the final effector system is 
activated. The formation of inositol(1,4,5)-trisphosphate (InsP3) is the focal point for two major 
pathways, one initiated by a family of G protein-linked receptors and the other by receptors linked 
by tyrosine kinases either directly or indirectly [26]. These separate receptor mechanisms are 
coupled to the energy-requiring (GTP or ATP) transducing mechanisms which activate 
phospholipase C (PLC) to hydrolyse the lipid precursor phosphatidyl inositol (4,5)-bisphosphate to 
give both diacylglycerol (DAG) and InsP;. DAG activates protein kinase C (PKC). InsP; binds to 
an InsP3-receptor (IP4R) to mobilize stored calcium and to promote an. influx of external calcium, 
perhaps, working in conjunction with InsP, [26]. 

InsP; is (һе second messenger that controls many cellular processes by generating internal 
calcium signals. A Ca release from intracellular stores is mediated via two separate, but related, 
types of channels: the ryanodine receptors (RYRs) and the inositol(1,4,5)-trisphosphate receptors 
(ТР.Вз). These two types of receptors will be considered together because they share considerable 
structural and functional homologies [26]. The remarkable structural similarity that exists between 
IP;Rs and RYRs probably reflects a common evolutionary origin. Upon binding InsP; the IP;R 
undergoes a large conformational change which is perhaps related to the coupling process leading 
to channel opening. RYR can be activated by Са" and cyclic ADP-ribose or by depolarization via 
the direct link from the plasma membrane dihydropyridine receptor [26]. The plant alkaloid 
ryanodine opens the channels at low (nanomolar) concentrations but closes them at higher doses 
(micromolar). Caffeine can also open the RYR channels. These RYRs contribute to calcium 
signalling in many different cell types. Recent evidence suggests that the sensitivities of IP4R to IP; 
and RYR to Ca^' can be modulated. IP;R and RYR-controled Са?" stores vary in their number’ 
and intracellular distribution in different cells [26]. 

This family RYR of intracellular calcium channels displays the regenerative process of 
calcium-induced calcium release. Such a signalling pathway controls many cellular processes [10- 
11]. The mechanism of agonist-evoked Ca oscillations has been extensively discussed and a number 
of very different hypothesis are the subject of intense debate [11]. 

The possibility of chaotic behaviour has been discussed for the dynamics of the electrical 
activity of the human heart [27] and brain [8] on the basis of the theory of dynamical systems and 
quantitative parameters of deterministic chaos. In this paper we use the following terms: 
deterministic chaos and aperiodic (irregular) motion. The term ‘aperiodic’ motion (regime) is used 
for the description of any irregullar (aperiodic) oscillations of calcium signal. This aperiodicity can 
be caused by different intrinsic effects: transient character of process, quasiperiodicity, deterministic 

chaos, random motion [5—6]. We say that aperiodic signal is caused by deterministic chaotic 
processes only when this signal is analyzed with the use of some algorithmes for verification of 


*Calcium influx can be triggered in some types of cells by inositol phosphates [26]. 
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deterministic chaos [5—6, 28—30]. Today the aperiodicity was found in many oscillating 
physiological processes, but, only in some cases the deterministic chaos is verified [8, 9]. 


SOME PRINCIPLES OF THE THEORY OF DYNAMIC SYSTEMS 


The main successes in the analysis of the properties of chaotic states and determination of 
their distinguishing characteristics derive from the mathematical analysis of dynamical Systems. 
Quite extensive material has been accumulated in this field. This material is presented in a number 
of books and reviews [4—6, 30—35]. 

The object of this paper is to describe the characteristic features of deterministic chaos in 
cytosolic calcium oscillations. We will attempt to differ chaos from other irregular oscillations (e.g., 
stochastic noise). 

The concept of chaotic oscillations in physiological systems became possible due to the 
development of the theory of dynamical systems that works with such concepts as limit cycle, 
*strange" attractor, etc. However, they have not been widely employed in physiology, and for this 
reason we believe it is necessary to define the basic concepts of the theory of dynamical systems 
that are employed for describing chaotic physiological processes [32]. 

Attractor - trajectory of the system in phase space after transient processes are completed. 

Bifurcation - qualitative change in the character of the trajectory of the system. 

Bifurcation parameter - controlling parameter; when this parameter reaches a definite value 
(bifurcation point), the character of the trajectory of the system changes qualitatively (for example, 
regular behavior of the system is replaced by chaotic behavior). 

Bistability — existence of two stable stationary points under identical conditions. 

Deterministic chaos - irregular, unpredictable behavior of a dynamical system, observed in a 
definite range of values of system parameters, which is characterized by high sensitivity to the initial 
conditions. 

Dissipative ИСЕ - stable and spatially nonuniform structures arising as a result об 
development of instabilities in a uniform nonequilibrium dissipative medium. 

Dissipative system - a dynamical system whose state variables are being constantly consumed. 

Dynamical system - a system whose evolution is determined uniquely by its initial state. 

Fractal dimension - a quantitative characteristic of an attractor showing how completely the 
attractor fills the phase space (a more detailed discussion is given in [4—7]). 

Largest Lyapunov exponent (4) - for a system, consisting of N independent variables 
(Ci, С,..., Су), is defined by the following expression: 

АЦ!) 


1 
A= li lim -1 1 
Adone "E AL(0) (0 


N A 

where ДЦ) = УСС Cy) - СС Со + 8Co у, | is the distance between neighboring trajecto- 
i=l 

ries of a system in the phase space of the system, C; (1, C) is the value of the i-th variable at time / 

under the condition that its initial value was Су, and SC) is an infinitesimal variation of the value 

(8Q,/C << 1). 

Hopf bifurcation — bifurcation of the formation of a limit cycle from a point, characterized 
by a transition of the real part of the complex-conjugate pair of eigenvalues through zero at the 
bifurcation point. | 

Intermittency — a type of chaotic behavior in which long time intervals of regular motion are 
replaced by bursts of irregular behavior. 

Limit cycle — attractor, comprising a closed trajectory in the phase space of the system. 
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Multistability — the existence of more than two stable stationary points under identical 
conditions. 

Phase portrait — trajectory of the system in phase space. Phase space - space whose 
coordinates are the values of the quantities characterizing the behavior of a dynamical system. 

Random noise - internal and/or external fluctuations of the medium in which the dynamical 
system evolves. 

Return map - plot of the value of a successive maximum of a measured parameter (for 
amplitude — next amplitude map) as a function of the preceding value. 

Sustained oscillations - undamped oscillations in a dissipative nonlinear system which are 
supported by influx of reagents. 

Strange attractor – attractor which is neither a point nor a limit cycle. 

Trajectory - line depicting the change of the variable or vector of the system variables in 
time. 

Transient chaos - a deterministic chaotic state, which is observed for an oscillatory system 
when it evolves in time from the initial state into the position of the end state. 

In the most general case the evolution of a physiological system, neglecting spatial effects 
and fluctuations, is described by a system of ordinary differential equations [29]: 


dC, f dt e FAC Оки), E93 s, (2) 


where C; is the value of the i-th variable of the system, и is the controlling parameter (the 
concentration of stimulator, temperature, etc.), and М is the dimension of the System. The 
properties of this system cam be analyzed by several methods, in particular, in terms of the theory 
of the phase space of dynamical system. The dynamical approach is most convenient if the phase 
space of the system is used for its description. It turns out that the diverse time dependencies, 
observed for stationary regimes in phase space, correspond to a very small number of qualitatively 
different forms of phase portraits (for example, a point corresponds to a stationary state, a limit 
cycle corresponds to regular oscillations). The most important concept employed in the analysis of 
the properties of the phase portrait of a system is the concept of stability. It is intuitively obvious 
that a system is stable, if it is capable of restoring the character of its motion under insignificant 
external disturbances. In application to the systems described by Eq.(2) this means that the phase 
portrait of the system with such disturbances does not change qualitatively. In the case when a 
qualitative change occurs in the phase portrait of the system accompanying a change in some 
parameter of the system (ie. the character of the time dependency of intracellular calcium 
concentration), it is said that a bifurcation has occurred. Thus the existence of a bifurcation is due 
to a change in the stability of different states and/or appearance of new states. 

In order to consider the possible deterministic behavior of physiological systems, we describe 
some types of attractors by using simple examples. If the behaviour of a system can be described in 
terms of the rate of consumption or accumulation of a single variable, then the dimension of the 
phase portrait of such a system is 1, and only one type of attractor — a point — is possible. In the 
case when two different variables are present, two types of attractors — a point and a limil cycle can 
be realized in the phase space of the system. The limit cycle corresponds to regular temporal 
oscillations of the system. For example, in [36] it is shown that a simple two-dimensional model 
has a limit cycle and can describe hormone-induced calcium oscillations: in liver cells. In the 
general case, oscillations in a physiological system can be only observed far from thermodynamic 
equilibrium if some autocatalytic processes occur in the system, i.e. when some substance 
accelerate its own production [31]. If the system is bistable, then so-called relaxational sustained 
oscillations will be realized in it. Such sustained oscillations correspond to rapid transition of the 
system from one stationary state into another, after which the system slowly relaxes back. The 
behavior of such a system can be described with the help of a cross-shaped bifurcation diagram, 
which is characterized by regions of stabilit$ of one stationary state (SSI) and of the other stationary 
state (SSID, region of bistability (BS), and region of sustained oscillations (RO) (Fig.1) [32]. 
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Fig.2. Bifurcation diagram of the logistics map 


Hy 


Fig. І. Cross-shaped bifurcation diagram illustrating behavior of the system becomes more complicated as 


It follows from the above discussion that the 
asians a и. зако Gs Е oe eia the number of its components increases. When the 
parameters) number of variables is equal to 3, there appears a new 
type of motion, which in the phase space of the system does not correspond to a point or a limit 
cycle. This motion is caused by the appearence of a so-called strange attractor. In accordance with 
the hypothesis of Ruelle and Takens [4] a strange attractor corresponds to the presence of 
deterministic chaos in the temporal behavior of the system. The concept of deterministic chaos is 
inherently connected with the concept of order. Thus, regular sustained oscillations give rise to 
definite order іп the system - periodicity. Perhaps, in physiology, order and chaos are related with 
one another by complicated relations: ordered regimes alternate with chaotic regimes. Such 
intermittency is observed, for example, in many chemical reactions which can be considered as the 
simplest models of a living systems [37—39]. 

When irregular oscillations are observed in an experiment, there arises the question: what 
type of sustained oscillations should they be referred to? In literature, often without additional 
investigation, these oscillations are called irregular, complicated, composite, unusual, chaotic, etc. 
In this case it is reasonable to determine such an oscillations as irregular oscillations, which can be 
of different nature (for example, chaotic and stochastic). Additional investigations must be 
performed in order to determine whether such oscillations are chaotic (i.e., deterministic) or 
stochastic (i.e., owing to random noise). 

Among experimentators the most popular quantitative characteristics of deterministic chaotic 
oscillations are the positive value of the largest Lyapunov exponent (LLE) and the fractional value 
of the fractal dimension. Efficient algorithms for calculating both these and a number of the other 
parameters, with the help of which the observed chaotic regimes can be described quantitatively, 
have been developed [28—30]. Such quantitative description has been done for cardiac oscillations 
[8, 9]. 

Scenarios of the deterministic chaos appearance, i.e., the sequence of bifurcations of 
different regimes, as a result of which a chaotic process develops, can serve as a qualitative 
criterion for deterministic chaos. In particular, an unusually simple path to chaos, called 
Feigenbaum's sequence, was discovered for the logistics gquation (Fig.2) [4—6] 


X, = Hk (= 3): (3) 
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Fig.3. The return maps of the irregular oscillations (the data are taken from [12]) showing the effect of increasing 
[ACh] on the calcium signal (a, b, c - the same pancreatic acinar cell, d – another pancreatic acinar cell) 
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Fig.4. The return maps of the calcium spiking in a single rat hepatocyte. The cell was stimulated by perfusion 
with vasopressin (the data are taken from [14, 21]) 


The law discovered by Feigenbaum is applicable to any system whose behavior is 
characterized by a very general property: in a definite range of values of the parameters a system 
operates in the periodic state with period T; on passing through the threshold the period doubles 
and becomes 2T; on passing through the next threshold the period doubles again and is now 4; etc. 
Thus the system is characterized by a sequence of period doubling bifurcations which, finally, lead 
to the appearance of deterministic chaos. The quantitative relations between the values of the 
bifurcation parameters are strictly determined [4—6]. We can notice that the description of chaotic 
behaviour by using the maps like those in (3) is very simple and fruitful. Such a map can be 
obtained from any experimental data, for example, as a function of the value of a successive 
maximum of a measured parameter from the previous one. It appears that for many “typical” 
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chaotic regimes this next amplitude map is a function with single maxima [4]. Other reasons for the 
verification of chaotic oscillations: 


I  — Irregular behavior is observed with strictly determinate values of the controlling parameters 
and initial conditions. 

П  — Typical chaotic regimes are characterized by a definite type of attractor (Rossler, Lorenz, 
etc.) and return map (with one maximum). 

Ш — Deterministic chaos arises as a result of a definite sequence of bifurcations: 1) period 
doubling; 2) period addition; 3) destruction of a quasiperiodic regime; 4) intermittency. 

IV. — Deterministic chaos is characterized by a positive value of the largest Lyapunov exponent 
and fractional fractal dimension. 

У — Deterministic chaos manifests high sensitivity to the initial conditions. 


AN APPROACH TO INVESTIGATE IRREGULAR OSCILLATIONS OF INTRACELLULAR 
CALCIUM SIGNALS 


The attractive goal of the possible application of the deterministic chaos theory lies in the 
possibility to investigate an irregular behaviour of the oscillating intracellular calcium signals. These 
signals control many cytosolic processes, including fertilization, cell growth, transformation, 
secretion, muscle contraction, sensory perception, and neural signalling [26]. The number of 
reports about these oscillations is increasing, and the oscillations of [Ca^] have been recognized as 
a general phenomenon by which cells transduce the extracellular signals to intracellular effector 
system. Now we describe the results of the investigation of whether or not irregular oscillating 
calcium signals can be characterized in terms of the low-dimensional deterministic chaos. 

Let us consider the possible application of using the return maps for the analysis of aperiodic 
calcium oscillations in secretory cells [12, 14]. The inset oscillations are presented in Fig.3—6. 
These data shows the time dependencies of the intracellular calcium spikes for different agonists at 
their different concentration. It follows from data presented in Fig.3—6 that in all cases the 
oscillating process has transient character. The amplitude of these oscillations is changed in time, 
but the time interval between adjoined maxima ("period") remains unchanged. These calcium 
oscillations have an irregular character due to the aperiodicity of the amplitude. The quantitative 
description of the data presented on Fig.3—6 is impossible by using the reconstruction of phase 
space, calculation of the largest Lyapunov exponent or fractal dimension, due to the small number 
of experimental points (it is needed 1000— 10000 points for such an analysis). 

We consider the possible approach to verification deterministic chaotic character of calcium 
signals by using the next amplitude map. In order to construct such maps, let us introduce the 
following definitions: 

E; - i-th amplitude of calcium signal; 

T; - transient component of E; 

$; - aperiodic component of Е. 

The following relationship between these values exists: 

E =T, +8; 

In the attempt to take into account the transient character of the irregularity of calcium 
oscillations, one can find the aperiodic component by the standart trend analysis in the following 
manner: 


Ті; = СЕ ЕЕ) 3 
5; e Ti 


This procedure leads to a decrease in the number of points, i.e. if calcium signal has М 
maxima then next amplitude map has (N-2) points. The maps reconstructed from the values of 5 
are shown in the Fig.3—6. The analysis of these maps can be based on the notion that deterministic 
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chaos is characterized by a smooth single-peak return map, whereas stochastic noise is 
characterized by a random map like the "sky with stars". 

Fig.3 shows the effect of increasing of the acetylcholine concentration ([ACh]) on the 
calcium signal in pancreatic acinar cells. Within a narrow range of [ACh] (Fig.3.a, b, c) the effect 
of increasing agonist concentration is to convert the response from chaotic ocsillations, which are 
characterized by single-peak return map (see Fig.3.a, b), to practically stochastic noise (see 
Fig.3.c). The character of signal at higher [ACh] for another cell (Fig.3.d) can be estimated as a 
mixture of deterministic chaos and noise. These data allow to make a decision that irregularity in 
the oscillating calcium signals is caused by appearance of deterministic chaos at low [ACh]. At high 
[ACh] this irregularity is also connected with the stochastic noise. 

The influence of the another agonist - vasopressin (VP) on the calcium signal іп hepatic cells 
is presented in Fig.4 for different concentrations [VP]. In this case the irregular part of calcium 
ocsillations can be characterized as a deterministic chaos for each [VP]. The differences between 
these oscillations at different values of [VP] can be described in terms of quantitative parameters of 
chaos only (like the largest Lyapunov exponent and fractal measure). 

The other kinds of return maps are presented in Fig.5 and Fig.6 for carbachol (CCh) and 
gonadotropin-releasing hormone (GnRH). These two maps have the similar appearance, but the 
‚ amplitude of the calcium oscillations is higher for GnRH with respect to the ССВ. The comparison 
of these two maps with that presented in Fig.3-4 can give the idea that in this case the irregularity 
is caused by deterministic chaos and random noise in the case of GnRH. For CCh the irregularity 
is caused by the random noise only. | 


GnRII inm 


10 CCh 10ит 


5 + Ман А 


205 


-50 -30 -10 10 30 50 
5» nM У» nM 
Fig.5. The return map of the irregular oscillations (the data | Fig.6. The return map of the irregular oscillations (the data are 


are taken from [41]) showing the effect of carbachol on the taken from [42]) showing then effect of gonadotropin-releasing 
calcium signal in rat parotid gland cell `. hormone on the calcium signal in gonadotrophs 


So, in many cases the aperiodic calcium oscillations are verified on the basis of the approach 
developed as low-dimensional deterministic chaotic process. The proposed approach will be a 
usefull instrument in the analysis of cytosolic calcium signals. 


CONCLUSIONS 


In conclusion we discuss the basic problems arising during the study of irregular oscillations 
in physiological systems. Probably the most important one is the problem of the difference in the 
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deterministic chaos and stochastic noise in these systems starting from experimental data. In order 
to solve this problem simple predictive criteria must be developed for the appearance of 
deterministic chaos and the methods of quantitative analysis of irregular oscillations must be further 
elaborated. The simplest and most promising path is the search for and classification of sequences 
of bifurcation of the appearance of deterministic chaos. Moreover, the cell is a spatially organized 
system, so there appears the problem of experimental investigation and quantitative description of 
spatiotemporal chaos. Another problem related to calcium oscillating signals is connected with the 
information transfer in living cells. The deterministic chaotic signal brings much more information 
than the periodic one [40]. For example, chaotic systems can presumably allow optimal adaptation 
to fluctuating environmental conditions. It has been argued recently that chaotic signals can be 
controlled by applying judiciously chosen minute perturbations precisely because chaos is highly 
sensitive to initial conditions. Perhaps, the key to the problem “what is an intrinsic aim of chaotic 
calcium oscillations in living cells?" can be obtained in terms of the information production by 
living cells. 
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ДЕТЕРМІНОВАНИЙ ХАОС Y ЦИТОЗОЛЬНИХ КАЛЬЦІЄВИХ ОСЦИЛЛЯЦІЯХ 


П.С.СТРИЖАК, І.С.МАГУРА, К.Б.ЯЦИМИРСЬКИЙ 

З метою аналіза можливості прояви детермінованого хаосу в фізіологічних системах 
запропоновано підхід щодо відокремлення хаотичної складової нерегулярних осцилляцій 
концентрації внутріклітинного кальцію. Цей метод базується на аналізі властивостей 
зворотних мап, які сконструйовано, виходячи з деяких експериментальних результатів, які 
описано в літературі для секреторних клітин. Доведено, що зміна концентрації ацетил- 
холину веде до зміни детермінованого хаотичного режиму на стохастичний шум. Інший 
агонист - вазопресин - завжди веде до появи детермінованого хаосу у цитозольних 
кальцієвых осцилляціях. Карбахол та гонадотропін призводять до появи стохастичного 
шуму в системі. 


ДЕТЕРМИНИРОВАННЫЙ ХАОС B ЦИТОЗОЛЬНЫХ КАЛЬЦИЕВЫХ ОСЦИЛЛЯЦИЯХ 


П.Е.СТРИЖАК, И.С.МАГУРА, К.Б.ЯЦИМИРСКИЙ 


С целью анализа возможности  детерминированного хаотического поведения в 
физиологических системах предложен подход, позволяющий выделить хаотическую 
составляющую нерегулярных осцилляций концентрации внутриклеточного кальция. Такой 
метод основан на анализе свойств карт возврата, которые строились исходя из значений 
амплитуд кальциевых осцилляций. Такой подход позволяет отличить детерминированное 
поведение от стохастического шума. Возможности предложенного метода проиллюстриро- 
ваны на примере анализа некоторых экспериментальных данных, описанных в литературе 
для секреторных клетках. В частности показано, что изменение концентрации ацетилхоли- 
на ведет к смене детерминированного хаотического режима на стохастическийфшум. 
Введение другого агониста — вазопрессина - всегда вызывает появление детерминирован- 
ного хаоса во внутриклеточных кальциевых осцилляциях, в то время как карбахол и 
гонадотропин вызывают появление стохастического шума в системе. 
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MEDICAL AND BIOLOGICAL 
APPLICATIONS 


Human Skin Millimeter Wave Range Reflectometry Application for Diagnosis 
of Some Diseases 


I.A.IVANCHENKO, V.G.LIZOGUB, L.V.SVESHNIKOVA, YU.V.CHOVNIUK 
Scientific Research Center "Vidhuk", Kiev, Ukraine 


We observed some groups of patients by the method we proposed. The method consists in 
determination of the human skin local dielectric properties in the millimeter wave (MMW) range. 
It was performed through the measurement of the MMW differential reflectivity of the skin 
surface in a number of reflexogenic: zones. The results of investigation of the patients with heart 
rhythm disturbances showed that the frequency characteristics of AF max varied depending on the 
group of patients. Thus the technique we suggested not only can help to distinguish a healthy 
person from unhealthy one, but also to differentiate the patients between the groups of diseases. 


(Received June 6, 1994) 


INTRODUCTION 


Identification of diseases at their early stages is a problem of primary importance. Now 
various researchers conduct intensive investigations with aim of finding physical and technical 
means that can help to solve this problem. | 

It was shown in some works devoted to the MMW action on the human organism that the 
ММУ exerts substantial influence on the disturbed immune system, producing a therapeutic effect 
[1]. It appears thus that, on the contrary, when the immune system disturbances’are present, then 
some variations in the organism's characteristic are to be observed, which should be measured in 
MMW range. 

We worked out the MMW Differential Reflectometer [2] which in some cases allowed to 
determine latent pathologies by means of the differential reflectivity (DR) measurement in the 
acupuncture zones. 

In accordance with aforesaid, it can be assumed that the MMW Differential Reflectometer 
can be of essential use to the development of the early diagnosis. 

The method of differential reflectivity (DR) measurement in MMW band was described 
earlier [2]. It was shown that the sites with significant anisotropy were revealed on patients skin 
surface. The anisotropy was changing in time according to the functional state of the organism. 
During investigation of patients with neurology and gastroenterology [3], anomalous di€lectric 
properties of skin most frequently manifested themselves in the zones of “functional significance” 
according to the concepts of Oriental medicine. 

The equipment we elaborated [4] determined an indicatrix of the anisotropic tissue 
reflection, which was an ellipse with a ratio between long and short axes of A. The value of A was 
determined by : 

A = 10%, 


Where exponent AF, or anisotropy factor, was found from a pair of the measured data DR,,,, and DR min: 
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AF = DR, - DR, (in dB) 

The results of our investigations of skin reflection characteristics [4] showed the suggested 
methods to be accurate enough to detect anisotropic structures on skin surface, to determine their 
localization and to differentiate them in magnitude and extent. It provided the reliability of the 
obtained data on characteristics of the anisotropic structures on skin surface. 

The existence of the observed anisotropic structures was confirmed for a variety of 
pathologies and for a great number of patients (over 100). Thus, the method enabled us to discover 
anisotropic structures on skin surface and to determine their characteristics: spatial, temporal and 
frequency response. Anisotropy factor measurements along the acupuncture meridian of 
pericardium allow to establish the criterion of 0.5 dB (1296) in distinguishing between ill and 
healthy organisms. 

Correlation between the obtained values of parameters and some characteristics of 
pathological processes was found. Eventual change of DR profile was noted. The experiments also 
confirmed the fact that “Sources” and “Luo-points” located on the corresponding acupuncture 
channels offer the most substantial information on different pathologies [4]. 

We could draw from the experimental data the conclusion that the structures on the skin 
surface which manifested themselves as the areas with anisotropic properties were not directly 
related to morphological indications. These anisotropic properties were connected with the 
pathological processes in human organism. We suggested a mechanism of this connection. 

It is known that different pathological. processes result, on the one hand, in general 
metabolic disturbances, and, on the other hand, they are followed by local trophic and sensitivity 
changes in skin. Difficulties in determination of the above changes are due to the lack of adequate 
and sensitive methods of metabolic disturbance estimation. The ability of pathological processes to 
cause changes in the biochemical composition of blood and lymph in local sites may bring about 
changes in biophysical parameters of these liquids, for example in complex dielectric permeability. 
Owing to these changes the dependence of this value on chosen orientation, i.e., the phenomenon 
of anisotropy may appear. 

In a recently published work [5] the data obtained from the measurement of tensor 
components of the dielectric permeability and. conductivity of human skin, which confirm the 
anisotropic nature of these parameters, were presented. 

We revealed the frequency dependence of the skin anisotropy factor (AF) while studying 
dielectric properties of human skin measuring reflection of MMW radiation [4]. This experimental 
result is of special interest in connection with the known works by the other authors [6-7], which 
described the resonance effects of MMW radiation on various biological objects. It is certain that 
the line includes the microwave resonance therapy [8], in which a biological effect manifests itself 
on the level of system or organism, so as in: sensor reactions, alpha-rhythm suppression tests, etc. 
Observation of resonance phenomena directly in the area of effect brought us to the idea of using 
MMW range reflectivity measurement of skin anisotropy factor as a direct method of diagnostic 
testing. This study was aimed to reveal the relation between AF frequency response function (FRF) 
type and some types of pathologies. 


INVESTIGATION METHODS 


Experimental technique, which was used in this study, was described in the previous work 
[4]. DR values were measured on the skin surface by means of MMW-reflectometer (Fig.1), which 
was equipped with a two-channel radiation source. The radiation source consisted of two 
rectangular waveguide antennas, whose polarization planes were mutually perpendicular, thus 
providing the determination of the anisotropy factor AF on a skin in a chosen frequency range and 
obtainment of the FRF spectrogram. This technique was elaborated to make a distinction between 
the healthy people and the cardiac patients [4]. The relative error of reflectivity measurement at a 
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fixed frequency was 5%. The total error, within limits of which the surface was regarded as 
anisotropic, was characterized by 0.3 dB (7% deviation in AF values). : 

After the relevant apparatus preparation and calibration were done, the radiation source was. 
placed on the chosen site of the surface under study at the chosen initial angle of the antenna pair 
orientation axis z (Fig.1). The first record-of the ЕКЕ spectrogram was made in this position. Then 
the spectrograms were recorded at different orientation angles with the step of 45 degrees in a 180 
degrees range. Thus, a set of the DR FRF was recorded for each zone at different orientation 
angles. 


Fig.l. Measuring device of differential reflectivity (DR): a) — block diagram: GEN - generator, IND - 
indicator, MEAS - measuring circuit, SOUR - radiation source, 1 and 2 - rectangular waveguide antennas; 
b) - sections view through antennas: Ij, I, - incident waves, R;, В, - reflected waves; c) - configuration of 
antenna pair ends. The arrow Z shows antenna pair axis direction 


A single DR FRF curve within a chosen frequency range was the spectrogram with maximum 
and minimum points (Fig.2). The frequency difference between two adjacent peaks was equal to 
0.3 GHz. The figure clearly demonstrated that the difference between single FRF curves of the set 
is most distinct at the top points of the peaks. AF value was determined for each peak as a 
difference between the maximum (triangle) and the minimum (cross) DR values at the peak in dB. 
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We carried out two clinical 
investigation runs at the 12th 
Hospital in Kiev. The aim of the 
lst run was to obtain and analyze 
the spectrogram set, to discover the 
possible symptoms and to determi- 
ne their probabilities for every 
group. The aim of the 2nd run was 
to verify the diagnosis probabilities 
in accordance with conclusions of 
the Ist run. 

In addition, we studied the 
temporal characteristics of AF in 
the 2nd run. 

57 58 59 60. — 61 The first run was conducted 

Frequency, GHz on 34 patients with various patho- 
logies and on 2 subjects of practi- 
cally healthy students of the Medi- 
cal Institute. The distribution of patients according to pathologies is shown in Tab.1. 


Table 1 
Patient distribution over pathology groups 


Differential reflectivity, dB 


Fig.2. A typical DR spectrogram measured in the patient's acupuncture zone 


Total over | Had values 
group AF,20.5 dB 


Inflamm.proc. (IP 


Total 


In the group of 9 patients with coronary disease (CD) 6 patients had stable angina pectoris of 
II-IV functional class. 4 of them had attacks of angina on exertion, and 2 of them had attacks of 
angina at rest. 3 patients had acute myocardial infarctions (the investigation was carried out when 
patients were in subacute stage of Myocardial infarction). 3 of the patients had the typical for 
Myocardial infarction dynamics of electrocardiogram and an increase in enzymes concentration in 
blood plasma. Diagnostic cycle-energometry was conducted to all ischemic patients with angina on 
exertion. Tolerance for physical effort of all the 4 patients was reliably (P < 0.05) less than normal · 
relatively to their age and amounted to 87.6+5.9 W. Zones of hypo- and akinesia were detected by 
echocardiography among 3 patients. More than 1 mm depression of ST segment was observed 
among 4 CD patients. 

The CD group of patients also included 4 patients with hypertension disease (the age of 
patients varied from 25 to 70 years). Systolic arterial blood pressure was 168.4+6.3 mm Hg diastolic 
arterial blood pressure was 103.8+10.7 mm Hg. Patients with symptomatic hypertension were not 
included in the study. Complaints about headache, vertigo, photopsia, sensitivity disturbances 
prevailed! in the clinical picture. of hypertension patients. 2 patients had sustained hypertension. 
Less often cardialgia of different extent was observed. Investigations were carried out along with 
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the patients taking of hypotensive preparations (adelphan, solutio clophelini, cordaphen). Main 
indices of blood and urine of the patients under study did not differ from the age norm. 

The group of gastroenterology (GE) consisted of 6 patients, 4 of them had hepatitis (which 
were complicated with sugar diabetes, gastritis, pancreatitis and nephrite), while 2 had chronic 
cholecystitis. The age of the patients varied from 45 to 64. 

Separate group of mixed pathology was composed of 4 patients who had coronary disease 
complicated by gastroenteric pathology (CD+GE): peptic ulcer, hepatitis, diabetes mellitus. 

The other groups were: 6 patients with acute bronchial pneumonia (P), 6 patients with acute 
respiratory diseases (ARD) and 3 patients with different inflammatory processes (erythema nodular, 
infective allergic polyarthritis) (IP). 

The patients of the P group not only had the clinical signs of the disease (increase in 
temperature, moist cough, dullness and bluntness of percussion sound, the presence of small moist 
bubbling rales), but also the typical laboratory changes: left shift of leukocyte formula, increase in 
the rate of erythrocyte accumulating, increase in the amount of gamma-globulin in blood. 
Analogous changes were observed among the patients of the ARD group. Pathological changes in 
lungs, throat and nose were confirmed by radiography and thermography. ; 

Experiments, conducted on 2 persons from control group (C) in four symmetrically located 
points showed that the fluctuation value of the AF was less than 0.4 dB (10%). 

The main stage of the investigation was conducted as follows. Acupuncture points were 
chosen symmetrically on the both upper limbs of the patients. The frequency range, in which 
investigations were carried out, was 58.8—60.8 GHz, i.e., over 2.0 GHz. There was the same for all 
patients row of 8 peaks of the FRF in this range (Tab.2). 


Table 2 
AF, value number distribution over peak frequencies for different pathology groups (run 1) 
AF, total 
number over Peak frequency, GHz 
freq.range 


The maximum value of the AF,,,, was determined for each spectrogram curve by means of 
the technique described above (Fig.2). This value corresponded to a peak with some frequency 
Е (GHz). For example, in Fig.2, Е = 57.20 GHz, AF, = 2.2 dB. Measured pairs of the values 
AF naxy Е made up a set, which was subjected to further analysis. | 

In the second run 15 patients with disturbances in heart rhythm (ciliary arrhythmia) were 
under study. The patients were divided into two' groups according to etiology of the disease. In 
addition, 5 healthy persons were under study (C). 

The first group of patients (IP) consisted of 8 middle-agers (45—55 years), whose diagnoses 
according to the clinical investigation were: rheumatism, active phase, I-II levels of activity, 
rheumocarditis. All of the 8 patients had mitral coronary deficiency. 2 of them had also stenosis of 
the left atrioventricular orifice. One of the patients had mitral-aortal disease. 5 of the patients had 
polyarthritis. Rheumatotests were positive for all the members of the group. The degree of blood 
circulation deficiency was: the I-st in 3 patients, II-A in 3 patients, what was confirmed by the 
existence of oedemae of shins and feet. 
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The second group (CD) comprised 7 patients of age 45—59. 3 of them had the diagnosis of 
Myocardial infarction (the investigation was conducted in recovery period), which appeared on the 
background of coronary atherosclerosis. 4 of the patients had ciliary arrhythmia which appeared on 
the background of postinfarction and atherosclerotic cardiosclerosis. 

All the patients under study had ciliary arrhythmia of a permanent nature. Pulse deficiency 
was an average 8 + 2 beats a min. in the group. 

FRF spectrograms of the AF value were obtained in acupuncture points EH 7 Daling on the 
both upper limbs of the patients. \ 


RESULTS 


As the results of the experiments MMW spectrograms of the AF values were obtained from 
the upper limbs of the patients in the acupuncture points LI 13 Shouwuli (27 patients) and He 7 
Shenmen (9 patients). The total number of the measured spectrograms (on the right and left hands) 
is presented in Tab.1. The number of the patients and the number of the spectrograms, where the 
AF, value is equal or more than 0.5 dB. 

Comparison between the AF values measured in two acupuncture points showed that the 
data obtained in the He 7 points give a higher value of the disturbance signal/normal signal ratio 
than that of the LI 13 point in the groups CD and CD- GE, whereas measurement in the point 
LI 13 is more preferable for the group GE. 

Tab.2 shows the distribution of the total (for all of the spectrograms) number of the AF, 
values, which pertain to one frequency of the range for different groups of pathologies. One can 
see from Tab.2 that there is a correlation between the group of diseases and the frequency (GHz) 
on which falls the AF, value in the patient's spectrogram. While in the groups CD, GE and 
CD+GE the AF,, values was observed all over the frequency range under study D (58.8— 
60.8 GHz), in the groups P, IP, ARD the AF,,,, values were detected іп low subrange D, only 
(58.8—59.4 GHz), whereas in the upper О. (59.6—60.8 GHz) no value was noted. 

Subsequent measurements of the patients in the P and IP groups after a two day's interval 
showed that during common treatment of the patients with the disease in acute phase the AF, 
values typical for this group had decreased, approaching the norm, which was determined earlier 
for the C group. 

The case of latent pathology revelation should to be mentioned. Patient G. was primarily on 
the healthy group list. But the investigation of the AF showed significant deviations from the norm. 
It turned out thus that G. suffered from chronic antritis. 

The results of the second run: FRF spectrograms of the AF values were obtained in the 
acupuncture points EH 7 Daling on the both upper limbs of the patients. Tab.3 shows the total 
number of AF,,,, values at the peak frequencies for pathology groups studied. Tab.4 and 5 show the 
data of the repeated AF nay measurements with the interval of 5 min. 


Table 3 
АБ, ах value number distribution over peak frequencies for different pathology groups (run 2) 


АБ ах total 


number over Peak frequency, GHz 
freq.range 
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Table 4 
AF nax variation during 5-min MMW irradiation 
a) coronary disease 


Patient's limb No. AF na dB 
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ТаЫе 5 
AF max variation during 5-min MMW irradiation 
b) inflammatory processes 


DISCUSSION 


The experimental data presented here show that using our technique to investigate the 
patients with aforementioned diseases there are some indications, that may help to distinguish 
between the diseases. 

We divided all the patients under study into 3 groups according to Bayes method [9]: 

5; — CD, GE and CD+GE groups; S, — P, IP and ARD groups, S, — the control group 
(healthy persons). We chose the condition of AF, > 0.5 dB as E, symptom. The condition of 
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getting AFma into subrange D; was chosen as E; symptom. Then the probability of the absence of 
this symptom, i.e., the probability of E, symptom was determined. Г 
Using the data given in Tab.1 and 2, we obtained the diagnostic table (Tab.6), where P(S;) is 


a priori probability of the S; diagnosis, KE, ;) is the probability of the Е, indication in case of S; 


diagnosis, P(E, |S j) is the probability of the absence of the E; symptom in case of S; diagnosis. 


Table 6 
Bayes diagnosis chart of symptom probability for pathology groups 


10 
E 1.08 1.0 
0.06 0 


The comparison between the results of AF measurement in the points He 7 and LI 13 
enables us to draw the conclusion that the manifestations of the anisotropic structures on the 
patients’ skin detected by MMW reflection are related to the state of the corresponding 
acupuncture channels. This confirms the idea of an electromagnetic nature of the acupuncture 
diagnosis [6]. 

The AF,, value decrease during the patients’ treatment confirms the objectivity of Фе 
method as well as allows us to suppose its employability for noninvasive monitoring of patients' 
state during the treatment and for predicting its success. 

The case of patient G. showed that the proposed technique may be useful for the diagnosis of: 
latent pathologies. 

The data obtained in the second run showed that the conclusions of the previous work about 
the correlation between the disease groups and the AF, distribution with frequency were- 
confirmed. Namely: if one chose among the spectrograms those in which in the D» subrange 
(59.6—60.8 GHz) the maximum AF values (for the right or the left hand) were present, those 
patients were to belong to the CD group. 

Using the indications mentioned in Tab.3, 8 of 10 CD patients and 4 of 5 IP patients were 
identified correctly. 

Tab.6 provided the probabilities of $, (CD), S; (IP) and S; (C) diagnosis in case of the 


presence of the symptom complex E і.е., both E; and Е, symptom intersection: P(S 8) = 3896; 


P(S, E) = 62%; P(S, E) = 0%. This data make evident that the probability of the second diagnosis 


is higher than the probability of the first one. This result allows us to draw the conclusion at first 
about the frequency dependence of a human organism response to microwave electromagnetic 
radiation and, at second, about the connection between the pathology type and the frequency 
characteristic of the response. 

Thus the technique we suggested not only can help to distinguish the healthy person from the 
sick one, but also to differentiate patients according to the groups of diseases. 
The results of the data analysis on the AF temporal changes showed that during a 5 minutes 
MMW irradiation of the point under study the changes had various features. Permanent ,values 
above the norm were seldom detected. The cases when an increase in the values above the norm 
took place on the both upper limbs of a patient were also rare. Most often a change either in the 
normal range or with the later return to the normal value took place in one of the limbs of a 
patient. An increase in the AF above the normal value took place in the other limb of the same 
patient whose initial value was normal as well as above normal. 

These observed features allow us to draw the conclusion about the asymmetry of the MMW 
effect on the acupuncture points EH 7 of the right and the left limbs of a patient. The tendency of 
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the AF value to the normal limits revealed the normalizing effect of the radiation on the patient, 
i.e., that the irradiated point on the limb was useful for the MMW therapy. 
An increase in the AF. above the normal value or the permanent value above the norm 
showed that the irradiated point was useful for the functional state and latent disease diagnosis. 
Thus the application of the technique suggested for investigation of the patients with heart 
rhythm disturbances showed that the frequency characteristics of the AF,,,, differed depending on 
the pathology group. 


CONCLUSION 


Differential reflectivity DR in MMW range was measured in the acupuncture points for the 
set of patients’ groups and the control group. The table of the AF, frequency distribution for 
different groups was obtained. It was found that the frequency distribution of the AF maximum 
values was correlated with the pathology type of the patients’ groups. It was noted that the AF 
value decreased to the norm during the treatment of the patients with acute. phase of inflammatory 
processes. 

The suggested technique may be helpful for noninvasive diagnostic testing. 

The conducted investigations convincingly showed that the pathological changes of 
generalized nature such as rheumatism or atherosclerosis, which cause the development of coronary 
disease, cause not only heart changes but also affect the biologically active points. It is very 
important for diagnosis, that the patterns of mentioned changes are dissimilar for different 
pathologies. This allows us to use those differences for the disease diagnosis and the choice of an 
adequate treatment. 


ЗАСТОСУВАННЯ РЕФЛЕКТОМЕТРІЇ ШКІРИ ЛЮДИНИ B КВЧ-ДІАПАЗОНІ З МЕТОЮ 
ДІАГНОСТИКИ ДЕЯКИХ ЗАХВОРЮВАНЬ 


ТА ІВАНЧЕНКО, В.Г.ЛИЗОГУБ, Л.В.СВЕШНИКОВА, Ю. В.ЧОВНЮК 


Обстежено деякі групи пацієнтів за запропонованою нами методикою. Вона полягала у 
визначенні діелектричних властивостей людської шкіри в КВЧ-діапазоні електромагнітних 
хвиль шляхом вимірювання диференціального коефіцієнту відбиття хвилі від поверхні 
шкіри в рефлексогенних зонах. Результати дослідження різних груп пацієнтів показали, що 
частотні характеристики розподілу позиції максимуму коефіцієнту анізотропії розрізнялися 
залежно від групи патології захворювань. Таким чином, запропонована нами методика 
дозволяє не тільки відрізняти хворих від здорових осіб, але також диференціювати групи 
захворювань за допомогою КВЧ-тестування. 


ПРИМЕНЕНИЕ РЕФЛЕКТОМЕТРИИ КОЖИ ЧЕЛОВЕКА В КВЧ-ДИАПАЗОНЕ ДЛЯ 
ДИАГНОСТИКИ НЕКОТОРЫХ ЗАБОЛЕВАНИЙ 


И.А.ИВАНЧЕНКО, В.Г.ЛИЗОГУБ, Л.В.СВЕШНИКОВА, Ю.В.ЧОВНЮК 


Обследовались некоторые группы пациентов с помощью предложенной нами методики 
определения локальных диэлектрических свойств кожи человека в КВЧ-диапазоне 
электромагнитных волн путем измерения дифференциального коэффициента отражения 
волны от поверхности кожи в рефлексогенных зонах. Результаты исследования различных 
групп пациентов показали, что, в зависимости от группы патологии заболеваний, частотные 
характеристики распределения позиции максимума коэффициента анизотропии различны. 
Таким образом, предложенная нами методика позволяет не только отличать больных от 
здоровых, но и проводить дифференциацию между группами заболеваний с помощью КВЧ- 
тестирования. 
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Efficacy and Safety of the Microwave Resonance Therapy in Rehabilitation of 
Patients with Cerebral Palsy: a Double-Blind, Comparative Study 


M.I.SKOPYUK, A.A.SOLOVYEVA 


Scientific Research Center "Vidhuk", Kiev, Ukraine 
Center for Medical & Social Rehabilitation of Children with Infantile Cerebral Palsy "Tonus", Kiev, Ukraine. 


Microwave Resonance Therapy (MRT) is a nonmedicinal, noninvasive method of treatment based 
on the local action of millimeter band electromagnetic radiation on selected zones of skin surface. 
The output power of electromagnetic radiation does not exceed 3 mW/cm?. 


A double-blind comparative study of MRT efficacy and safety was carried out during the 
treatment of 51 patients with infantile cerebral palsy (ICP). The main purpose of the study was to 
prove the applicability of electromagnetic radiation with multifrequency output power specum to 
motor disorders rehabilitation in children with ICP. 


After a neurological assessment, a two-week’s treatment was given. All patients were randomly 
assigned to three groups according to the treatment provided. The first group was treated with 
MRT (according to SRC "Vidguk" protocols) in combination with the conventional physical 
therapy. The second group received only physical therapy. The third group received "placebo" 
(the output signal of radiation devices was off in operating range of frequencies) and physical 
therapy. The results of the treatment given to the second and third groups were used for 
comparison. 


The treatment of the first group of patients resulted in improvement in 6096 of the patients while 
in the second and the third groups the similar result was reached only in 3096 of the patients. A 
significant decrease in spasticity was clinically and instrumentally measured in patients of the first 
group what enabled them to improve functional abilities (attainment of a new skill). In addition, 
no negative changes, side-effects or complications were observed in the first group. 


The conclusion is drawn that MRT in complex treatment is a safe and effective method for 
rehabilitation of motor disorders in children with ICP and it can be recommended for broad 
clinical rehabilitation practice with no complications and side-effects. 


( Received September 19, 1994) 


INTRODUCTION 


One. of the most convincing applications of physics of the alive is the establishment of 
quantum medicine currently based on microwave resonance therapy whose wide spread has been 
preceeded by numerous investigations into millimeter field action on biosystems of various levels of 
hierarchy of the alive up to human [1]. In particular, it is proved that MRT is an effective method 
for the treatment of motor disorders in children with ICP and does not cause any side effects and 
aggravations. 

According to literature data, there are three types of responses in the functional state 
dynamics: (i) those being formed immediately during the exposure, for example, sensory response 
[2]; (ii) those formed as a result of a single exposure, e.g., the change in the immune system [3] 
and (iii) those developing after a prolonged course of treatment, e.g., healing of ulcer [4]. The 
electromagnetic exposure was also proved to affect tonus state of patients with ICP [5]. 

According to the modern concepts, ICP is understood as a symptomatic complex of motor 
and tonus disorders [6] which are caused by the cerebral ontogenesis disturbances of the brain 
during antenatal and early postnatal periods. 

Children with the most common ICP form, that is spastic diplegia (SD) were enrolled in the 
study [7]. To control MRT efficacy and safety during the treatment, we used neurological indices 
which characterized patients' motions. 
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PATIENTS AND METHODS 


Patients screening 


The trial was organized as a randomized double-blind study in parallel groups. The study was 
set on in September, 1991. The last patient was enrolled in the study in May, 1992. The study 
ended before June, 1992. 

The groups included children aged 4—15 years with ICP in the spastic diplegia form 
(tetraparesis affecting primarily the lower extremities) in final residual state. The form of disease 
was diagnosed according to the effective classification of Ministry of Health approved in 1973 [8]. 

The patients lived in the city of Kiev, Zhytomyr, and Chernihiv regions of Ukraine in a 
comparatively similar environment and had similar nutrition patterns. 

After the treatment the children were followed up for nine months to observe long term 
clinical efficacy and side effects. 

Children: (i) with light form of spastic diplegia (having only cosmetic defects); (ii) with other 
ICP forms; (iii) with accompanying mental disability (moderate or severe mental retardation; (iv) 
with remoteness of the last fit of epilepsy less than three years ago; (v) at the age under four and 
over fifteen; (vi) living in other areas than those mentioned above were not included in the study. 

Patients participated in the study on voluntary basis. Every patient or patient's tutor gave 
his/her written consent to participate in the study, record all the necessary information and carry 
out the planned investigation. 


Study design, tasks and purposes of the study 


The present study was carried out to prove the acceptability of electromagnetic radiation with 
multifrequency power spectrum for physical rehabilitation of the children suffering from ICP. 

According to the purpose of the study the following problems were solved: 

1. Study of clinical efficiency and safety of MRT. 

2. MRT efficiency assessment as compared to the results of the conventional conservative 
treatment. 

To solve these problems the patients were assigned to groups in the consecutive order of their 
case histories numbers. 

The treatment was carried out as follows: (i) in the first group the МЕТ and the 
conventional physical therapy were used; (ii) in the second group only the physical therapy was 
used; (iii) in the third group “placebo” as well as physical therapy was used. 

One day before the beginning of the treatment the patients underwent a clinical examination, 
neurological tests were made and an angulation measurement (range of motion) was carried out for 
each of them. The treatment schemes were selected. 

During two five-day periods of the active treatment, with two days of rest between them the 
patients were given the treatment. At the end of the treatment all patients underwent clinical and 
laboratory re-examinations. 

MRT was applied using the electromagnetic radiation generator *POROG". The treatment 
was administered fromg9 to 13 o'clock, its duration being 20—25 minutes with the obligatory rest 
for 30 min. after treatment. 

The number of exposure zones and their compoundings were individually chosen proceeding 
from the patient's neurological state. 

Methods the conventional conservative treatment included physical exercises and, in case of 
necessity, a thermal treatment (paraffin applications to the upper and lower extremities joints). The 

drug therapy during the trial was not given and its unrestricted use was controlled. 
з During “placebo” procedure output radiation of our devices was not in operating 
frequencies range. 
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Patients' compliance 


Regimen and the treatment scheme were controlled by records made in case sheets. The 
patients who- missed more than one day of the treatment were excluded from the subsequent 
analysis of the treatment results. 


Clinical and laboratory tests 


Clinical examinations included anamnesis and ICP history, assessment of the patients' 
general condition (sleep, appetite, etc.), and a neurological examination with determination of 
motor disorder indices characterizing child's statics and locomotion. The posture, carriage, 
ambulation (walking), limb contractures and hands manipulation function were assessed in marks 
[9]. Apart from the clinical examination, an instrumental assessment of active and passive 
movements amplitude of the lower extremities was made by means of angulation measurements. 
The assessment was also made in marks. The system of marks and its correspondence to clinical 
state of the patients is given in Appendix I. 

Evaluation of the treatment efficacy was done as follows: 

— “unsatisfactory” — no improvement in clinical symptoms or aggravation was found; 

— "satisfactory" — improvement in clinical symptoms according to the pre-treatment 
evaluations; 

— "improvement" — improvement in clinical symptoms, enabling the patient to attain new 
skills. ` 
Statistical analysis 


The clinical and laboratory examination data were assessed by the method of statistical 
variation. Treatment results were compared using Student's test and y’-test. For every index its 
distribution was determined, mean deviation and variation limits were found. 


Initial data 


There was no significant difference between the three groups of patients by age, sex, severity 
of the disease and indices of motor disorders and lower extremities range of motion (0.5 « P « 1.0) 
(see Table 1). 


Patients cohort 


Fifty one patient, 24 boys and 27 girls, aged 4—15 years, with SD in final residual stage were 
enrolled in the trial. 

All of them were treated at SRC “Vidhuk” medical department in 1991-1992. Two patients 
(one from the first group and one from the second group) were excluded from the study (acute 
respiratory disease during treatment). 

Thus, 49 patients (23 male and 26 female) were included in the final efficacy analysis. 


Patients characteristics 


According to the severity of the motion defect, the patients were conventionally assigned to 
two subgroups: those with (a) severe and (b) moderate motor disorders (30 and 19 patients, 
respectively). Average age and sex distributions are shown in Table 1. 

The patients of (a) group had problems with sitting and were unable to stand and walk. 

High muscular tonus of the patients hampered their voluntary motion. Active motions were 
made along with general rigidity, slowly and with small amplitude, accompanied by overall and 
coordination synkinesis. 

21 patients had combined bending positioning in hip, knee and ankle joints (sole bending). 
Abnormal positioning were involuntary enforced in the patients' vertical position. 10 patients had 
muscles of neck and trunk included in bending synergias. 9 patients demonstrated straightening 

, Positioning in knee and ankle joints. 


LI 
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Table 1 


Groups of patients participating in MRT safety and efficacy treatment of infantile cerebral palsy 
(principal indices prior to the treatment) 


Control Placebo 
Sex 
р 
на рте БРА Bee СИ, 


Clinical severity of the disease А 


рва uL sab exci d 
mu о ee I odo 8 lu s от 


State of motor sphere in marks MO*/(limits) 
23/1 
19/0-3 
0.6/(0-2 


N 
бл 
T 
ES 
t2 
a 
— 


> 


sitting 
standing 
walking 


r2 
CA 
— 
> 


+ 


left 3.1/(1-5 2.7/0—5 


Range of motions in joints in marks MO*/(limits) 


hip joint 

left 3.8/(2—5 
knee joint 

left 3.9/(3—5 4.0/(2—5 3.6/(2—5 


| 


foot joint 
2.6/(2—4 3.3/(2—5 Г 3.3/(2—5 


e: 
e 


2.5/(2—4 3.3/(2-5 3.3/(2—5 


* — Average meaning of indices in the sample 


The patients’ carriage feature was kyphosis or kypho-scoliosis in thoracic spine. 13 patients 
had the total kyphosis spreading over the lumbar spine. 

All the patients had a selective increase in muscular tonus of adductor groups, resulting in 
intrarotating and adducting hips positioning. Due to acute spasm of hip adductors, bending and 
pronation patients kept their feet crossed. 

15 patients had equinovarus. 

Most of the patients demonstrated major hand function disorders preventing them from self- 
care (0—2 marks). 18 patients had stable flexion radial pronating contractures in radial and 
radiocarpal joints, with pronation at the angle of 90 degrees and more. The patients’ hands and 
fingers were permanently half-bent, the first finger adducted to the palm, the hand is put apart in 
elbow. 

Passive correction of abnormal positions had failed or was minimal and extremely unstable 
one. 
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_ Five patients had a distinctly limited range of hand motions, hampered by pronation 
positioning of arms and hands, flexion positioning of fingers. Passive correction was possible only 
to the medium position. Supination was less than 10 degrees. Children took things awkwardly, witk 
pronatory positioning of the arms and hands. Six patients had a limited adduction of the first fingei 
and could not sufficiently put aside the third and fifth fingers what impeded fine movements 
requiring differentiated work of fingers. The speed, adroitness, hand movements were disturbed. It 
was extremely painful for the children to lace shoes, button themselves up as well as to draw and 
write. 

Considerable alterations were noticed in reflexes, namely enhancement of tendon and 
periostal reflexes with extended areas of their induction in 28 patients. Twenty one patients had 
abnormal reflexes. 

Excessive athetoid and chorea-athetoid movements took place in 8 patients. Hyperkineses 
appeared in case of making an attempt to do an active movement and were not seen at rest. 

Coordination disorders in nose-finger test and dysmetria were seeh in 7 patients. 

Patients of (b) group were sitting confidently, but their ambulation was to some degree 
disordered in every patient. Seven of them had additional swingings in sagittal or frontal plane. 
Step movements were slow, the phase of foot carrying was shorter than usual, the step itself was 
shorter, the rest was predominantly on the fore portions of feet. The fore push and rolling were 
severely distorted. Eleven patients limped, their resting legs was straightened, slightly adducted and 
turned inside. The second foot was bent in hip and knee joints, hip turned inside, the rest was on 
the fore edge of foot. 

Abnormal bending positioining of lower extremities was seen in 11 patients, extension 
positioning of knee and foot joints — in 8 patients. 

Actually every patient had spasms in adducting groups of hip muscles (from light to 
moderate) which did not result in feet crossing. 

Six patients had equinovarus. 

3 patients performed hand manipulation function. The rest of the patients had distorted hand 
functions due to expressed and stable flexor positioning, hampered active foot and arm supination, 
disordered resting and grasping functions. Eight patients had fine movements disorders along with 
insufficiently differentiated movements of the third-fifth fingers as well as limited drawing of the . 
first finger. 

Tendon. and periostal reflexes were enhanced in 10 children with moderate severity of the 
disease. They were combined with abnormal reflexes in 6 children. 

Four patients had a hyperkinetic syndrome in the form of choreo-athetoid hyperkinesis 
appearing in case of active movements. 

Atactic disorders in the form of intention tremor during finger-to-nose test and dysmetria 
took place in 4 patients. | 

According to instrumental findings, relating to the lower extremities position and mobility by 
means of angulometry, it was found in all of the selected patients, that passive movements were 
limited in hip joints of 3 patients, in knee joints of 3 patients, and in foot joints of 16 patients. 
Active movements were restricted or impossible for the majority of the patients. Limitations of both 
active and passive movements in foot joints of 47 patients were especially evident. 


Statistic processing of the investigation results 
The results of MRT efficacy and safety for all groups are shown in Table 2. 
The results of the Ist group treatment 


According to the results of the first group treatment, as a whole, an improvement was gained 
ћи patients who demonstrated positive changes іп most of the neurologic status indices. 
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Table 2 


MRT efficacy and safety investigation results in treatment patients with infantile cerebral palsy 
(statistical analysis of treatment results) 


MEME —— 
Number of patients 


28/0-4 
24/0-4 

1.8/(0—4 
gai 0-2 


А gen 
unsatisfactor 
6. 7 
n 2 


* — Average meaning of indices in the sample 


Sitting was improved in 5 patients. One of the patients developed an ability to sit with 
support, one began to sit without support, two were sitting more confidently and one patient 
reduced kyphotic deformation of the back. 

Two patients developed the ability to stand and walk, walking improvement was found in 7 
patients. Equinovarus was eliminated in 3 patients and feet crossing became less observable in one 
of the three patients. 

Improvement in hands movement was noticed. Manipulation function of hands was 
improved in 4 patients, one of them significantly reduced pronation and flexion pathologic 
positioning of arms and the possibility of their partial both passive and active improvement. Three 
patients developed an improvement in fine motility of hand movements. One patients of the four 
showed a reduction in athetosic hyperkineses in his hand when he tried to do active movements. 

No changes in tendon and periostal reflexes were found. 

Coordination disorders were reduced in the two of the six patients. 

Active movements volume in knee joints were extended in 4 patients. Eight patients gained a 
full range of passive movements in feet joints and one patient demonstrated a considerable increase 
in active movements of both feet joints. No considerable changes in the range of hip joint motions 
were found. No clinical aggravation in the patients group was found. 


мо 
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The results of the second group treatment 


As a treatment result, improvement was recorded in five patients. 

Two patients were sitting better, one of them began to sit without support, and the other 
became more confident in sitting. 

Standing function did not improve. Two patients improved their walking. One of them began 
to move using one hand or walking sticks as a support, while the other's walking became more 
confident. The gait did not change in any of the patients. 

Hand movements improved in 1 patient at the expense of pronatory hand positioning 
reduction. 

No other changes in neurological indices under study were found. 

One patient showed an insignificant increase in active and passive movements range in knee 
joint and two patients improved motions in feet joints. 

Three patients did not demonstrate any changes after the treatment. 

One patient had an increase in equinovarus and enforced athetoid face muscles movements 
in cases of emotional tension. No other worsening of the patients' clinical condition was found. 


The results of the third group treatment 


As a treatment result, two patients demonstrated improvement. 

One patient improved his sitting — he began to sit more steadily 

Standing, walking and gait did not change. 

One patient became able to grasp things with one of his hands. 

No changes in the other neurological indices under study were found. 

One patient showed an insignificant increase in the range of active and passive motions in 
knee joints. 

No clinical effect was observed in two patients. 

No aggravation of the patients' clinical state was noticed. 


Safety control 


Control of operating condition of the used generators was done daily, no breakdowns were 
recorded. The treatment was carried out under the constant surveillance of a physician in charge. 
Every day patients were asked to report any changes in their physical cqniition, any unfavourable 
changes in their neurological status. 


Side effects 


Patients of the Ist and 3rd groups did not report any adverse reactions and aggravation of 
their state. 

One patient from the 2nd group reported a worsening. However, she was ruled out neither of 
the trial nor of the sample, so far as the similar aggravations, according to her parents’ report, were 
observed earlier during treatment at other hospitals. 


DISCUSSION © 


When the treatment results with use of different methods were compared, it had been found 
out that the treatment efficacy greatly depended on the treatment used. 

When locomotive sphere indiges prior and after treatment were compared in every group, 
their analysis showed a Statistically reliable improvement in knee joints (0.1 > P > 0.05) and feet 
joints (P < 0.01) only in the first group wheré MRT was given in combination with conservative 
therapy (Table 3). 

Higher values of principal indices of static and motor functions were mainly seen in the first 
group as compared to groups 2 and 3. 
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Table 3 
State of patients with infantile cerebral palsy prior to and after microwave resonance therapy 


(comparative analysis) 
Group I (MRT) 
State of the locomotion sphere in marks 
2.5/(04) 2.9/ (14) 


3.0/(1-4) 


Range of movements in joints in marks 


hip joint 
left : 

knee joint : 
right 3.8/(3—5) 


Thus, 4 patients of the first group improved manipulatory hand function, while one patient 
only in both the second and third groups demonstrated such an improvement. Five patients of the 
Ist group improved their sitting function, versus two and one patient in the second and third 
groups, respectively; two patients of the first group attained an ability to sit and walk independent- 
ly, while in other patients these functions actually did not change. It should be noted that the 
patients of group 1 improved their walking, going more confidently, while the other groups 
patients' function remained on the same level. 

A general assessment of the treatment results with application of integral "treatment efficacy 
criteria" supports the above particular tendencies. Thus, in the group of patients treated with MRT, 
improvement was recorded more frequently than in those treated with conventional therapy and 
*placebo". Besides, there were no negative treatment results, side or adverse reactions. Nobody of 
the patients treated with MRT complained of general condition or neurological status aggravation. 


CONCLUSIONS 


The investigation results show that the employment of MRT in complex treatment of ICP 
patients, ensures positive changes of principal indices of neurological status characterized by statics 
and locomotion. At the same time MRT is a safe means in treatment of motion disorders in ICP 
patients, which does not produce side and adverse réactions and complications. 
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EFFICACY AND S 


Appendix 
System of the disease severity assessment in ICP patients 


зро ee UI cm Trl NACER E 


A 


Function 


Sitting Patient is unable Sits with support Sits independently Sits steadily, but with a | Sits without 
to sit irrespective but unstably, with a | carriage defect expressed deviations 
of support carriage defect from the norm 
Standing Patient cannot Sits in an abnormal | Stands near a Stands independently Sits without 
maintain standing | way with a support with preservation of an expressed deviations 
position, considerable with preservation of | abnormal posture from the norm 
irrespective of support an abnormal 
support posture 
Walking Pace movements Moves with a Moves with a Moves, but very shaky Walks independently, 
are not formed, considerable support on one with a moderate abnormal positioning 
cannot walk even bilateral support hand abnormal positioining of | o lower extremities is 
with a support lower extremities expressed 
inconsiderabl 
Gait Not formed Unsteady, non- Unsteady, non- Steady, with pace Without expressed 
rhythmic, with rhythmic, with pace | phases violations deviations from the 
additional swings in | phases violations norm 
saggital and frontal 
planes 
Abnormal Persistent bending | Passive correction Physiologically An extremity can be 


of an abnormal 
positioning is 
possible, active one 
is limited 


contracture 
permitting partial 
passive correction; 
active correction is 
impossible 


positioining 
of lower 
extremities 


normal positioning 
of an extremity is 
possible but it is 
not maintained 
under a patient's 
body weight 
Grasps things with 
difficulty, with 
perverted 
positioning of 
upper extremities, 
passive correction is 
impossible or little. 
Active movements 


placed in a position 
close to the norm 


Hand 
manipulation 
function 


Hands 
functions 
approach the 
norm 


Fine hands 
movements area 
distorted, motor 
clumsiness 


Grasps things with 
other person's help, 
but. does not hold 

them 


Grabs things in a 
clumsy way with 
pronatory hands 
positioning; partial 
active correction of 
abnormal hands 
positioning is possible 
Active movements - less 


Active and No active or No active or passive Volume of active Volume of 


passive lower. | passive movements | movements in - less than 25%, than 50% of their movements 50 to active 
extremities joints, passive joints | passive movements | physiological volume, 75% of physiological movements 50 
movement are possible less - 50 to 100 % of full volume of passive volume, passive ones | to 75% of 
range than 50% their physiological movements are in full volume physiological 
range volume, 
passive ones 
are in full 


volume 
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ЕФЕКТИВНІСТЬ I БЕЗПЕКА МІКРОХВИЛЬОВОЇ РЕЗОНАНСНОЇ ТЕРАПІЇ B 
ЛІКУВАННІ ДИТЯЧОГО ЦЕРЕБРАЛЬНОГО ПАРАЛІЧА: ПОДВІЙНЕ СЛІПЕ 
ПЕРЕХРЕСНЕ ДОСЛІДЖЕННЯ 


M.LCKOIIIOK, A.A.COJIOBHOBA 

Мікрохвильова резонансна терапія (МРТ) - безмедикаментозний неінвазійний метод 
лікування, який здійснюється шляхом локальної дії електромагнітного випромінювання 
міліметрового діапазону (КВЧ) на виділені зони поверхні шкіри. Потужність електромаг- 
нітного випромінювання не перевищує 3 мВт/см. 

Подвійне сліпе порівняльне дослідження ефективності МРТ провадилося під час лікування 
51 пацієнта, хворого на дитячий церебральний параліч (ДЦП). Метою дослідження було 
визначення придатності джерел електромагнітного випромінювання з полічастотною спек- 
тральною характеристикою вихідного сигналу для фізичної реабілітації хворих дітей. 

Після оцінки неврологічного статусу пацієнтам призначалося двотижневе лікування. В 
залежності від проведеного лікування їх було розподілено на три однорідні за складом та 
тяжкістю захворювання групи. Першу групу лікували за загальноприйнятою методикою 
фізіотерапевтичного лікування в поєднанні з мікрохвильовою резонансною терапією 
(методика MHIL, “Відгук”). Хворі другої групи дістали лише фізіотерапевтичне лікування. У 
лікуванні третьої групи фізіотерапія доповнювалась процедурою "плацебо" (під час лікуван- 
ня був відсутній вихідний сигнал джерела електромагнітних випромінюваань у робочому 
діапазоні частот). 

У результаті лікування в першій групі досягнуто "поліпшення" у 6096 хворих, тоді як у 
групах 2 і 3 одержано аналогічний результат лише у 3096 хворих. Клінічно та інструменталь- 
но зареєстровано істотне зниження спастики у хворих першої групи, що дозволило розши- 
рити їхні можливості у розробці рухових навичок. 

Доведено, шо МРТ, ефективний метод лікування порушень рухових функцій у дітей з ДЦП, 
не викликає побічних ефектів і ускладнень. 


ЭФФЕКТИВНОСТЬ И БЕЗОПАСНОСТЬ МИКРОВОЛНОВОЙ РЕЗОНАНСНОЙ 
ТЕРАПИИ В ЛЕЧЕНИИ ДЕТСКОГО ЦЕРЕБРАЛЬНОГО ПАРАЛИЧА: ДВОЙНОЕ. 
СЛЕПОЕ ПЕРЕКРЕСТНОЕ ИССЛЕДОВАНИЕ 


М.И.СКОПЮК, А.А.СОЛОВЬЕВА 


Микроволновая резонансная терапия (МРТ) — безмедикаментозный неинвазивный метод 
лечения, который реализуется путем локального воздействия электромагнитным излучени- 
ем миллиметрового диапазона (КВЧ) на выделенные зоны кожной поверхности. Мощность 
злектромагнитного излучения не превышает 3 мВт/см. 

Двойное слепое перекрестное исследование эффективности МРТ было проведено при 
лечении 51 пациента, болеющего детским церебральным параличем (ДЦП). Целью исследо- 
вания было определение пригодности источников электромагнитных излучений с поличас- 
тотной спектральной характеристикой выходного сигнала для физической реабилитации 
больных детей. і 

После оценки неврологического статуса пациентам назначалось двухнедельное лечение. B 
зависимости от проводимого лечения их распределили на три однородные по составу и 
тяжести заболевания группы. Первую группу лечили по общепринятой методике 
физиотерапевтического лечения в сочетании с микроволновой резонансной терапией 
(методика МНИЦ “Видгук”). Больные второй группы получали только физиотерапевтичес- 
кое лечение. При лечении третьей группы физиотерапия дополнялась процедурой 
“плацебо” (при лечении отсутствовал выходной сигнал источника электромагнитных 
излучений в рабочем диапазоне частот). 
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В результате лечения в первой группе достигнуто “улучшение” у 60% больных, тогда как в 
группах 2 и 3 получен аналогичный результат только у 30% больных. Клинически и инстру- 
ментально зарегистрировано существенное снижение спастики у больных первой группы, 
что позволило расширить их возможности в выработке новых двигательных навыков. 
Доказано, что МРТ — эффективный метод лечения нарушений двигательных функций у 
детей с ДЦП, не вызывает побочных эффектов и осложнений. 
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EVENTS 


On June 20-24th, 1994 the Ist congress of the Ukrainian Biophysical Society was held at the 
premises of Department of Biology, Kiev Shevchenko University. Over 200 scientists from Ukraine 
representing scientific centers and institutions participated. Scientists from Russia, Canada and 
Germany took part in the work of the congress. 

The keynote address *Conditions and development prospects of biophysics in Ukraine" was 
given by Prof. M.F.Shuba, academician of NAS of Ukraine. During the first plenary session 
lectures of academician P.G.Kostyuk “Structural апа functional organization of a calcium 
channel”, of academician D.M.Grodzinsky “Contemporary biophysical problems of ionizing 
radiation action on biological systems", of Prof. S.P.Sit'ko "Biophysics and quantum physics of the 
alive” were delivered. 

Work of the congress was split in the following sections: biophysics of nucleic acids, 
biophysics of proteins, biophysics of sensory cells, membrane transport of ions, biophysics of 
mitochondria, biophysics of muscle, photobiology, radiation biophysics, peculiarities of bioobjects 
reactions to the action of physical factors, medical biophysics, theoretical biophysics, biophysical 
education. The plenary and sessional reports given as well as the numerous poster presentations 
presented had shown a high scientifical level of the conducted research. The congress results 
witness that in Ukraine fundamental scientifical research is currently conducted in all directions of 
contemporary biophysics — molecular and cellular biophysics, biophysics of complex systems and 
applied biophysics. 

It had been planned to hold the 2nd congress of the Ukrainian Biophysical Society in 
Kharkov in 1997. 


SS 
NU 


PHYSICS OF THE ALIVE 


Synergetics and Phase Transitions: Mounting the Quantum Ladder of Nature — 
A.V. Chaly, I.S. Dobronravova, S.P.Sit'ko 


Corrective Action of Millimeter Wave upon Biological Systems of Different 
Hierarhial Levels — L.S. Bundyuk, А.Р. Kuzmenko, N.N. Ryabchenko, G.S. Litvinov, 
S.A. Borovnik, L.I.Matseiko, І. У. Gorbenko, О. К. Fedorovich 


Effects of Ion- Conformational Interaction in Biomembrane Channels — 
С.Е. Weinreb 


Nonlinear Effect of Dynamical Self-Organization in Macromolecular Systems 
Which Is Induced by Photocontrolled Electron Flux — A.O.Gushcha, 
A.A. Dobrovolsky, М.Т. Kapustina, A.V. Privalko, V. №. Kharkyanen 


BIOPHYSICS 


Calcium Signal: Registration in Cells and Quantitative Estimation — V.L.Zyma, 
О.М. Dyachok 


Biological Mechanoluminescence (Review article) — V.E.Orel 


Deterministic Chaos in Cytosolic Calcium Oscillations — Р.Е. Strizhak, 
I.S. Magura, К.В. Yatsimirskii 


MEDICAL AND BIOLOGICAL APPLICATIONS 


Human Skin Millimeter Wave Range Reflectometry Application for Diagnosis of 
some Diseases — Г.А. Ivanchenko, V.G.Lizogub, L.V.Sveshnikova, Yu.V. Chovniuk 


Efficacy and Safety of the Microwave Resonance Therapy in Rehabilitation of 
Patients with Cerebral Palsy: a Double-Blind, Comparative Study -- M.I.Skopyuk, 
А.А. Solovyeva 


EVENTS 


The 1st congress of the Ukrainian Biophysical Society 


